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and  each  set  corresponded  to  a  different  configuration  of 
the  test  section.  The  test  section  configuration  of  each 
of  the  four  sets  of  experiments  were:  free  annular 
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determined  the  behavior  of  the  mean  flow  profiles  at 
several  axial  locations  for  the  four  sets  of  tests  with 
the  most  extensive  investigation  performed  in  the  fourth 
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set  of  tests.  Five-hole  probes  were  used  to  determine  the 
flow  property  variations,  and  their  development  and 


calibration  are  described  in  detail.  The  variation  of  the 


axial  velocity,  tangential  velocity,  total  pressure,  and 


static  pressure  profiles  and  their  behavior  are 


presented.  The  measurements  reveal  the  highly  complex 


flowfield  and  the  influence  of  the  swirl  on  the  mixing 


behavior.  A  short  analytical  model  of  an  ideal  mixer  with 
swirl  is  presented  to  estimate  the  effects  of  swirl  on 


"ideal”  mixing. 
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CHAPTER  1 


INTRODUCTION 

The  turbulent  mixing  of  two  concentric  incompressible 
streams  with  swirl  present  in  the  inner  stream  has 
received  an  increasing  amount  of  attention  during  recent 
years.  The  primary  motivations  for  this  have  been  the 
application  to  combustion  and  the  swirling  wake  behind  a 
propller  for  propulsion  or  a  windmill  (Reference  1) . 
Swirling  flows  have  been  classified  by  King,  Rothfus,  and 
Kermode  (Reference  2)  according  to  the  wall  boundary 
conditions  in  the  following  fashion: 

1.  unconfined  swirling  flows  -  where  the  wall 
effects  are  negligible. 

2.  Small  L/D,  confined  swirling  flows  -  in  short, 
large-diameter  chambers  where  sidewall  effects  strongly 
interact  with  the  swirl  to  produce  significant  secondary 
flows. 

3.  Large  L/D,  confined  swirling  flows  -  in  tubes 
where  circumferential  wall  effects  interact  strongly  with 
the  swirl  flow. 

To  date,  most  investigations  of  confined  mixing  with 
swirl  present  have  considered  the  mixing  of  fuels  and  air 
in  cylindrical  combustors  where  a  secondary  reverse  flow 


provides  the  necessary  flame  stability  and  greatly 
influences  the  mixing  process  (References  3,  4,  5,  6  and 
7).  Turbulent  mixing  with  swirl  has  many  engineering 
applications  where  a  centerbody  may  be  present  and  the 
secondary  reverse  flow  is  absent.  Such  an  application  is 
present  in  modern  turbofan  engines  where  mixing  of  the 
engine  core  stream  and  the  fan  stream  is  utilized  to 
improve  the  engine  thrust  and/or  to  reduce  the  nozzle 
exhaust  noise. 

Experimental  data  and  analytical  methods  are  available 
for  unconfined  swirling  flows  and  small  L/D  swirling 
flows  in  combustion  chambers  (References  l  through  8) .  A 
limited  amount  of  experimental  data  is  available  for  the 
turbulent  mixing  of  two  concentric  incompressible  streams 
in  a  constant  cross-sectional  area  channel,  with  swirl 
present  in  the  inner  stream.  The  only  work  in  this  area 
known  to  the  author  was  done  by  Launder  and  Morse 
(Reference  9)  and  was  concerned  with  a  swirling  jet  with 
and  without  an  external  coaxial  stream.  No  known  work  has 
been  done  with  a  centerbody  present. 

The  focal  point  of  this  study  was  experimentally 
obtaining  the  variation  of  mean  flow  properties  for  the 
turbulent  mixing  of  two  concentric  Incompressible  air 


streams,  with  swirl  present  in  the  inner  stream,  flowing 


through  a  constant  area  annulus,  in  addition, 
experimental  data  was  obtained  for  the  turbulent  mixing 
without  an  outside  wall  and  with/without  a  centerbody. 
Data  was  also  obtained  for  turbulent  mixing  of  a  swirling 
annular  jet.  Four  sets  of  experiments  were  performed  with 
a  swirl  ratio  (tangential  to  axial  velocity  ratio)  of  0.8 
in  the  inner  stream.  The  test  configuration  for  each  test 
is  listed  in  Table  1-1.  The  most  extensive  experiment  was 
performed  for  test  #4,  the  focal  point  of  this  study. 


TABLE  1-1 
EXPERIMENTS 


Test 

Flow  present 
in  Outer 
Annulus 

Center- 
Body  in 
Place 

Outer 
Wall  in 
Place 

a 

(m0/m,  ) 

1 

NO 

No 

No 

0.00 

2 

yes 

No 

No 

1.00 

3 

yes 

Yes 

No 

1.00 

4 

yes 

yes 

yes 

0.00 

4 

yes 

yes 

yes 

0.47 

4 

yes 

yes 

yes 

1.00 

4 

yes 

yes 

yes 

2.13 

4 

yes 

yes 

yes 

3.91 

where  m,  is  the  inner  stream's  mass  flow  rate 
and  mo  is  the  outer  stream's  mass  flow  rate 


A  five-hole  probe  was  used  for  mean  flow  measurements. 
The  experimental  measurements,  test  apparatus,  five-hole 
orobe.  data  reduction,  and  discussion  of  results  are 


k 


contained  in  Chapters  3,  4,  5,  6  and  7,  respectively. 

A  one-dimensional  analytical  model  of  a  constant  area 
ideal  incompressible  mixer  with  swirl  was  developed  and 
compared  to  a  mixer  without  swirl.  This  analytical  work 
is  presented  in  Chapter  2. 
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CHAPTER  2 


ANALYSIS  OF  AN  IDEAL  SWIRL  MIXER 


I.  INTRODUCTION 


Turbulent  mixing  of  two  concentric  incompressible  streams 
with  swirl  present  in  the  inner  stream  is  of  interest  in 
several  engineering  applications.  Research  is  presently 
being  conducted  to  obtain  a  better  understanding  of  the 
mixing  process.  To  date,  most  investigations  of  mixing 
with  swirl  present  have  considered  the  mixing  of  fuels 
and  air  in  combustors  where  a  secondary  reverse  flow 
provides  the  necessary  flame  stability  and  greatly 
influences  the  mixing  process  (References  l  and  2) . 
Turbulent  mixing  with  swirl  has  many  engineering 
applications  where  a  centerbody  may  be  present  and  the 
secondary  reverse  flow  is  absent,  such  an  application  is 


present 

in 

modern 

turbofan 

engines  where  mixing  of  the 

engine 

core 

stream 

and  the 

fan 

stream  is  utilized  to 

improve 

the 

engine 

thrust 

and/or 

to  reduce  the  nozzle 

exhaust  noise,  present  day  turbofan  engines  remove  most 
of  the  swirl  from  the  engine  core  stream  with  the  turbine 
exit  nozzle  row  and  then  mix  the  two  concentric  parallel 
jets  in  a  mixer  which  contains  devices  (e.g.  paddle 
mixers)  to  enhance  the  mixing  process.  These  forced 
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mixing  devices  also  reduce  the  total  pressure  of  the 
mixed  stream  and  thus  reduce  the  available  thrust  from 
the  turbofan  engine.  It  is  possible  that  the  mixing 
losses  of  the  engine  core  stream  and  the  fan  stream  can 
be  reduced  by  leaving  off  the  turbine  exit  nozzle  row  so 
that  the  engine  core  stream  enters  the  mixer  with  swirl. 
The  swirl  of  the  combined  stream  would  subsequently  be 
removed  by  guide  vanes  at  the  exit  to  the  mixer.  The 
presence  of  swirl  in  the  inner  stream  entering  the  mixer 
makes  the  boundary  between  the  inner  and  outer  streams 
unstable  because  the  fluid  angular  momentum  decreases  in 
the  radial  direction  (Reference  3) .  This  instability  in 
the  radial  direction  enhances  the  turbulent  mixing  of  the 
two  streams. 

in  the  following  section,  an  ideal  incompressible 
constant  area  mixer  with  swirl  is  modeled  and  its 
performance  determined.  The  swirl  velocity  profile  of  the 
inner  stream  entering  the  mixer  is  assumed  to  be  that  of 
a  free  vortex.  The  swirl  velocity  profile  of  the  mixed 
stream  is  assumed  to  be  that  of  a' solid  body  rotation. 
Both  streams  entering  the  mixer  and  the  departing  stream 
are  assumed  to  have  uniform  axial  velocity. 


The  performance  of  an  ideal  incompressible  constant  area 


mixer  without  swirl  is  also  determined  so  that  the 
effects  of  swirl  can  be  determined  by  comparison  of  the 
two  solutions. 

II.  CONSTANT  AREA  MIXER  WITH  SWIRL 

In  this  section,  the  ideal  incompressible  constant  area 
mixer  with  swirl  is  considered.  The  performance  equations 
are  developed  in  Appendix  B.  The  configuration  of  this 
ideal  mixer  with  swirl  is  shown  in  Figure  2-1.  The 
following  conditions  are  prescribed: 

a.  An  incompressible  stream  enters  the  annular  channel 

"0",  formed  between  the  outer  radius  of  rt  and  the  inner 

radius  of  rh  ,  with  uniform  density  (  pQ)  ,  uniform  total 

pressure  (p_  ),  uniform  axil  velocity  (w.)  and  no  radial 

To  0 

or  swirl  velocity. 

b.  An  incompressible  stream  enters  the  annular  channel 

"2",  formed  between  the  outer  radius  of  Rt  and  the  inner 

radius  of  rt  with  uniform  density  (p2),  uniform  total 

pressure  <p_  ),  uniform  axial  velocity  (w.)  and  no  radial 
t2 

or  swirl  velocity. 

c.  The  incompressible  stream  in  annular  channel  "0" 
passes  through  stator-i,  an  ideal  free  vortex  stator,  and 
departs  with  a  free  vortex  swirl  velocity  profile  and  no 
radial  velocity  at  station  "1".  The  total  pressure  of  the 
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inner  stream,  p_  ,  and  the  cross-sectional  area  are 

o 

unchanged  through  stator-1. 

d.  The  ratio  o£  the  static  pressure  of  the  inner 
stream  with  swirl  at  its  outer  wall,  pd ,  to  the  static 
pressure  of  the  outer  stream  without  swirl,  P2 ,  meets  a 
modified  Kutta  condition  at  the  end  of  the  cylindrical 
wall  of  radius  rt  dividing  the  two  incompressible 
streams.  This  modified  Kutta  condition  is  specified  by 
the  ratio  Pd/Pr 

e.  The  inner  stream  with  swirl,  "1",  and  the  outer 
stream  without  swirl,  "2",  enter  a  constant  area  mixing 
region  in  which  the  two  incompressible  streams  mix.  The 
mixing  process  is  assumed  to  occur  with  no  wall  friction. 
As  a  result,  the  total  angular  momentum  of  the  flow  and 
the  total  "stream  thrust"  are  conserved  (Appendix  G).  The 
stream  departing  the  mixing  region  is  assumed  to  have 
uniform  axial  velocity,  w3  ,  and  a  solid  body  rotation 
swirl  velocity  profile. 

f.  The  mixed  incompressible  stream,  "3",  passes 
through  stator-2,  an  ideal  stator  that  removes  the  swirl 
velocity  component.  The  total  pressure  is  assumed  to  be 
constant  along  streamlines  through  this  ideal  stator. 

g.  The  mixed  stream  without  swirl,  stream  "4",  passes 
through  an  ideal  nozzle  to  the  exit,  "e”,  which  exhausts 
to  uniform  static  pressure  of  p#  .  Analysis  is  performed 
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for  both  incompressible  flow  through  the  nozzle  and 
isentropic  flow  of  a  perfect  gas  through  the  nozzle. 


Special  attention  was  used  in  the  selection  of  input 
variables  for  calculation  of  the  performance  of  the  ideal 
incompressible  constant  area  mixer  with  swirl  so  that  the 
same  input  variables  could  be  used  for  calculation  of  the 
performance  of  the  mixer  without  swirl.  The  input 
variables  chosen  are  PTq#P0,  wq  ,  PT2,p2,a,  Rt/rh,  P.  »  Pd /P2 
,  Y  (for  isentropic  nozzle  calculations)  and,  St  ,  the  tip 
swirl  ratio  departing  stator-1  (  St*  v^/wg  ).  This  choice 
of  input  variables  fixes  the  axial  velocity  stream  "2" 
and  the  tip  to  hub  ratio  of  the  inner  annulus,  f,  as 
given  by  the  following  equations: 


w. 


V77T 

*d  '  *2 


f  =  - 


R? 

— + 

p  W 

0  Q  a 

ft  a 

rh2 

p,w. 
r2  2 

rh 

l  + 

pw 

-2-2  a 

P2w2 

+  ^_i«02  ( i+sf ) 

p„  /Pj 


1/2 


1/2 


(2-1) 


(2-2) 


Equations  2-1  and  2-2  show  that  for  increases  in  the 
magnitude  of  the  tip  swirl  ratio  departing  stator-l,  s,  , 
with  all  other  input  variables  constant,  the  value  of 


axial  velocity  of  stream  "2",  w2 ,  increases  and  the  value 
of  the  tip  to  hub  ratio  of  the  inner  annulus,  f, 
increases.  w2  increases  with  increasing  |st|  (  Pd  decreases 
as  St  increases)  in  order  to  satisfy  the  modified  Kutta 
condition  (ratio  of  pd  /P2  )  at  the  end  of  the  splitter 
plate  by  decreasing  P2 .  The  value  of  the  tip  to  hub  ratio 
of  the  inner  annulus,  f,  increases  with  increasing 
magnitude  of  St  to  give  the  ratio  of  outer  annulus  area 
to  inner  annulus  area  required  for  the  densities,  pQ  and 
P2,  the  axial  velocities,  w0  and  w2  ,  and  the  mass  flow 
ratio,  a . 

The  tip  swirl  ratio  of  the  mixed  stream  departing  the 
mixer,  v.  /w, ,  is  fixed  by  the  conservation  of  angular 
momentum.  The  relationship  between  the  two  tip  swirl 
ratios  is: 

v  2  g  f  (  g2  -  1  ) 

-3  *  s.  - - - - - -  (2-3) 

w3  (l+a  )  (i+a^/p2  )  (g2+l)  (f2-l) 

where  "g"  is  the  overall  tip  to  hub  ratio,  Rt/rh 

Equation  2-3  shows  that  vt3/w3  increases  with  increasing  st 
for  constant  value  of  all  other  variables.  Equations  2-1 
and  2-2  show  that  f  increases  with  increasing  St  , 
however,  this  change  is  small,  resulting  in  the  quantity 


(f2-l)/f  increasing  less  than  the  increase  in  st . 


The  axial  force  on  stator-l,  F1  ,  divided  by  the  total 
mass  flow  rate  and  the  axial  velocity  of  the  inner  stream 
(see  Appendix  B)  is 


(m0+m2 )  w0 


,  f  ln(f) 

• 2  _ _ 

(l+  a  )  C f 2— l) 


(2-4) 


A  plot  of  the  variation  of  F1  with  st  is  shown  in  Figure 
2-2  for  typical  input  values.  The  axial  force  on 
stator-2,  F2  ,  divided  by  the  total  mass  flow  rate  and  the 
axial  velocity  of  the  inner  stream  (see  Appendix  B)  is 

(m0+m2)w0  w2  /  pQw0  \  \  g2/  w2  \  g2 
24  (1+^—  a)  l 

\  P2W2  / 

The  thrust  on  stator-2  increases  with  an  increase  in 
magnitude  of  st  ,  however,  the  increase  in  the  thrust  of 
stator-2  is  less  than  the  increase  in  the  magnitude  of 
the  thrust  of  stator-l  and  the  combined  thrust  on  the 
stators,  F1  +  F2 ,  is  negative  and  its  magnitude  increases 
with  increases  in  the  magnitude  of  st  .  Plots  of  the 
variations  of  F2  and  F1  +  F2  with  St  are  shown  in  Figure  2-2 
for  the  same  input  values  as  used  for  the  calculation  of 


The  variations  o£  the  static  pressure  at  stations  "0", 
*1",  "2",  "3"  and  "4"  across  the  annulus  for  values  of  S( 
equal  to  0.48  and  0.68  are  shown  in  Figures  2-3  and  2-5, 
respectively.  Also,  the  corresponding  variations  of  the 
axial  velocity  and  swirl  velocity  at  these  same  stations 
are  shown  in  Figures  2-4  and  2-6,  for  St  equal  to  0.48 
and  0.68,  respectively,  input  variables  utilized  in 
Figures  2-,2  through  2-6  are: 

P,  *  101,000  Pa 
o 

PQ  *P2»  1.23  kg/m 

w  =»  30  m/sec 
o 

a  »  l 

Comparison  of  Figures  2-3  and  2-5  and  Figures  2-4  and  2-6 
show  that  an  increase  in  st  results  in: 

a)  v,  increases  with  a  resultant  decrease  in  p, 

b)  P2  decreases  with  a  resultant  increase  in  w2 

c)  w3  increases,  v3  increases  and  p3  decreases 

d)  w4  increases  and  p4  decreases 


The 

total 

pressure  profiles  for 

station  "4",  PT^  ,  for 

both 

values 

of  St  ,  0.48 

and  0.68, 

are  plotted  in  Figure 

2-7. 

This 

figure  shows 

that  as 

St  increases,  then 

PT/PT  »  1.003 
T0  t2 

Rt/rh  »  2.0 
P2  /Pd  -  1.0 


with  respect  to  the 


decreases  and  the  change  in  PT^  , 
dimensionless  radius, r  /rh ,  increases. 


The  exhaust  velocity,  w# ,  resulting  from  flow  through  an 
incompressible  nozzle  to  uniform  static  pressure,  pa  ,  is 
given  by  the  following  equation: 


and  is  the  incompressible  stream  function 


which  is  related  to  the  dimensionless  radius, 
r/rh ,  at  station  "4"  by  the  equation 
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2w| 


r  a 

r2 

< 

+  — * 

-  r4 

1-h2  — 

1-h4  — 

2w§ 

(2-7) 


The  axial  velocity  at  station  "4",  w4  ,  is  plotted  in 
Figure  2-8  versus  the  dimensionless  radius,  r  /rh ,  for 
the  two  values  of  st  ,  0.48  and  0.68.  The  exhaust 
velocity,  wa  ,  is  also  plotted  in  Figure  2-8  versus  the 
dimensionless  radius  at  station  "4",  r  /rh ,  for  the  two 
values  of  s  and  for  the  two  values  of  P# /PT  ,  0.9  and 
0.99.  Although  w4  is  higher  for  the  higher  value  of  St  ,  w9 
is  reduced  for  increased  st  for  either  value  of  P#  /P^ 


because  the  total  pressure  at  station  ”4",  PT  ,  decreases 


as  s(  increases. 


The  axial  thrust  o£  this  nozzle,  fn ,  is  given  by  the 
following  relationship  for  incompressible  flow: 


w3  -  w3 

•t  *h 


(m0+ft2)w0  3  w0v2(i-h2) 


(2-8) 


where  w#{  is  the  exhaust  velocity  at  the  tip,  <j»  *0 
and  w#  is  the  exhaust  velocity  at  the  hub,  4>»(l-h  )  vt 

h  * 

Figure  2-8  shows  that  as  St  increases,  so  does  the 

difference  between  the  exhaust  velocity  at  the  tip,  w#<  , 

and  the  exhaust  velocity  at  the  hub,  w#  . 

h 

The  net  thrust  of  the  incompressible  swirl  mixer  and 
nozzle,  Fsw  ,  is  the  sum  of  the  thrust  of  stator-l, 
stator-2  and  the  nozzle. 


Fsw-  F,+  FJ+  Fk 


(2-9) 


As  described  above,  the  sum  of  the  thrust  of  both 
stators,  F1  +  F2 ,  decreases  as  st  increases,  as  does  the 
thrust  of  the  nozzle,  FN  .  Thus,  the  net  thrust  of  the 
incompressible  swirl  mixer  and  nozzle,  F$w ,  decreases  as 
the  swirl  is  increased  with  the  other  input  varialbes 


held  constant 


The  exhaust  velocity,  w#  ,  for  isentropic  flow  of  a 
compressible  (perfect)  gas  to  uniform  static  pressure,  ^ 
,  is  given  by  the  following  equation: 


where  $ 


2  »  vt2a 

w3  R? 


and  the  stream  function,  V , 
dimensionless  radius,  r/rh , 
Equation  2-7. 


is  related  to  the 
at  station  “4"  by- 


The  exhaust  velocity  for  compressible  isentropic  flow 
through  a  nozzle  is  greater  than  that  for  incompressible 
isentropic  flow  through  a  nozzle  operating  with  the  same, 
inlet  conditions  and  the  same  exhaust  pressure.  The 
energy  required  for  the  higher  exhaust  velocity  of  the 
compressible  flow  comes  from  the  internal  energy  of  the 
gas.  Comparison  of  Equations  2-6  and  2-10  reveals  the 
different  velocities  calculated  for  the  two  cases.  For 
the  isentropic  flow,  the  quantity  within  the  large 
brackets  is  increased  by  the  quantity  (p#  /P4  )£  1-  (P4 
and  the  quantity  within  the  small  square  root  symbol  is 
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multiplied  by  the  quantity  (Y/(Y-1))#  both  increasing  the 
numerical  value  of  the  exhaust  velocity. 


The  axial  thrust  of  the  nozzle  alone,  with  compressible 
flow,  ,  is  given  by 


fhc  Y-l 

-  3 

(m0+m2)w0  3  Y 


w  v£ 

0  l3 


(1-h2  ) 


(2-11) 


where  w,{  is  the  exhaust  velocity  at  the  tip,  <p  »0 

and  w.  is  the  exhaust  velocity  at  the  hub,  <M(l-h2)v2 
*h  ‘3 


This  nozzle  thrust  increases  slightly  with  increasing  S, 
for  some  values  of  the  input  variables,  decreases  for 
some  others  and  shows  a  maximum  for  others.  The  upstream 
interactions  are  complex,  indeed,  for  variations  in  swirl 
parameter,  so  it  is  not  simple  to  explain  such  results. 
The  results  for  the  overall  performance  follow  intuitive 
reasoning,  however,  as  evidenced  in  Figures  2-10  through 
2-13  and  the  following  discussion. 


The  nozzle  thrust  is  plotted  in  Figure  2-9  versus  st  for 
different  values  of  the  total  pressure  ratio,  PTo/pT2  and 
the  following  values  of  the  other  input  variables: 


19 


101,000  Pa 


P,  /P_  =  0 . 8  a=  1 
To 


P0  =  P2  =  1.23  kg/ra  P2  /Pd  =  1.0  Y=  1.4 


w0  =  30  ra/sec 


Rt/rh  =  5.1 


Although  the  nozzle  thrust  varies  slightly  as  st 
increases,  the  net  thrust  of  the  swirl  mixer  plus  nozzle 
decreases  with  increasing  St  .  This  decrease  in  the  net 
thrust  with  increase  in  St  is  due  to  the  dominance  of 
change  in  the  net  thrust  of  the  two  staters,  stator-1  and 
stator-2,  over  the  slight  variation  in  the  thrust  of  the 
nozzle. 


III.  CONSTANT  AREA  MIXER  WITHOUT  SWIRL 


The  analysis  of  an  ideal  incompressible  constant  area 
mixer  without  swirl  is  contained  in  Appendix  B.  The 
physical  configuration  of  this  mixer  is  the  same  as  shown 
in  Figure  2-1  with  the  two  stators  removed.  Without 
stator-1  and  stator-2,  states  "0"  and  "1"  are  the  same 
and  states  "3"  and  "4"  are  the  same.  The  Kutta  condition 
at  the  end  of  the  splitter  plate  requires  that  the  static 
pressures  of  both  entering  streams,  "0”  and  "2" ,  must  be 
equal.  Using  the  same  input  variables  as  are  used  for  the 
mixer  with  swirl  (  PT^,prt,  wft  ,  PT  ,  p  ,  a  ,  Rt/rh,  P,  and  Y  )  , 
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the  axial  velocity  of  stream  "2"  and  the  tip  to  hub 
radius  ratio  are  given  by  the  following  relationships: 


Wi  ■  JPt>  "  Pt»]+  t  v° 


(2-12) 


f  -  — 


••2  A  (.f 


P  W 

1  +  -2-la 


(2-13) 


The  total  pressure  at  station  "3"  is  given  by  the 


following  relationship: 


t3  to  2 


PoWQ2  (1+a)  (1+a  Pq/Qj) 

p  -  — —  . 


r  p  w 

,  J.  0  0  „ 

1  +  —  a 


F*  p  w_ 

_  o  0 
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,  +P°W°a 
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(PT  -PT  ) 
T0  t2 


(2-14) 


The  thrust  for  incompressible  flow,  F-T  ,  is  given  by 


Fjt  3  ( +  m2 )  w# 


(2-15) 


where 


W  s 


2(PT3-  P.  )/P3 


and  p3  *  PQ(1+ o  )/(l+apo/p2) 


The  thrust  for  compressible  flow,  FSTC  ,  is  given  by 
Equation  2-15,  above,  but  with  the  exhaust  velocity,  w#  , 
given  by  the  following  relationship: 


( 2Y  p_  r 


p_  /p-V/y! 


2-16 


Y-l  P. 


where  P3  is  given  above  and 
1  (1+ot)  (l+apQ/p2  ) 

p  3  p  —  —  P  ■  I  I,-  — » 

3  T3  2  0  0  (l+ap0w0/(p2  w2  )  ) 

The  exhaust  velocity  and  the  thrust  of  the  compressible 
flow  are  higher  than  those  for  incompressible  flow. 
Higher  nozzle  thrusts  are  obtained  with  the  highest  value 
of  the  total  pressure  of  the  mixed  stream,  pt  .  Equations 
2-12,  2-13  and  2-14,  above,  show  that  for  a  fixed  value 
of  p_  ,  w„  ,  p„ ,  p.  and  a ,  PT  increases  with  increasing 

Tq  ’  0  0  2  t3 

values  of  p_  . 

t2 


IV.  COMPARISON  OF  CONSTANT  AREA  MIXERS 


This  section  compares  the  performance  of  an  ideal 
incompressible  constant  area  mixer  with  swirl  to  one 
without  swirl.  The  ratio  of  the  thrust  of  the  mixer  with 
swirl,  Fsw  ,  to  one  without  swirl,  FST  ,  was  calculated  for 
different  values  of  S  with  the  same  values  of  the  input 
variables  (  p  ,  P0 ,  w0  ,  P^  ,  P2  ,  o  ,  Rt  /r„  and  p#  ). 


The  working  equations  listed  in  Section  XI  of  Appendix  B 
were  programmed  on  the  University  of  Washington's  CDC 
computer.  The  results  for  several  examples  are  plotted  in 
Figures  2-10  through  2-13  for  the  following  constant 
values  of  input  variables: 
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Pt0  = 

101/000  Pa 

P./PT  *  0.8 
•  To 

P  =1.23  kg/m 

wo  ■ 

30  m/sec 

Rt/rh  =2.0 

P2/Pd  =1.0 

Figures  2-10  and  2-12  are  for  alpha  ( a  )  equal  to  one  and 

Figures  2-11  and  2-13  are  for  alpha  ( a  )  equal  to  five. 

Figures  2-10  and  2-11  are  for  p  equal  to  1.23  kg/m 

(p0  /p,  =1)  and  Figures  2-12  and  2-13  are  for  p0  equal  to 

0.615  kg/m  (p  /p  =0.5).  Each  figure  contains  a  plot  of 
0  2 

the  thrust  ratio  of  the  swirl  mixer  to  a  straight  mixer 
(Fgyy  FSJ)  versus  the  tip  sw i r  1  ratio  (  St  =  v,/wo  ),  a  plot 
of  the  tip  to  hub  ratio  of  the  inner  channel  (  f  =  rt  /rh) 
versus  S,  and  a  plot  of  the  axial  velocity  at  station 
"2"  (  w3  )  versus  St  .  Each  plot  has  curves  plotted  for 

different  values  of  the  total  pressure  ratio,  p_  /p  . 

'o  T2 

Any  of  these  four  figures  shows  that  an  increase  in  St 

reduces  the  thrust  ratio  and  increases  w2  ,  resulting  in 

an  increase  of  f.  These  figures  also  show  that  an 

increase  in  the  total  pressure  ratio,  P_  /p_  ,  will 

*0  *2 

increase  the  thrust  ratio  and  decrease  w  ,  with  resulting 
decrease  in  f. 

This  comparison  of  the  ideal  ir compressible  constant  area 
mixer  with  swirl  with  the  ideal  incompressible  constant 
area  mixer  without  swirl  shows  that  the  introduction  of 
swirl  in  the  inner  stream  and  its  subsequent  removal 
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reduces  the  available  thrust.  Figures  2-10  through  2-13 

show  that  the  penalty  for  introduction  of  swirl  to  the 

inner  stream  can  be  reduced  by  decreasing  the  total 

pressure  of  stream  "2",  p_  ,  toward  its  minimum  (  w2  =  3  ). 

The  change  in  the  cross-section  of  the  entrance  (for 

mixers  having  the  same  mass  flow  rates)  resulting  from  a 

decrease  in  p  with  all  other  input  variables  fixed  is 
T2 

shown  below.  The  flow  cross-sectional  area  for  the  inner 

stream  is  the  same  in  both  cases,  but  the  area  for  the 

outer  stream  increases  with  decreasing  p  because  w2 

t2 

decreases.  I 


The  hub  radius  of  the  mixer  increases  with  decreasing  E> 

'3 

to  maintain  the  same  overall  hub  to  tip  ratio  and  the 
swirl  velocity  (  vt  )  of  the  mixed  stream  decreases.  The 
swirl  velocity  (  v,  )  of  the  mixed  stream  is  expressed 
analytically  below.  It  is  clear  that  v.  decreases  with 
increasing  hub  radius  because  f  (tip  to  hub  ratio  of  the 
inner  stream)  decreases  as  indicated  in  the  drawings 
above  and  in  Figures  2-10  through  2-13.  The  lower  value 


2k 


of  v^  reduces  the  difference  between  the  tip  and  hub 
total  pressures  at  station  "4".  Further,  as  a  result,  the 
nozzle  thrust  is  penalized  less  because  of  the  reduced 
variation  of  stagnation  enthalpy  across  streamlines. 


vt(  1+a  g2+l 


Comparison  of  Figures  2-10  and  2-11  or  Figures  2-12  and 
2-13  shows  that  an  increase  in  the  mass  ratio ,  a  (»m, /m0 ) , 
reduces  the  loss  in  thrust  of  the  mixer  with  swirl  as 
compared  to  the  mixer  without  swirl.  This  is  as  expected 
because  an  increase  in  the  mass  ratio  reduces  the  portion 
of  the  total  mass  flow  to  which  swirl  is  added. 

Comparison  of  Figures  2-10  and  2-12  or  Figures  2-11  and 
2-13  shows  that  the  density  ratio,  P0/Oj  ,  of  0.5  has  less 
loss  in  thrust  than  does  the  density  ratio  of  1.0.  This 
again  is  expected  because  the  portion  of  the  total  mass 
flow  to  which  swirl  is  added  is  lower  for  the  density 
ratio  of  0.5. 

An  analysis  was  also  performed  for  overall  tip  to  hub 
ratios  other  than  2.0  and  this  analysis  shows  that  an 
increase  in  tip  to  hub  ratio  decreases  the  thrust  the 
thrust  ratio  of  the  swirl  mixer  to  a  straight  mixer.  This 


decrease  in  thrust  ratio  is  mainly  due  to  the  increase  in 
magnitude  of  the  thrust  of  stator-l  with  an  increase  in 
overall  tip  to  hub  ratio. 
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CHAPTER  3 


EXPERIMENTAL  MEASUREMENTS 

I.  INTRODUCTION 

The  mixing  of  two  coaxial  air  streams  (inner  stream  with 
swirl)  was  measured  in  four  set  of  experiments.  These 
sets  of  experiments  used  the  test  apparatus  and 
measurement  systems,  described  in  Chapter  4,  and  the  two 
five-hole  probes,  described  in  Chapter  5.  The  first  three 
sets  of  experiments  were  done  in  test  sections  that  had 
no  outside  wall  and  used  the  initial  measurement  system 
and  the  first  five-hole  probe,  probe  #1.  The  fourth  set 
of  experiments  was  performed  with  an  outside  wall  on  the 
test  section  and  used  the  improved  measurement  system  and 
the  second  five-hole  probe,  probe  #2. 

The  first  three  sets  of  experiments  were  run  to  get  an 
understanding  of  the  mixing  of  the  two  coaxial  air 
streams  when  no  outside  wall  is  present.  The  initial  sets 
of  experiments  also  provided  a  learning  environment  for 
the  final  set  of  experiments,  mixing  of  two  coaxial  air 
streams  in  a  constant  area  mixer. 

The  following  sections  of  this  chapter  describe  the  test 


arrangement  for  each  set  of  experiments  and  the  data 
gathering  procedures  used  in  the  experiments.  Reduction 
of  the  experimental  data  for  these  experiments  is 
described  in  Chapter  6  and  the  results  are  presented  and 
discussed  in  chapter  7. 

II.  FIRST  THREE  SETS  OF  EXPERIMENTS 

The  first  set  of  experiments  was  performed  with  the  test 
section  in  the  configuration  shown  in  Figure  3-1.  In  this 
configuration,  a  single  swirling  stream  of  air  from  a 
circular  annulus  in  a  flat  wall  flows  into  stagnant 
surroundings,  probe  #1  was  used  at  three  axial  locations 
to  make  experimental  measurements  for  one  mass  flow  rate. 
Table  3-1  lists  the  conditions  for  the  three  test  runs  in 
this  set  of  experiments. 

The  configuration  of  the  test  section  for  the  second  set 
of  experiments  is  shown  in  Figure  3-2.  In  this  set  of 
experiments,  two  streams  of  air  (inner  stream  with  swirl) 
from  concentric  annuli  in  a  flat  wall  flow  into  stagnant 
surroundings.  Experimental  measurements  were  made  at 
three  axial  locations  for  a  mass  flow  ratio  (a)  of  one 
using  probe  #1.  Table  3-1  lists  the  conditions  for  the 
three  test  runs  in  this  set  of  experiments. 


The  test  section  for  the  third  set  of  experiments  is 
shown  in  Figure  3-3.  in  this  arrangement,  two  streams  of 
air  (inner  stream  with  swirl)  from  concentric  annuli  in  a 
flat  wall  flow  into  an  open  region.  This  set  of 
experiments  differed  from  the  second  set  of  experiment  by 
the  addition  of  an  extension  of  the  inner  wall  of  the 
inner  annulus  (centerbody)  into  the  open  region,  probe  #l 
was  used  to  make  experimental  measurements  at  five  axial 
locations  for  a  mass  flow  ratio  ( a )  of  one.  The 
conditions  for  the  five  test  runs  in  this  set  of 
experiments  are  listed  in  Table  3-1. 

The  measurement  syst«m  used  for  the  first  three  sets  of 
experiments  is  shown  in  Figure  4-19  and  the  set-up  of  the 
data  gathering  system  is  presented  in  Table  3-2.  A 
Scanivalve  with  a  single  pressure  transducer  and  48  ports 
was  used  in  all  three  sets  of  experiments  to  measure  the 
pressure  at  the  five  ports  of  probe  #1,  the  total 
pressure  of  both  plenums,  and  the  static  pressures  on  the 
outside  wall  of  the  inner  annulus,  in  the  third  set  of 
experiments,  the  static  pressure  along  the  centerbody 
were  also  measured.  The  scanivalve,  with  one  pressure 
transducer,  was  able  to  measure  all  the  static  pressure 
taps  on  the  centerbody  by  changing  the  pressure  taps 


being  measured  from  one  test  run  to  the  next. 

III.  DATA  GATHERING  PROCEDURE  FOR  THE  FIRST  THREE  SETS  OF 
EXPERIMENTS 

The  initial  measurement  system,  shown  in  Figure  4-19,  was 
used  for  the  first  three  sets  of  experiments.  The 
procedure  for  taking  data  was  also  the  same  for  the  first 
three  sets  of  experiments  and  consisted  of  the  following 
steps: 

1.  Required  choked  orifice  plates  were  installed  in 
each  supply  line. 

2.  Supply  voltage  from  each  signal  conditioner  was 
checked. 

3.  Reference  position  of  probe  fl  was  set  and  the 
output  voltage  corresponding  to  r,  z  and  ap  were 
recorded. 

4.  The  scanivalve  was  stepped  from  port  #1  through 
port  #48  while  the  output  of  signal  conditioner  "F"  and 
the  output  of  the  voltage  staircase  from  the  Controller 
were  monitored. 

5.  The  output  of  signal  conditioner  "D"  was  adjusted 
to  zero  for  a  pressure  difference  across  the  differential 
pressure  tranducer  of  zero. 

6.  pressure  regulating  valve  V-l,  shown  in  Figure  4-1, 


was  adjusted  to  give  the  required  supply  line  pressure. 

7.  probe  #l  was  positioned  in  the  radial  and  axial 
location  of  interest. 

8.  probe  #1  was  rotated  to  give  an  output  voltage  of 
about  zero  from  signal  conditioner  "D". 

9.  Time,  supply  line  pressure,  and  tank  pressure  were 
manually  recorded. 

10.  The  Controller  was  given  a  "start"  signal  and  the 
Controller  stepped  the  Scanivalve  through  the  48 
positions.  The  Datalogger  received  a  "scan"  command  from 
the  Controller  about  half  a  second  after  the  Scanivalve 
completed  each  step  and  then  all  eight  channels  were 
scanned,  measured,  recorded  on  magnetic  tape,  and 
displayed  on  the  CRT  terminal  (see  Figure  4-19).  A 
complete  scan  of  all  48  positions  took  about  60  seconds. 

11.  Steps  #7  through  #10,  above,  were  repeated  for  each 
desired  measurement  location  of  probe  #1  for  the  test 
run. 

IV.  FOURTH  SET  OF  EXPERIMENTS 

The  final  and  most  extensive  set  of  experiments  was 
performed  with  the  test  section  in  the  configuration 
shown  in  Figure  3-4.  In  this  configuration,  two  streams 
of  air  (inner  stream  with  swirl)  from  concentric  annuli 


flow  into  a  constant  area  annular  duct,  probe  #2  was  used 
to  make  experimental  measurements  at  seven  axial 
locations  for  mass  flow  ratios  of  0.47,  1.0,  2.13,  and 
3.91  and  at  eight  axial  locations  for  flow  from  only  the 
inner  annulus.  Table  3-3  lists  the  conditions  for  the  36 
test  runs  in  this  set  of  experiments.  The  improved 
measurement  system  used  for  this  set  of  experiments  is 
shown  in  Figure  4-18  and  4-20  and  its  set-up  is  presented 
in  Table  3-4.  All  static  pressures  in  the  test  section 
were  measured  by  the  Scanivalve  with  two  pressure 
transducers.  The  pressure  at  each  of  the  five  ports  of 
probe  #2  was  measured  by  individual  pressure  transducers 
as  described  in  Chapter  4. 

V.  DATA  GATHERING  PROCEDURE  FOR  THE  FOURTH  SET  OF 
EXPERIMENTS 

The  improved  measurement  system,  shown  in  Figures  4-18 
and  4-20,  was  used  for  this  set  of  experiments.  The  data 
taking  procedure  for  this  set  of  experiments  consisted  of 
the  following  steps: 

1.  probe  #2  and  its  section  of  outside  wall  (see 
Figure  4-17)  were  installed  at  the  desired  axial 
location. 

2.  Required  choked  orifice  plates  (see  Table  3-3)  were 


installed  in  each  supply  line. 

3.  Probe  #2  was  moved  to  a  known  radial  position.  The 
output  voltage  from  the  linear  variable  differential 
transformet  (LVDT)  corresponding  to  this  radial  position 
was  recorded. 

4.  The  output  of  signal  conditioners  "j",  "L",  and  "m" 
were  recorded. 

5.  The  length  of  tubing  connecting  Probe  #2  to  its 
block  of  five  pressure  transducers  was  disconnected  so 
that  atmospheric  pressure  was  applied  to  both  ports  of 
each  pressure  transducer. 

6.  The  Scanivalve  was  moved  to  its  "home"  position 
(location  #48) . 

7.  The  Datalogger  was  set  to  continuously  scan  channel 
"01"  and  then  the  pressure  regulating  valve  VI  (see 
Figure  4-18)  was  manually  positioned  to  give  the  desired 
line  pressure  (see  Table  3-3). 

8.  The  Datalogger  was  set  to  scan  channels  "00" 
through  "17"  on  receipt  of  a  remote  "scan"  signal  and 
send  the  measured  values  to  the  Data  Tape  Recorded  and 
CRT  monitor. 

9.  probe  #2  was  positioned  at  the  desired  radial 
location. 

10.  The  "scan"  button  on  the  front  panel  of  the 
Controller  was  pressed  to  initiate  scan  #l. 


11.  The  length  of  tubing  connecting  probe  #2  to  its 


block  of  five  pressure  transduccers  was  reconnected. 

12.  probe  #2  was  rotated  while  monitoring  the  output  of 
signal  conditioner  "F"  (using  a  digital  volt  meter)  until 
the  output  voltage  reached  a  minimum. 

13.  The  "start"  button  on  the  front  panel  of  the 
Controller  was  pressed  to  initiate  the  automatic  stepping 
of  the  Scanivalve  through  its  48  positions.  The 
Datalogger  received  a  "scan”  command  from  the  Controller 
about  a  half  second  after  the  Scanivalve  completed  each 
step  and  all  seventeen  channels  were  scanned,  measured, 
recorded  on  magnetic  tape,  and  displayed  on  the  CRT 
terminal  (see  Figure  4-20) .  A  complete  scan  of  all  48 
Scanivalve  positions  took  about  75  seconds. 

14.  At  the  completion  of  the  automatic  stepping  of  the 
Scanivalve,  the  length  of  tubing  connecting  probe  #2  to 
its  block  of  five  pressure  transducers  was  disconnected 
and  the  "scan"  button  on  the  front  panel  of  the 
Controller  was  pressed  to  initiate  scan  #49. 

15.  Steps  #9  through  #14,  above,  were  repeated  for  each 
radial  location  in  the  test  run. 

16.  Steps  #2  through  #15,  above,  were  repeated  for  each 
test  run  (mass  flow  ratio)  at  that  axial  position. 

17.  Steps  #1  through  #16,  above,  were  repeated  for  each 


axial  location  of  probe  #2  in  this  set  of  experiments. 
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TABLE  3-1 

TEST  RUNS  FOR  FIRST,  SECOND,  AND  THIRD 
SETS  OF  EXPERIMENTS 


Set 

Run 

No. 

l  Orifice  Size* 

Line 

pressure 

(psig) 

Probe 

Axial 

Location 

Name  of  Data  Files 

(inside 

Outside 

Raw 

Fine 

Output 

1/1 

2 

1 

500 

0.28" 

RAW11 

D1A 

PLT1A 

1/2 

2 

1 

500 

1.01" 

RAW  12 

DIB 

PLT1B 

1/3 

2 

1 

500 

2.01" 

RAW13 

D1C 

PLT1C 

2/1 

2 

2 

300 

0.24" 

RAW21 

D2A 

PLT2A 

2/2 

2 

2 

300 

1.00" 

RAW  2  2 

D2B 

PLT2B 

2/3 

2 

2 

300 

4.01" 

RAW23 

D2C 

PLT2C 

3/1 

2 

2 

500 

0.26" 

RAW31 

D3A 

PLT3A 

3/2 

2 

2 

500 

1.00" 

RAW32 

D3B 

PLT3B 

3/3 

2 

2 

500 

4.01" 

RAW33 

D3C 

PLT3C 

3/4 

2 

2 

500 

8.88" 

RAW34 

D3D 

PLT3D 

3/5 

2 

2 

500 

23.99" 

RAW35 

D3E 

PLT3D 

*  Orifice  size 
1 
2 


m  at  600  psia 
0.000  lb/sec 
0.493  lb/sec 


TABLE  3-2 


DATALOGGER/TRANSDUCERS/SIGNAL  CONDITIONERS 


USED  IN  FIRST  THREE  SETS  OF  EXPERIMENTS 


DATALOGGER 

1  SIGNAL  CONDITIONER 

TRANSDUCER 

CHANNEL  SCALE 

VARIABLE 

UNIT 

VOLTAGE 

00 

400  mV 

r 

LVDT 

Conditioner 

- 

LVDT 

01 

40  V 

Scanivalve 
Port  # 

Controller 

Output 

— 

Scanivalve 

Encoder 

02 

40  mV 

z 

A 

10  VDC 

10k  Ohm 
potentiometer 

03 

4  V 

T, 

B 

10  VDC 

Thermistor 

04 

4  V 

T0 

C 

10  VDC 

Thermistor 

05 

400  mV 

P  -P 

1  3 

D 

5  VDC 

Statham  PM- 
283TC+0 .15 

06 

4  V 

aP 

E 

10  VDC 

10k  Ohm 

07 

40  mv 

Px 

F 

5  VDC 

Statham  pm- 
131TC+2 • 5 

Variables  PT  and  PL  were  read  from  gauges  on  the  control 
panel  and  manually  recorded.  The  pressure  drop  across 
each  flow  orifice  plate,  ap,  and  ap0  ,  were  recorded  on 
the  strip  chart  recorder  in  the  control  panel. 


TABLE  3-3 


TEST  RUNS  FOR  THE  FOURTH  SET  OF  EXPERIMENTS 


Run 

No. 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 


Orifice  Size* 


Inside  Outside 


3 

3 

2 

2 

2 

3 

3 

2 

2 

2 

3 

3 

2 

2 

2 

3 

3 

2 

2 

2 

3 

3 

3 

2 

2 

2 

3 

3 

2 

2 

2 

3 

3 

2 

2 

2 


1 

2 

2 

3 

4 
1 
2 
2 

3 

4 
1 
2 
2 

3 

4 
1 
2 
2 

3 

4 
1 
1 
2 
2 

3 

4 
1 
2 
2 

3 

4 
1 
2 
2 

3 

4 


*  Orifice  Size 
1 
2 

3 

4 


Line 


Probe 


Name  of  Data  Files 


mm 

Location 

Raw 

Fine 

Output 

500 

0.25" 

RAW001 

DAT001 

PLT001 

333 

0.25" 

RAW002 

DAT002 

PLT002 

500 

0.25" 

RAW003 

DAT003 

PLT003 

333 

0.25" 

RAW004 

DAT004 

PLT004 

200 

0.25" 

RAW005 

DAT005 

PLT005 

500 

2.25" 

RAW021 

DAT021 

PLT021 

333 

2.25" 

RAW022 

DAT022 

PLT022 

500 

2.25" 

RAW023 

DAT023 

PLT023 

333 

2.25" 

RAW024 

DAT024 

PLT024 

200 

2.25" 

RAW025 

DAT025 

PLT025 

500 

4.25" 

RAW041 

DAT041 

PLT041 

333 

4.25" 

RAW042 

DAT042 

PLT042 

500 

4.25" 

RAW043 

DAT043 

PLT043 

333 

4.25" 

RAW044 

DAT044 

PLT044 

200 

4.25" 

RAW045 

DAT045 

PLT045 

500 

6.25" 

RAW061 

DAT061 

PLT061 

333 

6.25" 

RAW062 

DAT062 

PLT062 

500 

6.25" 

RAW063 

DAT063 

PLT063 

333 

6.25" 

RAW064 

DAT064 

PLT064 

200 

6.25" 

RAW065 

DAT065 

PLT065 

500 

8.25" 

RAW081 

DAT081 

PLT081 

500 

12.25" 

RAW121 

DAT121 

PLT121 

333 

12.25" 

RAW122 

DAT122 

PLT122 

500 

12.25" 

RAW123 

DAT123 

PLT123 

333 

12.25" 

RAW124 

DAT124 

PLT124 

200 

12.25" 

RAW125 

DAT125 

PLT125 

500 

18.25" 

RAW181 

DAT181 

PLT181 

333 

18.25" 

RAW182 

DAT182 

PLT182 

500 

18.25" 

RAW183 

DAT183 

PLT183 

333 

18.25" 

RAW184 

DAT184 

PLT184 

200 

18.25" 

RAW185 

DAT185 

PLT185 

500 

24.25" 

RAW241 

DAT241 

PLT241 

333 

24.25" 

RAW242 

DAT242 

PLT242 

500 

24.25" 

RAW243 

DAT243 

PLT243 

333 

24.25" 

RAW244 

DAT244 

PLT244 

200 

24.25" 

RAW245 

DAT245 

PLT245 

at  600  psia 

0.000  lb/sec 
0.493  lb/sec 
1.052  lb/sec 
1.928  lb/sec 


TABLE  3-4 


DATALOGGER/TRANSDUCERS/SIGNAL  CONDITIONERS 
USED  IN  THE  FOURTH  SET  OF  EXPERIMENTS 


_ DATALOGGER 

CHANNEL! 


SIGNAL  CONDITIONER 


SCALE 

VARIABLE 

UNIT 

VOLT AG 

4  V 

ap 

I 

10  VDC 

40  mV 

PL 

J 

5  VDC 

40  mv 

PT 

K 

5  VDC 

4  V 

r 

LVDT 

Conditioner 

- 

40  V 

Scanivalve 
Port  # 

Controller 

Output 

— 

J 

T, 

- 

- 

J 

T0 

- 

- 

Skip 

- 

- 

- 

40  mv 

P1 

G 

12  VDC 

40  mv 

P2 

H 

12  VDC 

40  mv 

P3 

A 

12  VDC 

40  mv 

P4 

B 

12  VDC 

40  mv 

P5 

C 

12  VDC 

4  0  mv 

Px 

D 

12  VDC 

40  mv 

Px 

E 

12  VDC 

^0  mv 

P,  -P3 

F 

5  VDC 

40  mV 

AE> 

L 

5  VDC 

40  mv 

*Po 

M 

5  VDC 

TRANSDUCER 


00 

01 

02 

03 

04 

05 

06 

07 

08 

09 

10 

11 

12 

13 

14 

15 

16 
17 


10k  ohm 
potentiometer 
Statham  pg- 
146TC-800 
Statham  pg- 
146TC-5M 
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CHAPTER  4 


TEST  APPARATUS 


I.  INTRODUCTION 

The  test  apparatus  consists  of  five  major  portions: 

1.  High  pressure  dual  air  supply  system 

2.  Concentric  dual  plenum  section 

3.  Swirl  generator  and  nozzles  section 

4.  Test  section 

5.  Measurement  equipment 

A  diagram  of  the  high  pressure  air  supply  system  is  shown 
in  Figure  4-1.  Figure  4-2  shows  the  next  three  major 
portions  of  the  test  apparatus  as  listed  above.  The 
following  sections  of  this  chapter  will  discuss  the  five 
major  portions  of  the  test  apparatus  in  the  order  listed. 

II.  HIGH  PRESSURE  DUAL  AIR  SUPPLY  SYSTEM 

The  high  pressure  air  supply  system  consisting  of  an  air 
compressor,  air  receiver  and  single  pressure  supply 
piping  was  modified  to  supply  two  separate  air  streams  at 
different  stagnation  pressures.  This  modification  was 
accomplished  by  connecting  a  new  supply  line  to  the 
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existing  air  receiver  and  installing  the  piping  and 
valving  shown  in  Figure  4-1.  Control  valve  VI  can  be  set 
to  maintain  a  downstream  pressure  of  up  to  600  psig.  The 
medium  pressure  air,  downstream  of  control  valve  VI,  is 
then  piped  into  two  separate  and  controllable  streams. 
Initially,  the  stagnation  pressure  of  supply  air  stream 
#1  was  controlled  by  control  valve  V2  to  maintain  a 
supply  pressure  of  up  to  39  psig  and  its  mass  flow  rate 
measured  by  orifice  plate  01.  Also,  the  stagnation 
pressure  of  supply  air  stream  #2  was  controlled  by 
control  valve  v3  to  maintain  a  supply  pressure  of  up  to 
30  psig  and  its  mass  flow  rate  measured  by  orifice  plate 
02.  This  supply  system  had  severe  stability  problems  at 
the  low  mass  flow  rates  required  for  these  experiments 
and  modification  was  required.  The  two  control  valves,  V2 
and  V3 ,  were  replaced  by  orifice  plates  sized  for 
specific  mass  flow  rates  and  operating  under  choked 
conditions.  Further  discussion  of  the  choked  orifice 
plates  is  included  in  the  "Measurement  Equipment"  section 
of  this  chapter. 

III.  CONCENTRIC  DUAL  PLENUM  SECTION 

The  concentric  dual  plenum  section,  shown  in  Figure  4-3, 
consists  of  a  30-inch  diameter  cylinder  supported  inside 
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a  40-inch  diameter  cylinder  and  this  section  is  90  inches 
long.  The  air  enters  the  inner  plenum  through  a  two-inch 
pipe  that  is  aligned  with  the  centerline  o£  the  plenum 
and  this  air  stream  is  dispersed  by  a  "core  buster".  The 
air  enters  the  outside  plenum  through  four  pipes  of 
two-inch  diameter  spaced  evenly  around  the  annulus  and 
each  entering  air  stream  is  dispersed  by  a  "core  buster". 
Both  the  air  stream  in  the  inner  circular  plenum  and  the 
air  stream  in  the  outer  annular  plenum  then  pass  through 
two  sets  of  screens,  three  sets  of  honeycomb  and  four 
sets  of  screens.  The  two  concentric  air  streams  leaving 
this  dual  plenum  section  are  essentially  uniform  and  have 
a  low  turbulence  level; 

The  concentric  dual  plenum  section  is  connected  to  the 
high  pressure  dual  air  supply  system  by  two  sections  of 
two-inch  diameter  schedule  40  pipe.  These  sections  of 
pipe  are  connected  to  the  dual  plenum  section  by  rubber 
hose  and  connected  to  the  high  pressure  dual  air  supply 
system  by  bolted  flanges.  Each  section  of  pipe  contains  a 
rupture  disc  to  prevent  overpressurization  of  the  dual 
plenum  section. 

The  dual  plenum  section  is  designed  to  handle  stagnation 
pressures  of  up  to  10  psig,  however,  in  these 
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experiments,  the  stagnation  pressures  were  below  1  psig. 
The  design  and  construction  of  the  dual  plenum  section 
was  influenced  by  the  space  limitations  for  the  test 
apparatus  and  the  size  of  the  doorway  leading  into  the 
space  that  houses  the  test  apparatus. 

The  "core  busters"  consist  of  1/2-inch  thick  aluminum 
plates  mounted  to  the  common  plenum  end  cap  by  threaded 
rods  (for  adjustment) .  The  "core  buster"  for  the  inner 
plenum  measures  7"xil"  with  75  holes  of  one-half  inch 
diameter  distributed  evenly  over  its  surface  and  was 
mounted  about  six  inches  away  from  the  plenum  end  cap. 
The  four  "core  busters"  for  the  outer  plenum  each  measure 
4.5"xll"  with  75  holes  of  one-half  inch  diameter 
distributed  evenly  over  each  surface  and  each  was  mounted 
about  four  inches  away  from  the  plenum  end  cap. 

A  set  of  screens  was  placed  upstream  of  the  honeycomb 
sections  to  break  up  the  large  eddies  and  provide  a  more 
even  flow  upstream  of  the  honeycomb.  Commercial  18x14 
mesh  brass  screen  fabric  was  used  for  the  screens  which 
gave  a  mesh  size  of  0.0625  inches.  Using  Figure  1.18  on 
Page  19  of  Reference  1,  this  mesh  size  and  a  mean 
velocity  of  5  ft/sec  give  a  Reynolds  Number  of  163  and 
the  value  of  u  /  u':  reaches  a  value  of  about  5,000  within 
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four  inches  downstream  of  the  screen.  This  dampening  of 
the  turbulence  intensity  is  acceptable  as  a  minimum  and 
the  minimum  distance  between  screens  was  selected  as  four 
inches.  This  agrees  with  the  suggestions  of  Schubauer,  et 
al ,  (Reference  2)  . 

Screens  are  capable  of  reducing  the  magnitude  of  the 
turbulence  intensity  in  the  direction  of  the  mean  flow. 
Honeycomb  is  capable  of  aligning  the  mean  flow  and 
reducing  the  lateral  components  of  turbulence;  however, 
the  longitudinal  component  of  turbulence  is  increased 
(Reference  3).  An  essentially  uniform  and  low  turbulence 
stream  is  obtained  by  having  the  air  stream  first  pass 
through  three  sets  of  honeycomb  to  align  the  flow,  reduce 
lateral  components  of  turbulence  and  reduce  longitudinal 
non-uniformity  in  the  mean  velocity.  The  air  stream  then 
passes  through  two  sets  of  screens  to  reduce  the 
longitudinal  components  of  turbulence.  The  resulting  air 
streams  can  then  be  accelerated  through  nozzles  to  the 
test  section  velocity  and  for  a  contraction  ratio  of  25 
and  45,  the  ratio  of  the  turbulence  intensity  to  the  mean 
velocity  is  reduced  to  0.16  and  0.12,  respectively 
(Reference  4) . 

The  effect  of  multiple  sections  of  honeycomb  on  reducing 


60 


longitudinal  non-uniforrai ty  in  the  mean  velocity  was 
analyzed  using  the  model  shown  in  Figure  4-4.  The  model 
is  composed  of  a  frictionless  duct  (section  1-2)  which 
allows  both  the  high  velocity  and  low  velocity  streams  to 
adjust  before  entering  the  honeycomb  (section  2-3) .  The 
honeycomb  is  assumed  to  be  composed  of  small  circular 
pipes.  The  two  air  streams  are  assumed  to  have  equal 
static  pressures  at  points  1  and  3.  The  sum  of  the  area 
of  the  high  velocity  flow  plus  the  area  of  the  low 
velocity  flow  is  constant  and  the  area  of  the  honeycomb 
is  constant.  For  1/8-inch  diameter  honeycomb,  a  Reynolds 
Number  of  195  corresponds  to  a  mean  velocity  of  3  ft/sec 
and  a  Reynolds  Number  of  1950  corresponds  to  a  mean 
velocity  of  30  ft/sec.  Since  the  design  overall  mean 
velocity  is  5  ft/sec,  the  flow  through  the  honeycomb  is 
assumed  to  be  developing  laminar  flow.  Application  of  the 
continuity  and  momentum  equations  to  the  model  give  the 
following  relationships  for  incompressible  flow: 


(4-1) 


(4-2) 


P,  -P3  =  1/2  P  V* 


64  LP 


d2  V. 


+  2.41 


-  1/2  pVH2  (4-3) 


[64L  V 

P,  -P  =  1/2  p  V,3  - -  +  2.41 

1  3  3  d2  V. 

L3 

The  frictional  loss  term  contained  within  the  brackets  in 
Equations  4-3  and  4-4  was  obtained  from  Equation  44  on 
page  308  of  Reference  5.  Since  the  pressure  difference 


-  1/2  p  (4-4) 


P.  -  P„  of  the  two  streams  are  equal,  Equations  4-2  and 

I  3 

4-3  can  be  equated  to  give  the  following  relationship: 
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Thus  Equations  4-1,  4-2  and  4-5  give  a  system  of 
equations  that  can  be  solved  to  determine  the  downstream 
velocity  discontinuity  and  distribution  given  the 
upstream  velocity  discontinuity  and  distribution. 


A  comparison  between  the  performance  of  three  two-inch 
thick  sections  of  honeycomb  and  one  six-inch  thick 
section  of  honeycomb  was  made  using  the  following  input 
data: 


v* 


Hi- 


18  ft/sec 
3  ft/sec 
9 


d  =  1/8  inch 
p  -  0.075  lb/ft 
V  =  0 .00016  ft/sec 


The  model  gives  the  following  results  for  the  two 
different  configurations: 

Three  -  2"  Sections  one  -  6"  Section 


^  =  4.544  ft/sec 
VL3  =  4  .435  ft/sec 
Al/A^5*  0.68311 

P  =  0.0397  in  H20 


6.o64  ft/sec 
3.840  ft/sec 
2.3687 

0.0129  in  H20 


From  the  above  example,  the  superiority  of  multiple 
sections  of  honeycomb  versus  one  section  of  the  same 
total  thickness  is  evident  in  its  ability  to  even  out 
mean  velocity  variations.  In  both  cases,  the  mean 
velocity  has  an  average  value  of  4.5  ft/sec  with  an 
incoming  standard  deviation  of  4.5  ft/sec.  The  leaving 
standard  deviation  for  the  three  two-inch  sections  of 
honeycomb  is  0.0535  ft/sec  whereas  the  standard  deviation 
for  the  air  stream  leaving  the  one  six-inch  section  of 
honeycomb  is  1.016  ft/sec.  The  multiple  sections  of 
honeycomb  produce  a  larger  pressure  loss  than  a  single 
section  of  honeycomb  with  the  same  total  thickness  of  six 
inches  due  to  the  multiple  entrance  losses. 


Because  only  a  limited  amount  of  two-inch  thick  honeycomb 
with  1/8-inch  cells  was  available,  the  honeycomb  was  used 
to  align  the  flow  and  reduce  the  mean  velocity  variations 


by  using  three  sections  of  honeycomb  in  each  plenum.  The 
position  of  the  three  honeycomb  sections  is  shown  in 
Figure  4-3. 

Downstream  of  the  honeycomb,  two  sets  of  screens  were 
placed  to  reduce  the  longitudinal  component  of  turbulence 
as  recommended  by  Scheiman  and  Brooks  in  Reference  3.  The 
position  of  these  sets  of  screens  in  the  dual  plenum  is 
shown  in  Figure  4-3. 

IV.  SWIRL  GENERATOR  AND  NOZZLES  SECTION 

The  swirl  generator  and  nozzles  section  is  shown  in 
Figure  4-5.  This  section  is  composed  of  the  swirl 
generator,  swirl  generator  housing,  inner  annulus  between 
the  four-inch  OD  center  body  and  the  six-inch  ID 
cantilevered  tube  and  outer  annulus  between  the  outside 
nozzle  and  inner  nozzle. 

The  swirl  generator  was  designed  to  give  a  swirling  air 
stream  with  uniform  axial  velocity  and  free  vortex 
angular  velocity.  This  type  of  velocity  profile  was 
obtained  by  bringing  the  air  stream  radially  into  a 
one-inch  slot,  passing  the  air  stream  over  adjustable 
stators  which  impart  an  angular  velocity  and  then  turn 


the  air  stream  into  an  annular  passage  which  converts  the 
radial  velocity  into  axial  velocity.  An  isometric 
cross-section  of  the  swirl  generator  is  shown  in  Figure 
4-6.  The  outside  diameter  of  the  swirl  generator  is  18 
inches  and  the  air  enters  the  radial  slot  around  the 
perimeter  through  a  bell  mouth-type  entrance  whose  length 
is  1/2  inch.  The  stator  blades  are  positioned  so  that  the 
leading  edge  of  each  blade  falls  on  a  circle  whose 
diameter  is  16  inches;  the  pivot  shaft  of  each  blade 
falls  on  a  circle  whose  diameter  is  14  inches  and  the 
trailing  edge  of  each  blade  falls  on  a  circle  whose 
diameter  is  12  inches.  The  stator  blades  were  made  using 
a  symmetric  NACA  0024  airfoil  shape  whose  centerline  was 
curved  to  give  the  desired  exit  angle.  The  geometry  of 
the  stator  cascade  is  shown  in  Figure  4-7.  Constant's 
rule,  which  is  given  below  in  empirical  form,  was  used  to 
determine  the  desired  departing  stagger  angle  for  a 
turning  angle  of  45.  Constant's  rule,  as  given  on  pages 
179  and  180  of  Reference  6,  is: 

6/0*=  0.25 IJa  (4-6) 

where  (refer  to  Figure  4-8) 

<p  =  8i  -  B2  =  turning  angle  Y1  =  stager  angle 

O  =  c /s  =  solidarity  Q‘=  Y1  -  Y2 

a  =  8,  -Y,  =  angle  of  attack  6=  82-Y2 
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The  following  values  apply  for  the  design  shown  in  Figure 
4-7: 

81  »0,Y1»0,  Y2»-60°,  0=2.32",  s-2.20" 


Thus  O  *1.055  and,  from  Equation  6,  B2=-45.4°.  This  value 
of  departure  angle  is  very  close  to  the  desired  value  of 
-45°  and  was  used  for  the  swirl  generator  design.  The 
stator  cascade  consists  of  twenty  stator  blades  which  can 
be  adjusted  in  unison  to  give  a  value  of  61  ,  from  -30  to 
30.  Each  blade  has  a  symmetric  NACA  0024  shape  about  its 
centerline  whose  radius  is  2.75  inches  (see  Figure  4-7). 
For  the  geometry  shown  in  Figure  4-7  and  an  entering 
radial  velocity  of  35  ft/sec  at  the  blade  leading  edge, 
the  radial  and  angular  velocities  at  the  blade  trailing 
edge  are  46.7  ft/sec  and  47.3  ft/sec,  respectively,  forB2 
=*-45.4°.  The  air  stream  leaving  the  stator  blades  then 
enters  a  channel  as  shown  in  Figure  4-9. 


Application  of  conservation  of  mass  and  angluar  momentum 
to  the  flow  in  the  channel  shown  in  Figure  4-9  gives  the 
following  relationships  for  invisid  incompressible 
axisymmetric  flow  in  radial  equilibrium: 

v3  *  v2  r2/f  ri  *  r  S  r0 
w3  -  2  u2r2/(r02  -r,2) 


(4-7) 

(4-8) 


66 


•  -  .•  •  -v- 


The  values  of  r2  ,  r0  ,  and  r,  for  this  design  are  6 
inches,  3  inches  and  2  inches,  respectively.  Thus  for  u2 
■  46.7  ft/sec  and  v2  *  47.3  ft/sec.  Equations  4-7  and  4-8 

give 

« 

w3  *  112  ft/sec  and 


v3  *  94.6  x  r0/r  ft/sec  2^r<3 


(4-9) 


The  angular  velocity  will  vary  in  accordance  with 
Equation  4-9  from  94.6  ft/sec  at  r  =*  3  inches  to  141.9 
ft/sec  at  r  *  2  inches.  The  angle  that  the  velocity 
vector  makes  with  the  centerline  will  vary  from  40.2° at  r 
=  3  inches  to  51. 7° at  r  ■  2  inches  for  the  above  design 
data. 


The  shape  of  the  channel  which  transforms  the  radial  slot 
into  an  annulus  was  chosen  to  be  circular  in 
cross-section  as  shown  in  Figure  4-9.  Although  other 
transitions  could  have  been  used,  the  circular 
cross-section  was  easiest  to  manufacture  and  gave  the 
required  gradual  transition. 


For  a  design  mass  flow  of  0.916  lb/sec,  the  velocities  at 
each  section  of  the  swirl  generator  are  as  follows: 

Entering  the  swirl  generator  housing  2.49  ft/sec 

Swirl  generator: 

at  entrance  to  bell  mouth  15.56  ft/sec 

at  leading  edge  of  stator  blade  35.00  ft/sec 
at  trailing  edge  of  "  "  u  =  46.7  ft/sec 

v  *  47.3  ft/sec 

at  exit  for  r*3  inches  w  *112.0  ft/sec 

v  *  94.6  ft/sec 

at  exit  for  r=2  inches  w  *112.0  ft/sec 

v  *141.9  ft/sec 

The  annular  channel  has  a  mean  length  of  32.5  inches. 
Because  of  this  length,  an  analysis  of  the  development  of 
the  boundary  layer  was  performed  using  the  waltz  method 
for  turbulent  boundary  layer  for  axisymmetric  flow  as 
described  by  Waltz  in  Reference  7.  The  mean  velocity 
along  a  line  half  way  between  the  channel  walls, 
corrected  for  the  displacement  thickness,  was  used  to 
determine  the  imposed  pressure  gradient  and  the 
properties  of  the  turbulent  boundary  layer.  The  variation 
of  this  mean  velocity  with  the  axial  distance  for  this 
annulus  is  plotted  in  Figure  4-10.  The  boundary  layer 
properties  at  the  exit  from  this  channel  for  the  above 
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design  values  are  as  follows  for  both  walls  of  the 
annular  channel: 


^3  0.098  inches  H 


12 


d2“  0.072  inches 


H 


32 


1.350 

1.763 


The  outer  portion  o£  the  swirl  generator  housing,  shown 
in  Figure  4-5,  is  an  extension  of  the  outer  plenum.  This 
extension  of  the  outer  plenum  is  connected  to  the  section 
containing  the  inside  nozzle  and  outside  nozzle.  The  two 
concentric  nozzles  were  designed  using  the  following 


empirical  relationship  due  to  vitoshinski  in  Reference  8. 

1/2 


.  L .  i\ 

\  rf)  (l  +  *a/LJ) 


(4-10) 


where  r0  ■  outlet  radius  (z  •  L) 
r,  >  inlet  radius  (z  >  0) 
L  *  nozzle  length 


For  design,  a  nozzle  length  of  36  inches  was  used  in 
Equation  4-10  and  the  actual  nozzle  length  was  limited  to 
the  first  29.75  inches.  The  following  values  for  r  and 
r.  were  used  for  the  nozzles: 

Inside  nozzle:  r(  -  15  inches  &  re  ■  3.003  inches 
Outside  nozzle:  r{  ■  20  inches  6  r0  ■  3.95  inches 


The  nozzle  profiles  generated  by  Equation  4-10  agree  very 


.vly. 


well,  over  the  majority  of  the  nozzle  length,  with  the 
profiles  generated  using  the  method  of  Tsien  contained  in 
Reference  9  and  the  profiles  generated  using  the  method 
of  Szczeniowski  contained  in  Reference  10.  as  shown  in 
Figure  4-5,  the  inner  nozzle  wall  and  the  outer  wall  of 
the  inner  annulus  are  machined  from  the  same  piece  of 
material  for  over  half  the  nozzle's  length.  This  common 
wall  ends  at  the  entrance  to  the  test  section  and  its 
thickness  at  this*  point  is  0.033  inches. 

Static  pressure  taps  were  placed  on  the  inside  of  the 
outer  wall  in  the  inner  annulus  to  get  a  measure  of  the 
adjustment  of  the  swirling  flow  as  it  approaches  the 
exit.  The  placement  of  these  static  pressure  taps  and  a 
detail  of  a  typical  pressure  tap  are  shown  in  Figure 
4-11.  stainless  steel  tubing  with  o. 021-inch  id  and 
0.010-inch  wall  thickness  was  used  and  placed  in  grooves 
machined  in  the  outer  wall  of  the  six-inch  id  cantilever 
tube.  The  grooves  in  the  outside  wall  were  filled  before 
final  machining  of  the  inner  nozzle  surface.  Holes  of 
0.021-inch  diameter  were  drilled  through  the  inside  wall 
of  the  six-inch  ID  cantilever  tube  at  all  positions 
listed  below,  except  the  end  positions: 


pressure 

Tap  Number 

Z 

Angular  position 

29  - 

35 

-  41 

end 

90, 

-30, 

-150 

30  - 

36 

-  42 

-0.10" 

70, 

-50, 

-170 

31  - 

37 

-  43 

-0.25" 

50, 

-70, 

170 

32  - 

38 

-  44 

-0.50" 

30, 

-90, 

150 

33  - 

39 

-  45 

-1.00" 

10, 

-110, 

110 

34  - 

40 

-  46 

-2.00" 

-10, 

-130, 

110 

An  analysis  of  the  development  of  the  turbulent  boundary 
layers  on  each  nozzle  was  performed  using  the  Waltz 
method  for  axisymmetric  flow  as  described  by  Waltz  in 
Reference  7.  For  an  entering  velocity  of  4.48  ft/sec,  the 
calculated  boundary  layer  properties  on  both  nozzle  walls 
at  exit  are  as  follows:  ft*  0.045  inches  h  »  1.39 

ft^»  0.032  inches  H32  *  1.74 

The  variation  of  the  mean  velocity  with  the  axial 
distance  is  plotted  in  Figure  4-12. 

V.  TEST  SECTION 


The  test  section  was  designed  so  that  it  could  be  set  up 
in  one  of  three  different  test  configurations.  The  three 
test  configurations  are  shown  in  Figure  4-13  and  listed 
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Configuration 


Three 


Description 


Flow  with  swirl  from  inner  annulus 

with  outside  flat  wall  and  with  or 

without  a  centerbody 

Flow  from  two  annuli  (axial  flow 

from  outer  and  swirl  flow  from  inner) 

with  outside  flat  wall  and  with  or 

without  a  centerbody 

Flow  from  two  annuli  (axial  flow 

from  outer  and  swirl  flow  from  inner) 

with  constant  radius  outer  wall  and 

centerbody 


The  test  section  has  a  length  of  48  inches.  Measurements 
of  the  total  and  static  pressures  and  flow  direction  for 
test  configurations  one  and  two  (without  outside  wall) 
are  possible  using  probe  fl  and  its  traversing  mechanism, 
as  shown  in  Figure  4-14,  which  is  mounted  on  rails 
parallel  to  the  centerline  of  the  test  section.  The 
probe's  traversing  mechanism  allows  for  rotating  the 
probe  face  to  align  the  probe  with  the  flow  direction  and 
radial  positioning  of  the  probe.  The  probe's  radial 
position  is  measured  by  the  linear  variable  differential 
transformer  (LVDT) .  The  probe's  angle  is  measured  by  a 


potentiometer  on  the  traversing  mechanism  and  the  axial 
position  is  measured  by  a  potentiometer  on  the  axial 
traversing  bed. 

An  outer  flat  wall  is  provided  by  a  36  inch  outside 
diameter  piece  of  pressboard  that  mounts  to  the  end  of 
the  outer  nozzle  for  test  configurations  one  and  two.  The 
outer  flat  wall  is  extended  inward  to  the  6-  inch  OD 
annulus  for  test  configuration  one  by  insertation  of  a 
large  disc  into  the  annulus  between  the  nozzle  walls. 

The  4-inch  OD  centerbody  of  the  swirl  generator  can  be 
ended  at  the  beginning  of  the  test  section  by  insertation 
of  an  end  cap  or  extended  through  the  test  section  by 
attaching  the  4-inch  centerbody  sections  shown  in  Figure 
4-15.  The  centerbody,  for  use  in  the  test  section,  was 
made  in  three  pieces:  one  section  of  two-foot  length  with 
twenty  wall  pressure  taps;  one  section  of  two-foot  length 
with  eight  wall  pressure  taps,  and  an  end  support.  All 
three  pieces  of  this  centerbody  are  interchangeable.  As 
shown  in  Figure  4-15,  the  wall  static  pressure  taps  are 
located  at  the  following  axial  positions: 
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V. 


Tap  # 

1 

2 

3 

4 

5 

6 

7 

8 

z  (in) 

0.125 

0.250 

0.375 

0.50 

0.75 

1.00 

1.25 

2.25 

Tap  f 

9 

10 

11 

12 

13 

14 

15 

16 

z  (in) 

3.25 

4.25 

5.25 

6.25 

8.25 

10.25 

12.25 

14.25 

Tap  #  1 

17 

18 

19 

20 

21 

22  23 

z  ( in) | 

16.25 

18.25 

20.25 

22.25 

24.25 

26.25  28.25 

Tap  # 

24 

25 

26 

27 

28 

z  (in) 

30.25 

34.25 

38.25 

42.25 

46.25 

The  static  pressure  tap  holes  in  the  wall  of  the  4-inch 
00  centerbody  have  a  diameter  of  0.021  inches  and  are 
connected  to  stainless  tubing  of  the  same  inside 
diameter. 


The  eight-inch  id  outer  wall  for  test  configuration  three 
is  made  from  plexiglass  tubing  of  l/2-inch  wall  thickness 
that  was  cut  into  five  interchangeable  sections.  Each 
section  has  wall  static  pressure  taps  of  0.021-inch 
diameter  which  are  located  every  two  inches  along  the 
length.  The  length  of  each  section  of  outer  wall  and  the 
pressure  tap  locations  are  v  shown  in  Figure  4-16  and 
listed  below. 


Length 


Tap  # 


47  1 48  1 49  50  51  52  53  54  55  56  57  58 


Locationl  l 


Section 


Length 


Tap  #  U59  60  61  62  1 63  64  65  1 66  1 67  1 68  69  70 


Locationl  1 


Measurement  of  the  total  and  static  pressures  and  flow 
direction  for  test  configuration  three  (with  outside 
wall)  were  performed  using  probe  #2  and  its  traversing 
mechanism  with  section  of  outer  wall,  as  shown  in  Figure 
4-17.  The  section  that  contains  probe  #2  and  its 
traversing  mechanism  has  a  length  of  2.375  inches  and  the 
head  of  the  probe  is  located  0.250  inches  from  the 
entrance  of  the  section.  The  probe  head  can  be  traversed 
radially  between  outside  and  inside  walls  on  the  side  of 
the  centerbody  opposite  the  traversing  mechanism. 
Measurement  at  different  axial  locations  is  accomplished 
by  interchanging  the  sections  of  the  outer  wall  and  the 
probe  section  of  the  outer  wall.  The  progressive  section 
le  -glits  9 1  low  location  of  the  probe  section  at  1/4  inch 
test  section  and  in  increments  of  two  inches 


through  the  test  section.  The  probe  section  of  the  outer 
wall  has  two  static  pressure  taps  1/4”  from  the  leading 
edge  as  shown  in  Figure  4-17. 

The  outside  wall  of  the  test  section  is  supported  at  its 
entrance  by  connection  with  the  end  of  the  outside 
nozzle,  along  its  length  by  supporting  rails  and  at  its 
exit  by  connection  to  the  outlet  section  of  the  outside 
wall  that  is  itself  supported  by  a  spider  connected  to 
the  centerbody  end  support.  The  outlet  section  of  the 
outside  wall  has  a  slight  nozzle  at  its  exit  with  an  exit 
to  entrance  area  ratio  of  0.84  to  provide  slight 
pressuration  of  the  test  section  above  the  pressure  in 
the  surrounding  room. 

VI.  MEASUREMENT  EQUIPMENT 

The  measurement  equipment  is  divided  into  two  systems. 
One  system  consists  of  the  measurement  equipment 
associated  with  the  high  pressure  dual  air  supply  system. 
The  other  system  consists  of  the  measurement  equipment 
associated  with  the  concentric  dual  plenum  section,  the 
swirl  generator  and  nozzles  section  and  the  test  section. 


A.  High  pressure  Dual  Air  supply  system 


The  measurement  equipment  associated  with  the  high 
pressure  dual  air  supply  system  is  shown  in  Figure  4-18. 
Control  valves  V2  and  V3 /  shown  in  Figure  4:1,  have  been 
replaced  by  orifice  plates  sized  to  supply  specific  mass 
flow  rates  for  choked  flow  conditions.  The  empirical 
relationship  for  choked  flow  through  an  orifice  plate  is 
given  on  page  21  of  Reference  6.  This  relationship  solved 
for  the  orifice  area  (A*)  in  terms  of  the  mass  flow  rate 
(m)  and  upstream  gas  properties  is  given  below  for 
isentropic  flow. 


A*  =  m  Jr  tt  /  (r  PT) 


(4-11) 


where  PT  and  Tt  are  the  upstream  stagnation  pressure  and 
temperature,  respectively,  R  is  the  gas  constant  and  T  is 
a  symbol  for  the  equation  written  below  in  terms  of  Y  , 
the  ratio  of  specific  heats. 
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Choked  flow  orifice  plates  were  sized  using  Equation  4-11 
for  mass  flow  rates  of  0.493  lb/sec,  1.052  lb/sec,  1.928 


lb/sec  and  3.15  lb/sec  with  an  upstream  stagnation 
pressure  of  600  psia  and  stagnation  temperature  of  500  R. 
The  actual  mass  flow  rate  was  calculated  using  the 
relationship: 

.  .  1/2 

m  =  m0(PT/PTo  )  (Tt<j/Tt)  (4-12) 

where  the  zero  subscripted  quantities  are  design  values, 
listed  above,  and  the  non-subscripted  quantities  are  the 
actual  values  measured  during  a  test. 

Flow  meter  orifice  plates  were  sized  using  the  procedures 
listed  in  Reference  11  for  a  pressure  drop  of  100  inches 
of  water  and  upstream  stagnation  conditions  of  600  psia 
and  500  R.  Identical  orifice  plates  with  1.000-inch 
diameter  bore  were  installed  in  both  supply  lines.  These 
orifice  plates  have  a  pressure  drop  of  100  inches  of 
water  for  a  flow  rate  of  1.09  lb/sec  of  dry  air  at  the 
design  conditions.  The  actual  flow  rate  in  each  supply 
line  is  calculated  using  the  relationship: 

.  fap  (in  h20)  P (psia)  500°R' 

m  (lb/sec)  *  1.0907  x -  -  -  (4-13) 

.100  in  H20  600  psia  T(°R)  . 


where  the  pressure  drop  across  each  orifice  plate  (aP)  is 


measured  and  recorded,  the  pressure  (p)  in  both  supply 
lines  is  measured  and  recorded  and  the  temperature  (T)  of 
the  air  is  obtained  from  the  measurement  of  the 
stagnation  temperature  and  pressure  in  each  plenum, 
thermodynamic  properties  of  dry  air  and  the  assumption  of 
adiabatic  flow  from  each  orifice  plate  to  the  plenum 
section.  The  pressure  drop  across  each  orifice  plate  is 
measured  and  recorded  by  two  separate  systems.  Each 
orifice  plate  is  connected  to  a  pneumatic  transmitter 
that  transmits  a  pressure,  proportional  to  the  measured 
pressure  drop,  to  the  strip  chart  recorder.  Each  orifice 
plate  is  also  connected  to  a  strain  gauge-type 
differential  pressure  transducer. 

The  two  types  of  orifice  plates,  choked  and  flow  meter, 
in  each  supply  line  allows  for  comparison.  The  choked 
orifice  plate  fixes  the  mass  flow  rate  for  a  constant 
supply  line  pressure.  However,  due  to  the  characteristics 
of  the  pressure  regulating  valve,  the  supply  line 
pressure  does  vary  slightly  during  a  test  run  and  the 
flow  meter  orifice  plates  measure  the  corresponding 
variation  in  mass  flow  rate. 
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B.  plenum,  Swirl  Generator/Nozzles  and  Test  Section 

The  measurement  system  associated  with  the  concentric 
dual  plenum  section,  the  swirl  generator  and  nozzles 
section  and  the  test  section  is  shown  in  Figure  4-19  for 
tests  without  the  8-inch  outside  wall  and  in  Figure  4-20 
for  tests  with  the  8-inch  outside  wall.  This  system  will 
be  discussed  by  considering  each  of  the  following 
subsystems: 

a.  Scanivalve 

b.  Sensors  and  Signal  Conditioners 

c.  Data  Measurement,  Recording  &  Monitoring 

d.  Controller 

1.  Scanivalve 

The  scanivalve,  shown  in  Figure  4-21,  is  a  Model  SGM 
having  three  valve  assemblies  with  48  pressure 
tabulations  in  each  valve  assembly,  solenoid  drive,  an 
odd-even  tube  marker  position  transmitter  and  pressure 
transducers.  The  scanivalve,  with  its  pressure 
transducers,  was  used  to  measure  all  static  pressures. 
The  stepping  of  the  scanivalve  and  determination  of  its 
position  were  performed  by  the  Controller,  which  will  be 
discussed  later  in  this  section. 
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2.  Sensors  and  signal  conditioners 

A  Model  PM131TC+2.5  Statham  pressure  transducer  was  used 
in  the  Scanivalve  for  all  tests  without  the  8-inch 
outside  wall.  During  these  tests,  this  pressure 
transducer  was  used  to  measure  wall  static  pressures  and 
the  five  pressures  of  probe  #1.  The  pressure  difference 
between  ports  #1  and  #3  of  probe  #1  was  measured  using  a 
Model  PM283TC+0 .15  Statham  pressure  transducer  which  gave 
the  capability  of  approximately  nulling  the  flow  angle 
tor)  of  the  air  stream  relative  to  this  probe. 

Seven  Model  PDCR22+1  ScanCo  differential  pressure 
transducers  were  added  to  the  measurement  system  for 
tests  conducted  with  the  8-inch  outside  wall,  in 
addition,  two  Model  PM80TC+100  Statham  pressure 
transducers  were  installed  for  measuring  the  pressure 
drop  across  the  two  flow  orifice  plates,  one  Model 
PG146TC-5M  Statham  pressure  transducer  was  installed  to 
measure  the  pressure  of  the  air  upstream  of  the  pressure 
regulating  valve  and  one  Model  PG146TC-800  statham 
pressure  transducer  was  installed  to  measure  the  pressure 
of  the  air  downstream  of  this  same  valve.  Two  of  the 
seven  Model  PDCR22+1  ScanCo  differential  pressure 
transducers  were  added  to  the  scanivalve  so  that  all  wall 


static  pressures  were  scanned  on  each  scan.  Five  of  the 
seven  Model  PDCR22+1  Scanco  differential  pressure 
transducers  were  mounted  in  a  special  block  for 
simultaneous  measurement  of  all  pressure  ports  on  the 
five-hole  probe  as  shown  in  Figure  4-22.  The  Model 
PM283TC+0 .15  Statham  pressure  tranducer  was  used  to 
measure  the  difference  between  the  pressures  at  port  #1 
and  port  #3  of  the  five-hole  probe  which  gave  the 
capability  of  approximately  nulling  the  flow  angle  (a,)  of 
the  air  stream  relative  to  the  probe. 

The  air  temperatures  inside  both  the  inner  and  outer 
plenums  were  measured  by  veco  thermistors  manufactured  by 
Victory  Engineering  Corp.  for  tests  without  the  8-inch 
outside  wall  and  were  measured  by  iron/constantan 
thermocouples  for  tests  with  the  8-inch  outside  wall.  The 
angle  of  the  probe  relative  to  the  centerline  of  the 
experimental  apparatus  (  ap  )  was  measured  by  a  10K  Ohm 
precision  potentiometer  that  was  connected  by  gearing  to 
the  shaft  of  the  probe.  The  axial  position  of  the  probe 
for  tests  without  the  8-inch  outside  wall  was  measured  by 
a  Model  WR8-1132/158  position  transducer  manufactured  by 
Lockheed  Electronics  company.  The  radial  position  of  the 
probe  was  measured  by  a  Model  H-3000-S3R  LVDT 
manufactured  by  Columbia  Research. 


All  the  sensors,  except  the  LVDT,  were  connected  using 
the  wiring  configuration  shown  in  Figure  4-23  to  separate 
Model  3511C  Strain  Gage  Power  Supply  and  signal 
Conditioner  units  manufactured  by  Systems  Research 
Corporation.  These  units  are  of  the  "plug-in"  type 
construction  with  designer  cards  for  bridge  completion 
and  can  be  mounted  in  a  rack,  six  of  these  units  were 
used  for  tests  without  the  8-inch  outside  wall  and  14  of 
these  units  were  used  for  tests  with  the  8-inch  outside 
wall.  The  sensors  that  were  not  a  full  bridge  with  a  full 
active  circuit  were  wired  as  shown  in  Figure  4-23  and 
their  bridge  circuit  completed  on  the  designer  card.  The 
type  of  bridge  and  active  circuit  for  the  four  types  of 
sensors  are  listed  below: 


Bridge 

Active 

Type 

Circuit 

Sensor 

pressure  Transducer 

Full 

Full 

Temperature  Thermistor 

Half 

Half 

probe  Angle 

Half 

Full 

Axial  position  Transducer 

Full 

Full 

The  LVDT  was  connected  to  a  Model  300B  Daytronic 
Differential  Transformer  Indicator  which  provided  the 
necessary  electronics  to  give  an  output  signal  that  was 
proportional  to  the  movement  of  the  core  in  the  LVDT. 
Calibration  of  the  LVDT  and  the  probe  Angle  potentiometer 
yielded  nonlinear  characteristics  near  the  ends  of  each's 
range.  Both  calibrations  fit  the  characteristic  equation 
for  the  half  bridge  type  with  full  active  circuit  given 
in  Reference  12  and  rewritten  below  in  the  form  used. 

V  -  V0+  A  <x-x0)/(l+B(x-x0)2  )  (4-14) 

where  v0  is  the  output  voltage,  x  is  the  variable  being 
measured,  A  and  B  are  calibration  coefficients  and  the 
subscript  " o  "  represents  the  center  of  the  range  of 
operation. 

Equation  4-14  was  used  to  represent  the  calibration  of 
all  sensors  that  had  linear  characteristics  over  their 
operating  range  by  setting  B  equal  to  zero.  The  results 
of  sensor/signal  conditioner  calibration  are  listed  below 
using  Equation  4-14. 


3.  Data  Measurement,  Recording  and  Monitoring 

The  output  of  all  signal  conditioners  and  the  Scanivalve 
position  voltage  from  the  controller  were  scanned  and 
measured  by  a  Model  2200B  Datalogger  manufactured  by  the 
John  Fluke  Co.  This  instrument  contains  an  A-D  converter 
that  measures  DC  voltages  on  four  ranges:  40mv,  400  mv, 
4V  and  40V.  Resolution  on  the  40mv  range  is  1  uv.  A 
measurement  speed  of  15  readings  per  second  was  selected 


for  all  test  runs.  The  Datalogger  was  equipped  with  the 
RS-232-C  Interface  option  and  outputted  the  measured  data 
at  a  baud  rate  of  4800  for  recording  on  an  external  tape 
recorder.  The  Datalogger  was  set  up  tp  scan  and  measure 
all  inputs  on  receipt  of  a  scan  signal  from  the 
Controller. 

A  Columbia  Data  Corporation  Model  300C  Data  Cartridge 
Tape  Recorder  was  connected  via  its  MODEM/CPU  cable  to 
the  Datalogger  for  receipt  and  recording  of  data  at  a 
baud  rate  of  4800.  The  terminal  port  of  this  tape 
recorder  was  connected  to  a  Tele-video  Model  920C  CRT 
terminal  for  monitoring  and  review  of  data  during  test 
runs. 

A  Preston  Model  723  Digital  Volt  Meter  (DVM)  was 
connected  to  specific  channels  for  monitoring  via  a 
selector  switch.  This  DVM  was  used  to  monitor  the 
following  channels: 

a.  Differential  pressure  between  ports  fl  and  #3  of 
the  five-hole  probe  during  nulling  of  these  two  pressures 
by  manually  rotating  the  probe. 

b.  Radial  position  of  the  five-hole  probe  as 
Indicated  by  the  output  of  the  LVDT  system  during  manual 
positioning  of  the  probe. 


c.  Axial  position  of  probe  transversing  mechanism 
used  with  probe  #l  during  manual  positioning  of  the 
probe's  axial  position. 

d.  Supply  line  pressure  during  the  manual  setting  of 
the  pilot  pressure  to  the  pressure  regulating  valve  shown 
in  Figure  4-18. 

4.  Controller 

The  Controller  was  designed  using  off-the-shelf 
integrated  circuits  (IC's)  and  References  13  through  16 
to  perform  the  following  functions: 

a.  step  the  Scanivalve  from  position  #l  through  #48 
at  a  preset  rate  after  receipt  of  a  start  command. 

b.  initiate  scan  of  Datalogger  for  each  step  of 
Scanivalve  after  a  preset  time  delay. 

c.  Count  steps  and  output  both  an  analog  voltage 
signal  proportional  to  the  scanivalve  position  and  a 
digital  display  of  the  step  count. 

The  Controller  was  constructed  in  two  physical  parts.  The 
first  part  houses  all  the  integrated  circuits  and 
requires  an  external  5V  DC  power  supply.  The  second  part 
houses  the  switching  circuit  for  the  high  current  20V  DC 
power  that  is  sent  to  the  stepping  solenoid  of  the 
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Scanivalve.  The  ic's  used  in  the  first  part  of  the 
Controller  are: 


IC# 

Device 

Description 

1,2,4,5,9610 

74LS00 

Quad  NAND  Gate 

3,667 

555  . 

Timer 

8 

74LS161 

4-Bit  up  Counter 

11612 

74LS90 

BCD  (Decade)  Counter 

13614 

74LS75 

4-Bit  Latch 

15616 

74LS47 

BCD  to  7-segment  Decoder/Driver 

17618 

- 

7-segment  LED  Single-Digit  Readout 

19  thru  25 

74LS90 

BCD  (Decade)  Counter 

IC's  #11  through  #18  were  mounted  on  two  Radio  Shack 
Catalog  Number  277-103  Single  Digit  Counter  Project 
Boards  to  give  the  two  digit  counter,  ic's  #19  through 
#25  were  mounted  on  a  portion  of  a  Radio  Shack  Catalog 
Number  277-115  Time  Base  Generator  project  Boards  as 
shown  by  the  top  circuit  board  in  Figure  4-24.  IC's  #1 
through  #10  were  mounted  onto  a  plug-in  experimentor  card 
as  shown  by  the  bottom  circuit  board  in  Figure  4-24. 

A  circuit  diagram  of  the  first  part  of  the  controller  is 
shown  in  Figure  4-25.  Quad  NAND  gates  (74LS00)  were  used 
to  provide  the  basic  digital  logic  switching.  The  pin 


connections  for  the  Quad  NAND  gate,  NAND  gate  symbol  and 
truth  table  are  shown  in  Figure  4-26.  Figure  4-25  shows 
the  23  NAND  gates  used  in  the  circuit.  The  particular 
NAND  gate  used  on  a  74LS00  is  designated  by  the  letter 
{A,  B,  C  or  D)  following  the  IC  number  (1,2, 4, 5, 9  or  10) 
contained  within  the  NAND  symbol  in  Figure  4-25. 

integrated  circuits  #3,  #6  and  #7  are  555  timers;  IC  #3 
is  set  up  to  generate  an  adjustable  clock  pulse.  The  user 
can  select  either  this  internal  clock  pulse  or  one  from 
an  external  source  with  a  selector  switch.  The  selected 
clock  pulse  is  then  inputted  into  a  series  of  seven 
decade  counters  (IC  #19  through  #25)  where  the  clock 
period  can  be  multiplied  by  10,  where  n  is  selectable 
from  zero  to  seven.  This  multiplied  clock  output  is  then 
inputted  into  the  logic  circuit  at  connector  "L"  and  also 
connected  to  a  light-emitting  diode  (LED)  for  the 
operator  to  view.  The  waveform  of  this  multiplied  clock 
pulse  is  shown  at  the  top  of  Figure  4-27. 

The  position  transmitter  of  the  scanivalve  is  connected 
to  pins  D,  E  and  F  in  Figure  4-25  and  during  operation 
gives  the  waveforms  shown  in  the  upper  portion  of  Figure 
4-27.  The  printed  circuit  pattern  of  the  position 
transmitter  module  is  shown  in  Figure  4-21  and  was  used 


with  a  5V  DC  common.  Four  NAND  gates  (1A,  IB,  1C  fit  ID) 
are  used  to  generate  an  output  pulse  which  goes  low  when 
the  Scanivalve  is  between  locations  (see  ic  ID  out 
wave£orm  in  Figure  4-27) .  This  output  is  used  to  trigger 
IC  #6  which  is  connected  as  a  single  shot  with  an 
adjustable  length  of  output  pulse.  The  output  waveform  of 
IC  #6  is  shown  in  Figure  4-27.  IC  #6  is  adjusted  so  that 
its  output  pulse  goes  low  during  the  time  that  the 
multiplied  clock  pulse  input  (pin  L)  is  also  low.  The 
rise  in  the  output  of  IC  #6  is  used  to  turn  off  the  power 
applied  to  the  solenoid  of  the  Scanivalve.  Thus  the 
signal  (IC  9C  out)  to  the  solenoid  is  turned  off  just  as 
soon  as  it  is  known  that  the  Scanivalve  is  at  the  new 
location  and  thus  reduces  both  the  power  consumption  and 
the  chance  of  skipping  a  location.  The  output  of  NAND 
gate  9B  is  used  to  initiate  a  scan  by  the  Datalogger.  The 
scan  signal  starts  at  the  beginning  of  the  second  half  of 
the  multiplied  clock  pulse  (low  signal  input  to  pin  L) 
and  ends  when  the  output  of  the  single  shot  ic  #6  drops 
to  low.  Thus  the  length  of  the  output  pulse  from  ic  #6 
determines  the  length  of  the  scan  signal  and  care  must  be 
taken  so  that  this  signal  is  long  enough  to  initiate  the 
scan  but  not  too  long  so  as  to  initiate  a  second  scan. 

IC  #7  is  connected  as  a  single  shot  that  is  triggered  by 


either  the  output  of  NAND  gate  ID  or  a  manual  count 
switch  depending  on  the  position  of  the  manual-auto 
selector  switch.  This  single  shot  gives  a  clear  count 
pulse  to  the  two-digit  counter  (iC's  #11  through  #18)  and 
the  voltage  staircase  (IC  #8) .  The  output  of  the  voltage 
staircase  versus  the  Scanivalve  position  is  shown  in 
Figure  4-28. 

The  zero  (reset  or  stop)  signal  from  the  scanivalve 
position  transmitter  is  connected  to  pin  E  in  Figure 
4-25.  When  pin  E  is  high  (zero  signal),  then  the  signal 
at  pin  K  remains  high  and  the  signal  at  pin  R  remains 
low.  A  cycle  of  the  stepping  of  the  Scanivalve  is 
initiated  by  bringing  the  reset  signal  at  pin  j  to  the 
low  state  and  the  waveforms  as  shown  in  Figure  4-27  will 
then  occur. 

The  circuit  of  the  second  portion  of  the  Controller  is 
shown  in  Figure  4-29.  The  output  signal  from  the  first 
portion  of  the  Controller  is  used  to  turn  on  the  RS-2039 
transistor  which  then  energizes  the  relay.  The  scanivalve 
can  be  stepped  by  either  pressing  the  step  button  when 
the  selector  switch  is  in  the  manual  position  or  by 
energizing  the  relay  when  the  selector  switch  is  in  the 
auto  position.  The  Scanivalve  can  be  reset  or  homed  to 


the  zero  or  48  position  by  pressing  the  reset  button  and 
holding  until  the  unit  has  stopped. 
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CHAPTER  5 


FIVE-HOLE  PROBES 

I.  INTRODUCTION 

A  five-hole  probe  is  essential  to  the  measurement  of  mean 
flow  in  the  incompressible  mixing  of  concentric  jets  with 
swirl  in  the  inner  stream.  This  type  of  probe  is 
constructed  with  a  single  forward  hole  surrounded  by  four 
holes  on  a  conical  surface  eminating  from  the  forward 
hole.  Erwin  (Reference  1)  describes  the  general 
characteristics  of  this  type  of  probe  and  gives  typical 
calibration  data  for  a  five-hole  probe.  The  procedure 
described  by  Erwin  to  use  this  type  of  probe  requires 
that  the  pressures  from  two  holes  on  opposite  sides  of 

the  conical  surface  be  nulled,  first  for  determination  of 

one  flow  angle,  and  then  the  other  flow  angle  is 
determined  by  the  other  three  pressures  and  the 
calibration  curve.  Baker,  Gallington  and  Minster 
(Reference  2)  describe  a  calibration  procedure  for  this 
type  of  probe  which  does  not  require  that  the  probe  be 
moved  to  a  null  position  during  measurement.  The  probe 
described  by  Erwin  is  not  practical  for  use  in  a 

blow-down  system  and  the  location  of  the  probe  head 
changes  when  the  probe  is  rotated.  The  calibration 

procedure  described  by  Baker  et  al  alleviates 
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these  drawbacks,  however,  the  accuracy  o£  measurement  is 
degraded  due  to  the  larger  calibration  regime.  Dring, 
joslyn  and  Hardin  (Reference  3)  describe  the  use  of  a 
modified  five-hole  probe  that  keeps  the  probe  head 
essentially  in  a  constant  position  when  the  probe  is 
rotated.  However,  this  modified  probe,  as  calibrated  by 
Dring  et  al,  requires  that  two  pressures  be  nulled  during 
measurement.  The  required  high  accuracy  and  the  limited 
time  of  measurement  are  satisfied  by  the  use  of  a  probe 
as  described  by  Drln^  et  al,  but  using  a  modification  of 
the  calibration  and  measurement  procedures  described  by 
Baker  et  al. 

Several  small  probes  were  constructed  using  the  methods 
described  by  Gallington  and  Hollenbaugh  (Reference  4). 
These  probes  were  calibrated  using  the  methods  of  Baker 
et  al  over  a  small  range  of  angle  of  attack  (a)  and  a 
large  range  of  angle  of  sideslip  (p) .  During  measurement, 
the  probe  is  rotated  to  minimize  the  pressure  difference 
between  the  two  pressure  holes  on  the  conical  surface 
that  are  directly  related  to  the  angle  of  attack  (a). 
Since  nulling  is  not  required  for  determining  the  exact 
angle  of  attack,  the  time  for  measurement  of  the  flow  by 
the  five-hole  probe  is  reduced  considerably.  The  accuracy 
of  measurement  still  remains  high  since  the  probe  is 
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calibrated  for  a  small  difference  in  the  flow's  angle  of 
attack  relative  to  the  probe's  axis  and  for  a  large  range 
of  angle  of  sideslip. 

II.  THEORETICAL  SECTION 

A.  Calibration  Equations 

The  coordinate  system  for  the  five-hole  conical  probe  is 
shown  in  Figure  5-1  and  the  port  numbering  system  is 
shown  in  Figure  5-2.  Figure  5-1  shows  the  end  of  the 
pressure  probe  and  only  the  forward  hole  -  port  5.  The 
velocity  components  in  the  radial,  tangential  and  axial 


directions 

(u. 

V 

and  w)  can  be  written  in 

terms  of 

the 

magnitude 

of 

the 

velocity,  v,  and  the  flow 

angles, a 

and 

0,  as  follows: 

u 

=  V 

sin0 

(5-1) 

V 

*  V 

sin  a 

COS  0 

(5-2) 

w 

»  V 

cos  a 

COS  0 

(5-3) 

As  recommended  by  Baker,  Gallington  and  Minster 
(reference  2),  four  pressure  coefficients  are  defined  in 
terms  of  the  pressures  measured  from  each  of  the  five 
pressure  ports  on  the  probe  and  the  measured  total  and 


static  pressures 


where  P1  through  pg  are  the  pressures  measured  at 
pressure  ports  1  through  5  shown  in  Figure  5-2;  P£  is  the 
sum  of  p1  ,  p2  ,  P3  and  P4 ;  PT  is  the  total  pressure;  and  pg 
is  the  static  pressure. 


For  a  fixed  flow  angle  8  and  small  changes  in  the  flow 
angle  a  ,  the  coefficient  Ca  is  roughly  proportional  to 
v/W.  Similarly,  for  small  changes  in  the  flow  angle  8  , 
the  coefficient  Cp  is  roughly  proportional  to  u/v.  The 
coefficient  cM  is  the  difference  between  the  actual  total 
pressure  and  the  pressure  measured  by  port  5 
nondimensionalized  by  the  difference  between  P5  and  1/4 
p£ .  The  coefficient  cq  is  a  ratio  of  that  portion  of  the 
dynamic  pressure  sensed  by  the  probe,  p#  -  1/4  P£  ,  to 
the  actual  dynamic  pressure,  Pr  -  P$  . 
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For  calibration  of  a  particular  probe,  the  following 
fourth-order  equation  for  C0  ,  Cq  ,  u/V  and  v/v  is  used 
as  recommended  by  Barker  et  al  (Reference  2) : 

k*.cI'-k»,cJ*k«cJcs*k.c0c|+ 

^  (5-8) 

KA,,<V3 

where  A  -  C0  ,  Cq  ,  u/v,  v/v  as  applicable  and  the  15 
subscripted  K's  are  calibration  coefficients  to  be 
determined  for  each  of  the  four  A's.  From  a  large  number 
of  data  points,  the  calibration  coefficients  can  be 
determined  by  using  matrix  operations  as  discussed  by 
Neter  and  Wasserman  (Reference  5)  and  the  resulting 
calibration  automatically  incorporates  a  least  squares 
fit.  Nine  matrices  are  required  in  this  calibration 
process: 

1  ^  Cp1  S  Cf1  '  ’ 

!  SS  S  SS  : : 

<4.  Cfc  <4  : ;  <*J 


I 


(5-9) 


where  the  subscripts  on  Ca  and  Cp  denote  the  particular 
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By  use  of  the  matrix  notation.  Equation  5-8  can  be 
written  for  all  data  points  as  follows: 

[a],  [d]  [k4]  (5-12) 

where  A  =*  u/v  ,  v/v  ,  C0  /  Cq  and 

Ka  »  Ku/y  ,  Kv/V  ,  Kc<j  ,  Kg  as  applicable 

Solving  Equation  5-12  for  the  calibration  coefficients 

that  give  the  least  square  error  between  the  actual  data 
and  the  value  predicted  by  the  calibration  can  be  done  as 
follows  for  the  u/v  matrix: 

1.  Write  Equation  5-12  in  terms  of  u/V 

tu/V]  -  CO]  [Ku,v]  (5-13) 

2.  Multiply  each  side  of  Equation  5-13  by  the 

transpose  of  the  [d]  matrix.  Equation  5-9 

C“/V][D]t  -CD]t  CD]  [Ku/V]  (5-14) 

3.  Multiply  each  side  of  Equation  5-14  by  the 

inverse  of  the  newly  formed  [dtd]  matrix 

cdtqi>/v]  cdt]  -  cdtd]4cdtd)  [ku/v]  (5-1 5) 

4.  Noting  that  the  product  of  a  matrix  and  its 


inverse  is  the  identity  matrix.  Equation  5-15  becomes 

[Ku  v]  -  CDt  D]'1  [u/v]  [DT]  (5-16) 
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The  same  derivation  is  followed  to  derive  the  following 
equations  for  the  calibration  coefficients  K,  ,  Kr 

v  V  >*0 

and  Kr  : 
co 

[K./J-CDTd]'1  [v/v]  [Dt]  (5-17) 

[KCJ  -CDtD]'1  [Co]  [Dt]  (5-18) 

[KCJ-CDt0]4[C<1]  [Dt]  (5.19) 


Thus,  from  a  given  set  of  data,  the  D  ,  u/v 
and  cq  matricies  are  calculated  and, 

operations,  the  calibration  coefficients 

Kc  ,and  kc  are  calculated,  using  Equations 

''O  q 

5-18  and  5-19,  respectively. 


,  v/v  ,  C0 , 
by  matrix 

Ku/V  '  ^v/v  > 
5-16,  5-17, 


B.  prediction  of  Static  and  Dynamic  pressures 

Once  the  calibration  coefficients  have  been  determined, 
then  the  static  pressure,  P$  ,  and  the  dynamic  pressure, 
PT-PS ,  can  be  predicted  from  the  five  pressures  measured 
by  the  probe,  P1  through  P5  ,  and  Equations  5-4  through 
5-8  as  follows: 

1.  From  the  pressures  p.  through  p  ,  use  Equations 
5-4  and  5-5  to  calculate  Ca  and  Cp  . 


'  ■ 
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2.  using  the  calculated  Ca  and  cp  and  the 

calibration  coefficients  Kr  and  Kr  ,  calculate  the 

v«0  '-q 

predicted  value  of  C0  and  cq  ,  using  Equation  5-8. 

3.  From  the  predicted  value  of  Cq  and  the  five 
measured  pressures,  p1  through  pg  ,  calculate  the  dymanic 
pressure,  p  -p  ,  from  Equation  5-7,  written  as  follows: 

T  O 

PT_PS  =  (PS"  1/4  Pr  )  /Cq  (5-20) 

4.  From  the  predicted  values  of  Cq ,  C0  and  PT -Ps 
and  the  five  measured  pressures,  P1  through  P5  ,  calculate 
the  static  pressure,  Ps  ,  from  the  combination  of 
Equations  5-6  and  5-7,  written  in  the  following  form: 

Ps  *  P5  -  (1/Cq-C0)  (P5-  1/4  P£  )  (5-21) 

The  predicted  value  of  the  total  pressure,  PT  ,  can  also 
be  calculated  from  Equation  5*6,  written  as  follows: 

PT=  P5  +  C0  (P5  -  1/4PZ  )  (5-22) 

C.  Prediction  of  Momentum 

in  axisymmetric  flow,  the  tangential  and  axial  momentum 
fluxes  are  quantities  of  interest.  By  aligning  the  x-axis 
of  the  five-hole  probe  with  the  radial  direction  of  flow 
and  the  z-axis  of  the  probe  with  the  axial  direction  of 
flow,  this  probe  can  be  used  to  determine  the  radial, 


132 


tangential  and  axial  velocities  and  the  static  pressure 
at  a  point.  Figure  5-3a  shows  this  type  of  alignment. 

These  probes  will  be  used  in  measurement  where  air 
streams  are  provided  by  a  blow-down  system.  For  accurate 
determination  of  the  flow  direction,  the  probe  will  be 
rotated  about  its  x-axis  to  essentially  null  the  flow 
angle  a  as  shown  in  Figure  5-3b.  when  the  probe  is 
rotated,  as  shown  in  Figure  5-3b,  the  u,  v  and  w  velocity 
components  can  be  determined  with  respect  to  the  probe's 
body  fixed  coordinate  system  (x,y',z')  and  then  the 
rotation  of  the  probe  can  be  taken  into  account  in 
determining  the  values  of  the  v  and  w  velocity  components 
with  respect  to  the  x,y,z  coordinate  system.  The 
resulting  u,  v  and  w  velocity  components  are  then  the 
radial,  tangential  and  axial  velocity  components, 
respectively,  in  the  cylindrical  coordinate  system  (r^ 
,z),  shown  in  both  Figures  5-3a  and  5-3b. 

The  axial  and  tangential  momentum  fluxes  can  be  written 
as  follows  for  incompressible  axisymmetric  flow: 


- 

Ps-  P 

+  p  w2 

(5-23) 

= 

o  (w2- 

v2/2) 

(5-24) 

m0 

= 

0  vw 

(5-25) 

where  M2 

and 

Mi  are 

the  axial  momentum  fluxes, 

Me  is 
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tangential  momentum  flux,  p0  is  a  reference  pressure,  and 
is  the  density  of  the  incompressible  fluid. 


Since  the  dynamic  pressure,  PT -Ps  ,  equals  pv  /2  for  the 
case  of  incompressible  flow,  then  the  equations  for  the 
axial  and  tangential  momentum  fluxes  can  be  written  in 
terms  of  the  values  predicted,  using  the  five-hole  probe 
as  follows: 

Mz  =  Ps-  P0  +  2(Pt  -  Ps  )  (w/v)2  (5-26) 

Mj  =  2  (PT  -  Ps  )  ((W/V)  -  (v/V)  /2)  (5-27) 

Me  =  2  (PT  -  Ps  )  (v/V)  (w/V)  (5-28) 

where  w/V  *  (l-(u/V)2  -(v/v)2  )Vz 


Thus  Equations  5-26,  5-27,  and  5-28  can  be  used  to 

determine  M2 ,  H' ,  and  M0  from  the  redicted  values  of  pg  , 
PT-PS,  u/v  and  v/v. 


D.  Measure  of  Accuracy 


The  calibration  procedure  consists  of  measuring  the  seven 
pressures  (P1  ,  P2  ,  P3  ,  P4  ,  P5  ,  PT  &  Ps )  and  the  two 

angles  (  a  and  B  )  at  each  point  of  interset.  From  this 
given  set  of  data,  the  D  ,  u/v  ,  v/v  ,  C0  and  Cq . 
matrices  are  calculated  using  Equations  5-16,  5-17,  5-18 
and  5-19.  Once  the  four  bets  of  calibration  coefficients 
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are  determined,  the  values  of  u/v,  v/v,  0o  and  Cq  can  be 
calculated  using  Equation  5-8  and  compared  to  the 
experimental  values  of  these  quantities  at  each  data 
point.  This  comparison  of  calculated  and  experimental 
values  of  u/v,  v/v,  CQ  and  Cq  constitute  a  validity  test 
of  the  calibration. 

as  discussed  by  Barker,  Gallington  and  Minster,  the 
accuracy,  as  defined  by  the  validity  test  of  the 
calibration,  is  dependent  on  the  way  the  calibration  data 
samples  are  taken,  the  manufacturing  anomalies  of  the 
probe,  and  the  associated  measuring  equipment  used  with 
the  probe,  since  the  calibration  procedure  utilizes  a 
curve  fitting  routine,  the  greater  accuracy  is  afforded 
in  the  region  where  the  calibration  data  samples  have  the 
higher  density.  Thus,  the  density  of  the  data  points 
effects  the  accuracy  and  the  distribution  of  the  data 
points  should  be  selected  in  such  a  way  as  to  insure  that 
the  accuracy  is  constant  over  the  probe's  range  of  use. 

When  calibrating  a  five-hole  probe,  the  range  of  flow 
angles,  ct  and  8  ,  to  be  measured  should  be  known 
beforehand.  The  range  and  density  of  the  calibration 
points  can  then  be  selected  to  insure  a  high  accuracy 
where  desired. 
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The  standard  deviation,  S,  of  the  calculated  values  of 
the  four  desired  outputs  (u/v,  v/v,  CQ  and  Cq  )  obtained 
from  the  experimental  values  of  these  quantities  gives  an 
overall  measure  of  the  accuracy  of  the  calibration.  This 
standard  deviation  can  be  estimated  using  the  following 
equation: 


where  A  *  u/v,  v/v,  CQ  and  cq  in  turn,  n  is  the  number  of 
data  points  being  compared,  and  e  and  c  are  subscripts 
that  refer  to  experimental  and  calculated  values, 
respectively. 

Assuming  a  normal  distribution  of  errors,  one  then  finds 
a  95%  probability  that  the  actual  value  will  be  within 
1.96  standard  deviations  of  that  calculated  from  the 
curve  fit. 

An  important  point  to  note  is  that  the  calibration 
process  for  the  five-hole  probe  is  limited  to 
incompressible  flow  because  the  calibration  depends  on 
Ca  and  Cp  ,  which  are  pressure  coefficients  that  are 
dependent  on  Mach  number. 
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III.  EXPERIMENTAL  SECTION 

A.  Apparatus  -  probe  #1 

The  calibration  experiment  was  conducted  in  the  30-inch 
hexagonal  subsonic  venturi  wind  tunnel  at  the  university 
of  Washington.  The  five-hole  probe  to  be  tested  (probe 
#1)  was  mounted  to  its  traversing  mechanism  which  allowed 
variations  in  the  flow  angle  a  and  radial  position.  This 
probe/transversing  mechanism  assembly  was  attached  to  a 
calibration  fixture  which  allowed  variations  in  the  flow 
angle  0  .  Both  flow  angles,  a  and  S  ,  were  changed 
manually,  pressures  P1  ,  p2  ,  p3  ,  P4  ,  P5  ,  pt  and  pg  were 
measured  using  a  Scanivalve  with  a  Statham  pressure 
transducer  and  a  digital  volt  meter  (DVM) ,  as  shown  in 
Figure  5-4a.  The  positioning  of  probe  #1,  in  the  tunnel 
for  calibration,  is  shown  in  Figure  5-4b. 

Figure  5-5  shows  the  shape  and  dimensions  of  Probe  #l. 
This  probe  has  a  hook  shape  so  that  rotation  of  the 
probe's  mounting  shaft  (a  part  of  the  probe/transversing 
mechanism)  only  rotates  the  probe's  face  and  does  not 
change  the  location  of  the  probe's  face.  The  probe  has  an 
outside  diameter  of  0.158  inches  with  each  of  the  five 
tubes  having  an  inside  diameter  of  0.021  inches.  The 


probe  tip  has  a  conical  half-angle  of  45  degrees. 


B.  initial  calibration  -  probe  #l 


The  probe  was  mounted  to  its  traversing  mechanism, 
leveled  on  a  surface  plate  and  the  resistance  of  the 
potentiometer  that  measures  the  angle  of  the  probe,  with 
respect  to  the  traversing  mechanism,  was  measured  to 
establish  a  zero  for  the  flow  angle  a  .  The  calibration 
fixture  was  mounted  to  the  wind  tunnel  and  aligned  with 
the  flow.  The  probe/traversing  mechanism,  with  the  probe 
attached,  was  mounted  to  the  calibration  fixture  for  a 
flow  angle  6  of  zero  and  the  probe  position  adjusted  for 
a  flow  angle  of  zero.  The  flow  angle  a  of  the  probe 
versus  the  output  voltage  of  a  bridge  circuit  that 
included  the  potentiometer  connected  to  the  shaft  of  the 
traversing  mechanism  was  used  to  set  the  flow  angle  a  . 
For  each  flow  angle  B  ,  the  flow  angle  a  was  varied 
manually  and  the  seven  pressures  and  temperatures 
recorded.  The  position  of  the  traversing  mechanism  on  the 
calibration  fixture  was  changed  manually  to  a  new  flow 
angle  6  and  the  measurements  repeated  over  the  range  of 
flow  angle  a .  The  wind  tunnel  was  operated  at  a  velocity 
of  approximately  82  ft/sec  during  this  phase  of  the 
calibration.  Data  was  taken  for  the  following  flow 
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angles: 


0,±1 

.9  ,±3.8 

,±5. 

7  ,±7.5  ,+11.3 

±15. 

1  ,±18.9 

,+23 

.6  ,+28.3  ,±32 

±36. 

2  ,±40.2 

,±48 

&  +51.5 

0  ,±2 

,±4  ,+6 

#±8 

,+10  ,+15  ,+20 

+25 

,+30  &  ■ 

*35 

This  phase  of  the  calibration  was  performed  over  a  larger 
range  of  flow  angle  a  than  the  actual  probe  will 
encounter,  so  that  the  complete  characteristics  of  this 
probe  could  be  determined  and  compared  with  other  probes 
whose  characteristics  are  documented.  The  pressure 
coefficients  CQ  ,  Cq ,  Ca  and  Cp  are  plotted  versus  the 
flow  angle  a  for  values  of  the  flow  angle  8  of  -10°, 
0°,  +10° ,  +25°  and  +35°  in  Figures  5-6,  5-7,  5-8,  5-9  and 
5-10,  respectively.  Note  the  near  linear  characteristic 
of  the  four  pressure  coefficients  versus  a  in  the  range 
of  a  from  -7.5°  to  +7.5°.  The  pressure  coefficients 
C0,  Cq  ,  Ca  and  Cp  are  plotted  versus  the  flow  angleS 
for  a»0  in  Figure  5-11  and  the  need  for  two  distinct 
calibration  regions  is  shown  as  the  characteristics  of 
the  pressure  coefficients  C0  ,  Cq  and  Cp  change  from 
nearly  constant  for  values  of  8  between  -10°  and  +10° to 
significant  changes  for  values  of  6  greater  than  4-10°. 
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The  calibration  data  was  separated  into  two  regions  for 
calibration.  The  first  calibration  region  being  for 
values  of  8  between  -10  and  +10  and  the  second 
calibration  region  being  for  values  of  8  greater  than 
+10  .  Since  probe  #1  would  be  subjected  to  only  small 
values  of  the  flow  angle  a  ,  the  probe  was  calibrated 
using  data  which  extended  from  a  equal  to  -7.5  to  +7.5  . 
The  first  calibration  region  has  99  data  points  (±<x  =»  0 , 
1.9  ,  3.8  ,  5.7  &  7.5  and  +8  -  0 ,  2  ,  4  ,  6  ,  8  &  10  ) 
and  the  second  calibration  region  has  54  data  points 
Cta  *  0  ,  1.9  ,  3.8  ,  5.7  &  7.5  and  8*  10  ,  15  ,  20  , 
25  ,  30  &  35  ) .  Both  calibration  regions  have  the  9  data 
points  at  the  flow  angle  of  8 *+10  in  common. 

C.  Calibration  Coefficients  -  probe  #1 


For  the  calibration  regions,  the  appropriate  set  of  data 

was  read  into  the  university  of  Washington's  CDC  computer 

for  reduction  into  the  five  matrices  J^oj  ,  j^u/vj  »  |v/vj  ,  |c0j. 

and  j^Cq  j  using  Equations  5-1  through  5-10.  Then  the 

calibration  coefficients  K,.,„  ,  Kw/„  ,  Kr  and  Kr  of 

u/v  *  */V  »*o  cq 

probe  #1  were  calculated  for  each  region  using  Equations 


5-16  through  5-19.  The  computer  program  for  calibration 
is  written  in  Fortran  and  includes  call  to  the  linear 
algebra  library  package  LINPACK  for  solution  of  Equations 


5-16  through  5-19.  Equation  9  was  used  to  calculate  the 
predicted  values  of  u/V,  v/v,  C0  and  Cq  ,  using  the 
measured  values  of  ca  &  Cp  and  the  specific  calibration 
coefficients  for  that  region  for  comparison  with  the 
calibration  values  of  these  variables. 

The  results  of  the  calibration  are  presented  in  Figures 
5-12  through  5-18.  Lines  of  constant  Ca  and  Cp  are 
plotted  on  a  u/V  versus  v/v  plot  in  Figures  5-12  and  5-13 
for  the  first  and  second  calibration  regions, 
respectively.  For  the  first  calibration  region,  constant 
contours  of  C0  plotted  against  Ca  and  Cg  are  shown  in 
Figure  5-14  and  constant  contours  of  Cq  plotted  against 
Ca  and  Cp  are  shown  in  Figure  5-16.  Contours  of  C„  and 
Cq  for  the  second  calibration  region  are  plotted  in 
Figures  5-15  and  5-17,  respectively.  The  value  of  the 
pressure  coefficient  Cp  is  used  to  determine  which  set 
of  calibration  coefficients  are  to  be  used  to  calculate 
the  flow  properties.  The  first  region's  calibration 
coefficients  are  used  for  values  of  Cp  less  than  0.35 
and  the  second  region's  calibration  coefficients  are  used 
for  values  of  Cp  greater  than  or  equal  to  0.35.  Figure 
5-18  provides  a  comparison  between  the  actual  and 
predicted  values  of  u/V  and  v/v.  The  difference  between 
the  two  values  represents  the  error  in  determining  the 
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flow  direction. 

D.  Apparatus  -  probe  #2 

The  five-hole  probe  for  measuring  the  mixing  with  swirl 
between  and  eight-inch  outer  tube  and  a  four-inch  inner 
tube  is  shown  in  Figure  5-19.  This  probe  also  has  a  hook 
shape  so  that  rotation  of  the  probe's  shaft  only  rotates 
the  face  of  the  probe  and  doesn't  change  the  location  of 
the  probe's  face.  The  probe  has  an  outside  diameter  of 
0.158  inches  with  each  of  the  five  tubes  having  an  inside 
diameter  of  0.021  inches.  The  probe  tip  has  a  conical 
half-angle  of  45  . 

probe  #2  and  its  tranversing  mechanism  was  removed  from 
its  section  of  eight-inch  tubing  and  mounted  to  a 
calibration  fixture  which  allowed  variations  in  the  flow 
angle  B  .  The  transversing  mechanism  for  this  probe 
allowed  for  radial  transversing  and  variation  of  the  flow 
angle  a  .  This  transversing  mechanism  is  shown  in  Figure 
5-20  and  its  potentiometer  was  used  to  determine  the  flow 
angle  .  Both  flow  angles, a  and  B  ,  were  changed  manually, 
pressures  p,  ,  P2  ,  P3  •  P4  •  ps  •  pt  and  ps  were  measured 
using  seven  separate  scanco  pressure  transducers.  The 
measurement  system  for  the  calibration  of  probe  #2  is 


shown  in  Figure  5-21a.  The  positioning  of  probe  #2  for 
calibration  in  the  30-inch  hexagonal  subsonic  venturi 
wind  tunnel  at  the  University  of  Washington  is  shown  in 
Figure  5-21b. 

E.  initial  Calibration  -  probe  #2 

probe  #2  was  manually  positioned  at  different  values  of 
the  flow  angles  a  and  S  in  the  same  manner  as  described 
for  probe  #1  in  Section  B.  The  wind  tunnel  was  operated 
at  a  velocity  of  approximately  75  ft/sec  during  this 
phase  of  calibration.  Data  was  taken  for  the  following 
flow  angles: 

a  =  0 ,+5 ,+10 ,±15, ±20 ,±25 ,+30  ,+35, ±40  , 

±45  &  +50 

8—  0  ,±2  ,±4  ,±6  ,±8  ,±10  ,  +  15  &  +20 
This  phase  of  calibration  was  performed  over  a  larger 
range  of  flow  angles  a  and  8  than  the  actual  probe  will 
encounter,  so  that  the  complete  characteristics  of  this 
probe  could  be  determined  and  compared  with  probe  #1  and 
other  probes  whose  characteristics  are  documented.  The 
pressure  coefficients  Ca  ,  Cp  ,  C0  and  cq  are  plotted 
versus  the  flow  angle  a  for  values  of  the  flow  angle B= 
-10,0, +10  and  +20  in  Figures  5-22,5-23,5-24  and  5-25, 
respectively,  pressure  coefficients  ca  ,  Cp  and  cq  are 
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nearly  linear  versus  the  flow  angle  a  in  the  range  of  a 
from  -20  to  +20.  The  pressure  coefficient  C0  has  a 
minimum  near  the  flow  angle  a  value  of  +5,  indicating  a 
degree  of  asymmetry  in  the  construction  of  probe  #2. 

probe  #2  was  to  be  used  in  a  diverse  environment  which 
would  extend  over  a  large  range  in  the  value  of  the  flow 
angle  a  .  ideally,  the  probe  would  be  adjusted  to  null 
pressures  p  and  p  ,  thus  giving  a  C_  value  of  zero. 
However,  the  probe  was  also  to  be  subjected  to  large 
values  of  the  flow  angle  a  that  could  be  aslarge  as  -50. 
The  four  pressure  coefficients  are  plotted  versus  the 
flow  angle  a  for  values  of  a  from  -50  to  +50  for  a 
value  of  the  flow  angle  8  equal  to  +4  in  Figure  5-26. 
The  pressure  coefficients  C„  and  Cp  are  nearly  linear 
versus  a  in  the  range  of  from  -50  to  +20.  However,  the 
pressure  coefficients  C0  and  Cq  vary  nonlinearly  with 
the  flow  angle  a  .  Two  calibration  regions  were  selected 
to  meet  the  operating  range  of  Probe  #2. 


Calibration 

regions: 

Reg  ion 
1 
2 


data  were  separated  into  the  following  two 

Range  of  Flow  Angles 
-20:50^  +  20  ,-10^6^+10 
-50^oi^-20,-  6^S^+  8 


The  range  of  the  flow  angle  8  for  calibration  of  the 
second  region  is  smaller  than  that  for  the  first  region 
because  the  probe  would  be  operating  in  the  second  region 
when  measuring  confined  flows  that  had  small  swirl  angle 
and  a  resultant  small  radial  angle  (  8  ) •  The  first 
calibration  region  has  99  data  points  (  o=0 ,+5 ,+10 ,±15  & 
_20  and  8=0 ,+2 ,+4 ,+6 ,+8  &_10)  and  the  second  calibration 
region  has  56  data  points  (  a =-20  ,-25  ,-30  ,-35  ,-40  ,-45  & 
-50  and  8=0, ±2, +4 ,+6  &  +8).  Both  calibration  regions 

have  7  data  points  at  the  flow  angle  of  ot=  -20°  in 
common. 

F.  Calibration  Coefficients  -  probe  #2 

The  calibration  coefficients  K,.,  ,  K„.  ,  K~  and  k~ 

u/v  w^v  c0  Cq 

of  probe  #2  were  calculated  for  each  region  using  the 
same  method  as  described  previously  for  probe  #1.  The 
predicted  values  of  u/v,  v/vr  C0  and  Cq  were  calculated 
using  Equation  5-8,  the  measured  values  of  Ca  &  Cp  and 
the  specific  calibration  coefficients  for  that  region  so 
that  these  values  could  be  compared  with  the  calibration 
values  of  these  variables. 

The  results  of  the  calibration  are  presented  in  Figures 
5-27  through  5-33.  Lines  of  constant  CQ  and  Cp  are 


plotted  on  a  u/v  versus  v/v  plot  in  Figures  5-27  and  5-28 
for  the  first  and  second  calibration  regions, 
respectively.  For  the  first  calibration  region,  constant 
contours  of  C0  plotted  against  CQ  and  Cp  are  presented 
in  Figure  5-29  and  constant  contours  of  Cq  plotted 
against  Ca  and  Cp  presented  in  Figure  5-31.  Contours  of 
C0  and  Cq  for  the  second  calibration  region  are  plotted 
in  Figures  5-30  and  5-32,  respectively.  The  value  of  the 
pressure  coefficient  Ca  is  used  to  determine  which  set 
of  calibration  coefficients  are  to  be  used  to  calculate 
the  flow  properties.  The  first  region's  calibration 
coefficients  are  used  for  values  of  Ca  greater  than  or 
equal  to  -0.65  and  the  second  region's  calibration 
coefficients  are  used  for  values  of  Ca  less  than  -0.65. 
Figure  5-33  gives  a  comparison  between  the  actual  and 
predicted  values  of  u/v  and  v/v. 

IV.  DISCUSSION 

It  is  convenient  to  separate  the  discussion  of  each 
probe's  behavior  into  sub-sections  regarding  flow  angles, 
total  and  dynamic  pressures  and  accuracy. 


A.  Flow  Angles  -  probe  #1 


A  comparison  between  the  actual  values  of  u/v  and  v/v  and 
the  values  predicted  using  the  polynominals  of  Equation 
5-8  are  given  in  Figure  5-18.  The  error  in  determining 
the  flow  angle  is  indicated  by  the  difference  between 
each  set  of  predicted  and  actual  values  of  u/v  and  v/v. 
This  figure  shows  that  the  fit  is  most  accurate  near  the 
origin  of  the  tangential  velocity  ratio,  v/v  (flow  angle 
«=0).  This  is  not  surprising  since  the  density  of  data 
points  was  greatest  in  this  region,  thus  forcing  a  higher 
accuracy  in  this  region.  This  figure  also  shows  that  the 
error  varies  in  an  unpredictable  way.  Both  of  the  above 
observations  concerning  the  comparison  between  actual  and 
predicted  values  of  u/v  and  v/v  were  also  noted  by  Baker, 
et  al,  in  Reference  2.  Figure  5-18  also  shows  that  the 
error  is  greater  in  the  second  calibration  region  (values 
of  the  radial  velocity  ratio,  u/v,  greater  than  0.15) 
than  in  the  first  calibration  region.  This  difference  is 
due  mainly  to  the  lower  density  of  data  points  in  the 
second  calibration  region  and  to  the  severe  flow  angle 
8  . 

Figures  5-12  and  5-13  are  plots  of  u/v  versus  v/v  showing 
lines  of  constant  Ca  and  for  the  first  and  second 


calibration  regions,  respectively.  The  lines  of  constant 
Ca  generally  tend  to  be  vertical  and  the  lines  of 
constant  Cp  generally  tend  to  be  horizontal.  It  is 
apparent  from  the  almost  vertical  Ca  lines  and  almost 
horizontal  Cp  lines  that  CQ  is  almost  independent  of  u/V 
and  Cp  almost  independent  of  v/v.  The  location  of  the  Ca 
equal  zero  line  in  both  Figures  5-12  and  5-13  gives  the 
value  of  the  flow  angle  a  at  which  the  pressure 
coefficient  Ca  is  nulled.  The  location  of  the  Cp  equal 
to  zero  line  in  only  Figure  5-12  gives  the  value  of  the 
flow  angle  8  at  which  the  pressure  coefficient  Cp  is 
nulled.  The  location  of  the  Cp  equal  zero  line  in  Figure 
5-13  does  not  give  the  value  of  the  flow  angle  at  which 
the  pressure  coefficient  Cp  is  nulled  because  the  null 
value  of  this  pressure  coefficient  is  outside  the 
calibration  region  and  plotted  from  extrapolation  to  show 
an  overlap  with  the  region  shown  in  Figure  5-12. 


B.  Total  and  Dynamic  pressures  -  probe  #1 


Constant  contours  of  C  is  plotted  against  C0  and  cp  are 
shown  in  Figures  5-14  and  5-15  for  the  first  and  second 
calibration  regions,  respectively.  As  indicated  in 
Equation  5-6,  C0  is  proportional  to  the  difference 
-  •*-.^een  the  actual  total  pressure,  PT  ,  and  the  total 
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pressure  measured  by  pressure  port  5  of  the  probe,  pg  . 

Figure  5-14  shows  that,  for  small  angles,  pg  gives  an 

accurate  reading  of  the  actual  total  pressure.  The  larger 

the  angle,  the  larger  the  difference  between  p  and  p  . 

*  5 

The  constant  contours  of  C0  are  not  perfectly  circular  in 
Figure  5-14  which  indicates  a  degree  of  probe  asymmetry. 
The  relatively  flat  section  of  the  constant  contours  of 
C0  in  Figure  5-15,  for  values  of  v/v  near  zero,  show 
relative  independence  of  p_  on  the  flow  angle  a  for 

9 

large  values  of  the  flow  angle  8  .  The  accuracy  of  the 
curve  fit  procedure  for  co  is  shown  by  the  relative 
position  of  the  labeled  data  points  with  respect  to  the 
lines  of  constant  0o. 

Constant  contours  of  Cq  plotted  against  Ca  and  Cp  are 
shown  in  Figures  5-16  and  5-17  for  the  first  and  second 
calibration  regions,  respectively.  Equation  5-7  defines 
Cq  as  the  ratio  of  the  dynamic  pressure  measured  by  the 
probe,  p  -  1/4P-  ,  to  the  actual  dynamic  pressure, 
PT  -  Ps  Figure  5-16  shows  that  Cq  has  a  maximum  value 
for  CQ-  0.25  &  Cp-  -0.35  and  that  Cq  has  values  greater 

that  0.475  for  values  of  Cp  less  than  -0.1.  This 
indicates  a  degree  of  physical  asymmetry  of  the  probe. 
Figure  5-17  shows  some  degree  of  physical  asymmetry  of 
the  probe  and  the  decreasing  trend  of  Cq  with  increasing 
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Cp  The  accuracy  of  the  curve  fit  for  C  is  shown  by  the 
relative  position  of  the  labeled  data  points  with  respect 
to  the  lines  of  constant  Cq . 

C.  Accuracy  -  probe  #1 

Table  5-1  shows  the  average  accuracy  achieved  over  all 
the  data  points  used  in  calibration. 


TABLE  5-1 

PARAMETER  ACCURACY 

PROBE  #  1 

Variable 

1.960 

C0 

0.036 

C<l 

0.017 

U/V 

0.49° 

V/V 

0.56° 

S  ince 

the  probe  was  used 

by 

approximately 

nulling  Ca 

before 

making 

measurements , 

a 

more  realistic 

measure  of 

the  accuracy 

is  obtained 

by 

comparison  of 

the  actual 

values 

and  the 

predicted  values 

of  all  parameters  for  the 

five 

central 

values  of  the 

flow  angle  a  (i.e 

.a  =  +3.8, 

«=  +1.9  and  .<x=  0).  Table  5-2  shows  the  average  accuracy 
achieved  over  this  set  of  data  points  for  the  four 
parameters  and  five  variables  of  interest. 
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TABLE  5-2 

CENTRAL  ACCURACY 

PROBE  #  1 

Variable 

1 .960 

C0 

0  .035 

Cq 

0.016 

u/V 

0.55° 

v/V 

0.51° 

PT 

.00083  psig 

ps 

.00194  psig 

V 

1.53  ft/sec 

6.2% 

.00092  psig 

With  95%  certainty,  the  above  table  gives  the  accuracy  in 
determining  the  listed  variable  from  the  five  pressures 
measured  by  the  probe.  As  indicated  above,  the  angles  are 
predicted  to  about  the  nearest  half  a  degree.  The  total 
pressure  is  indicated  to  3.5%  of  Pg  -  1/4  Pr  or  1.5%  of 
the  dynamic  pressure  and  the  static  pressure  is  indicated 
to  3.5%  of  the  dynamic  pressure.  The  velocity  is 
indicated  to  1.9%  of  the  actual  velocity.  The  axial 
momentum  flux,  Mz  ,  is  indicated  to  6.2%  of  the  actual 


axial 


momentum  flux  and  the  tangential  or  angular 
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momentum  flux,  Me  ,  is  indicated  to  1.7%  of  the  dynamic 
pressure. 


D.  Plow  Angles  -  probe  #2 

A  comparison  between  the  actual  values  of  u/v  and  v/v  and 
the  values  predicted  using  the  polynominals  of  Equation 
5-8  are  given  in  Figure  5-33.  The  error  in  determining 
the  flow  angle  is  indicated  by  the  difference  between 
each  set  of  predicted  and  actual  values  of  u/v  and  v/v. 
The  fit  is  most  accurate  near  the  origin  of  the 
tangential  velocity  ratio,  v/v  (flow  angle  a*  0).  The 
error  is  greatest  in  the  second  calibration  region 
(values  of  the  tangential  velocity  ratio,  v/v,  less  than 
-0.3)  than  in  the  first  calibration  region. 

Figures  5-27  and  5-28  are  plots  of  u/v  versus  v/v  showing 
lines  of  constant  Ca  and  Cp  for  the  first  and  second 
calibration  regions,  respectively,  it  is  apparent  from 
the  almost  vertical  Ca  lines  and  almost  horizontal  Cp 
lines  that  C0  is  almost  independent  of  u/v  and  Cp  almost 
independent  of  v/v  in  the  first  calibration  region.  The 
location  of  the  Cp  equal  zero  line  in  both  Figures  5-27 
and  5-28  gives  the  value  of  the  flow  angle  8  at  which  the 


pressure  coefficient  Cp  is  nulled.  The  location  of  the  C( 
equal  zero  line  in  Figure  5-27  gives  the  value  of  the 

flow  angle  a  at  which  the  pressure  coefficient  CQ  is 

nulled. 

E.  Total  and  Dynamic  pressures  -  probe  #2 

Constant  contours  of  C0  plotted  against  Cg  and  Cp  are 
shown  in  Figures  5-29  and  5-30  for  the  first  and  second 
calibration  regions,  respectively.  Figure  5-29  shows 
that,  for  small  angles,  p.  gives  an  accurate'  reading  of 
the  actual  total  pressure.  The  larger  the  angle,  the 

larger  the  difference  between  PT  and  pg  .  The  constant 
contours  of  C0  are  not  perfectly  circular  in  Figure  5-29 
which  indicates  a  degree  of  probe  asymmetry.  The  relative 
flat  section  of  the  constant  contours  of  C0  in  Figure 
5-30  for  values  of  v/v  near  -0.25  show  relative 

independence  of  p  on  the  flow  angle  6  for  large  values 
of  the  flow  angle  a  .  The  accuracy  of  the  curve  fit 
procedure  for  C0  is  shown  by  the  relative  position  of  the 
labeled  data  points  with  respect  to  the  lines  of  constant 

<v 

Constant  contours  of  Cq  plotted  against  Ca  and  Cp  are 


shown  in  Figures  5-31  and  5-32  for  the  first  and  second 
calibration  regions,  respectively.  Figure  5-31  shows  that 
Cq  has  a  maximum  value  for  Ca  =  -0.07  &  Cp  =  -0.08  and 
that  Cq  has  a  saddle  at  about  CQ  =  -0.3  &  Cp  -  0.24.  The 
elliptical  shape  of  the  contours  and  the  saddle  indicate 
a  degree  of  physical  asymmetry  of  the  probe.  Figure  5-32 
shows  contours  that  are  extensions  of  ellipses  for  values 
of  Cq  greater  than  0.5  and  asymmetric  contours  for  values 
of  Cq  less  than  0.5.  The  accuracy  of  the  curve  fit  for  Cq 
is  shown  by  the  relative  position  of  the  labeled  data 
points  with  respect  to  the  lines  of  constant  Cq  . 

F.  Accuracy  -  probe  #2 

Table  5-3  shows  the  average  accuracy  achieved  over  the 
148  data  points  used  in  calibration  for  the  four 
parameters  and  five  variables  of  interest.  With  95% 
certainty,  Table  5-3  gives  the  accuracy  in  determining 
the  listed  variable  from  the  five  pressures  measured  by 
the  probe.  The  total  pressure  is  indicated  to  3.0%  of  the 
dynamic  pressure  and  the  static  pressure  is  indicated  to 
11.8%  of  the  dynamic  pressure.  The  velocity  is  indicated 
to  5.2%  of  the  actual  velocity.  The  axial  momentum  flux, 
M2 ,  is  indicated  to  9.0%  of  the  actual  momentum  flux  and 


i5fc 
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the  tangential  or  angular  momentum  flux,  M0 ,  is  indicated 
to  4.5%  of  the  dynamic  pressure.  The  flow  angles  are 
predicted  to  about  the  nearest  1.5  degrees. 


TABLE 

5-3 

CENTRAL  ACCURACY 

PROBE 

#  2 

variable 

1.960 

C0 

0.054 

0.062 

U/V 

1.5 

V/V 

1.3 

PT 

.00136  psig 

ps 

.00539  psig 

V 

3.90  ft  /sec 

M 

9.0* 

Me 

.00206  psig 

The  accuracy  of  probe  #2,  shown  above,  is  significantly 
different  than  that  of  probe  #1,  shown  in  Table  5-2,  due 
to  the  very  large  range  of  flow  angles  that  probe  #2  will 
measure,  probe  #l's  accuracy  is  for  the  approximate 
nulling  of  the  flow  angle  a  and  probe  #2's  accuracy 
includes  values  of  the  flow  angle  a  as  large  as  -50 
degrees. 
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Figure  5 

-6  pressure  Coefficients  for  Probe  #1 

ALPHA  IN  DEGREES 
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O  ALPHA  VERSUS  CB  17  VALUES 
a  ALPHA  VERSUS  CO  37  VALUES 
+  ALPHA  VERSUS  CQ  17  VALUES 

Figure  5-7  pressure  coefficients  for  probe  #1,  B» 
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Figure  5-8  Pressure  coefficients  for  Probe  #1 
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PRESSURE  COEFFICIENTS  VERSUS  ALPHA  FOR  BETA  =  25  DEGREES 
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0  ALPHA  VERSUS  CB  17  VALUES 
a  ALPHA  VERSUS  CO  17  VALUES 
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Figure  5-9  pressure  coefficients  for  probe  #1,  B«  25 
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5-22  pressure  Coefficients  for  probe  #2,  8 
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Figure  5-23  pressure  Coefficients  for  probe  #2,  6 
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Figure  5-24  pressure  Coefficients  for  probe  #2,  6 
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PRESSURE  COEFFICIENTS  VERSUS  ALPHA  FOR  BETA  =  20  DEGREES 
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Figure  5-25  pressure  Coefficients  for  probe  #2,S»  20 
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Figure  5-33  Actual  and  predicted  values  of  u/V 
and  v/v  for  probe  #2 


CHAPTER  6 


DATA  REDUCTION 

I.  INTRODUCTION 

The  data  set  for  each  test  run  was  reduced  on  the 
University  of  Washington's  CDC  computer  using  the 
following  four  step  procedure: 

1.  The  recorded  raw  data  set  was  loaded  from 
magnetic  tape  onto  a  disc  file  from  a  remote  CRT  terminal 
via  campus  phone  lines. 

2.  The  raw  data  set  was  initially  reduced  by  time 
averaging  certain  data  points  and  removing  the  time 
dependent  zero  drift  of  each  low  pressure  transducer 
measurement. 

3.  The  flow  properties  for  the  test  run  were 
calculated  and  recorded  on  disc  file  for  later  plotting. 

4.  The  reduced  data  sets  with  the  same  mass  flow 
ratio  were  plotted  for  analysis. 

The  following  sections  describe  each  step  of  data 
reduction  as  shown  in  Figures  6-la  and  6-lb. 


II.  LOADING  RECORDED  RAW  DATA 


A  set  of  raw  data  for  each  test  run  was  recorded  on 
magnetic  tape  using  the  Fluke  Datalogger  and  the  Columbia 
300  Data  Tape  Recorder  as  shown  at  the  top  of  Figure 
6-la.  The  Columbia  300  Data  Tape  Recorder  was  then 
connected  in  series  between  a  Televideo  920C  CRT  Terminal 
and  a  1200  baud  MODEM  for  entry  of  each  set  of  raw  data 
into  the  computer.  Communications  was  then  established 
between  the  remote  terminal  and  the  CDC  computer  using 
the  campus  phone  system.  The  raw  data  set  recorded  on  the 
magnetic  tape  and  manual  entries  from  the  CCT  terminal 
were  read  by  the  computer's  ICE  Text  Editor  and  then 
saved  on  disc  as  a  RAW  DATA  FILE  (shown  in  Figure  6-la). 
The  ICE  Text  Editor  was  selected  to  read  the  data  because 
this  text  editor  used  a  carriage  return  to  indicate  the 
end  of  a  line  of  data  (the  same  format  as  the  data 
recorded  on  magnetic  tape  using  the  Fluke  Datalogger) . 
Also  this  text  editor  sent  a  sequential  line  number 
prompt  for  each  line  of  data  read  that  was  monitored  for 
the  progress  and  status  of  loading.  The  format  of  the 
initial  part  of  a  RAW  DATA  FILE  is  shown  in  Table  D-2 
where  the  first  two  lines  are  manual  entries  from  the  CRT 
terminal.  The  number  of  subruns,  the  subrun  number  and 
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the  number  of  scans  for  a  subrun  were  manually  entered 
from  the  CRT  terminal.  The  number  of  each  test  run  was 
incorporated  into  that  data  file's  filename  for  later 
identification.  After  loading  of  the  raw  data  file  for  a 
test  run,  this  file  was  available  for  the  first  step  of 
data  reduction. 

III.  INITIAL  DATA  REDUCTION  -  PROGRAM  "REDUCE" 

Each  raw  data  file  for  a  test  run  consisted  of  twenty  or 
more  subruns.  Each  subrun  corresponded  to  a  location  of 
the  five-hole  probe  and  consists  of  49  datalogger  scans 
of  the  18  measured  quantities  (the  time  and  17  voltages) . 
Scans  #2  through  #49  corresponded  to  positions  1  through 
48  of  the  scanivalve  and  scan  #1  corresponded  to  the 
initial  scan  of  position  48  of  the  Scanivalve. 

The  five  lengths  of  tygon  tubing  between  the  five-hole 
probe  and  its  six  pressure  transducers  were  disconnected 
during  scans  #1  and  #49,  allowing  zero  pressure 
measurements  of  these  six  pressure  transducers  and  the 
two  scanivalve  pressure  transducers  during  scans  #1  and 
#49.  Considerable  drift  in  the  output  from  each  low 
pressure  transducer  was  experienced  and  this  drift  is 
shown  in  Figure  6-2.  The  magnitude  of  the  time  variation 
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in  the  output  voltage  from  each  low  pressure  transducer 
was  found  to  be  independent  of  the  input  pressure's 
magnitude.  The  period  of  the  time  variation  was  found  to 
be  much  larger  than  the  length  of  time  required  for  a 
subrun,  tf  -t,  •  program  REDUCE  calculated  the  output 

voltage  for  a  zero  input  pressure,  Va(t),  for  each  low 
pressure  transducer  from  the  data  gathered  on  scan  #1 
(initial)  and  scan  #49  (final)  and  assuming  that  v0  (t) 
varied  linearly  between  tf  and  t(  using  the  following 
relationship: 

Va  (tf  >  ”  Vo  (ti  > 

VG  (t)  -  vo  (t,  )  +  - - x  (t  -t,  ) 

fci  “  fc. 

The  time  variation  of  the  output  voltage  from  the  low 
pressure  transducers  was  removed  by  calculating  the 
corrected  transducer  voltage,  V  ,  using  the  measured 
output  voltage,  V(t),  and  subtracting  the  output  voltage 
for  zero  input  pressure,  vo  (t),  at  the  same  time,  t, 
calculated  from  the  above  equation. 

program  REDUCE  averaged  the  following  inputs  over  fourty 
measurements  (scan  #6  through  scan  #45)  for  each  subrun 
of  a  test  run: 

1.  The  corrected  transducer  voltage  for  the  six  low 
pressure  transducers  connected  to  the  five-hole  probe. 

2.  The  voltages  from  the  tank  and  line  pressure 


transducers 


3.  The  voltages  from  the  two  flow  orifice  plate 
differential  pressure  transducers. 

4.  The  voltage  from  the  potentiometer  corresponding 
to  the  probe  angle  ap  . 

5.  The  temperatures  of  the  inside  and  outside 
plenums  and  the  LVDT  output  voltage  corresponding  to  the 
probe's  radial  position. 

The  output  from  REDUCE  for  each  subrun  consisted  of  the 
subrun  number,  the  hour  and  minute  at  the  start  of  the 
subrun,  the  corrected  transducer  voltage  for  scan  |l 
through  scan  #49  of  both  Scanivalve  tranducers  and  the 
averaged  measurements  listed  above.  Typical  output  from 
REDUCE  is  shown  in  Table  D-3  (FINE  DATA  FILE)  for  Run  #4. 
The  first  four  lines  of  the  FINE  DATA  FILE  in  Table  D-3 
are  test  conditions  which  were  manually  loaded  onto  the 
data  file  from  the  CRT  terminal  and  the  remainder  of  the 
FINE  DATA  FILE  is  output  from  program  REDUCE  for  two  of 
the  twenty-nine  subruns  in  this  test  run. 

This  initial  data  reduction  is  shown  in  the  bottom 
portion  of  the  flow  chart  in  Figure  6-la.  program  REDUCE 
and  a  description  of  input  and  output  variables  for  this 
program  are  contained  in  Appendix  d. 


IV.  CALCULATION  OF  FLOW  PROPERTIES  -  PROGRAM  "DATAR4 " 

Program  DATAR4  used  the  FINE  DATA  FILE  for  each  test  run 
as  the  input  data  file  as  is  shown  in  Figure  6-lb. 
program  DATAR4  also  required  the  following  data  files: 

1.  probe  calibration  Coefficients  of  Regions  1  and  2 

2.  Layout  of  Pressures  on  Scanivalve 

program  DATAR4  and  a  description  of  input  and  output 
variables  for  DATAR4  are  contained  in  Appendix  E.  Program 
DATAR4  has  seven  subroutines  to  calculate  specific 
properties  and  one  subroutine  to  plot  out  the  results  on 
the  printer.  The  main  program  portion  of  DATAR4  performs 
all  input  of  data  from  disc  files  and  output  of  results 
to  the  printer  (except  the  printer  plots)  and  to  the 
output  disc  file,  REDUCED  DATA  FILE.  The  main  program 
portion  of  DATAR4  also  averaged  the  wall  static  pressures 
over  all  subruns  within  the  test  run  and  calculated  the 
axial  momentum,  the  angular  momentum,  the  swirl  number, 
the  average  axial  velocity  and  the  mass  average  total 
pressure  from  the  five-hole  probe  measurements  at  the 
axial  location  of  the  test  run.  The  seven  subroutines 
called  by  the  main  program  portion  of  DATAR4  are 
explained  in  the  following  paragraphs  in  the  order  they 
are  contained  in  Appendix  E. 


Subroutine  PVPLOT  plots  out  the  results  on  the  printer, 
using  the  Numerical  plotting  system  (NPS)  software  of  the 
University  of  Washington's  CDC  computer.  This  subroutine 
produced  the  following  eleven  plots  for  the  axial 
location  of  the  test  run: 

1.  wall  static  pressure  versus  axial  location  on  the 
four-inch  diameter  centerbody. 

2.  Wall  static  pressure  versus  axial  location  on  the 
eight-inch  diameter  outer  wall. 

3.  wall  static  pressure  versus  axial  location  on  the 
inside  of  the  six-inch  diameter  tube  separating  the  two 
air  streams. 

4.  The  following  probe  measured  quantities  versus 
the  radial  position  of  the  probe  (r): 

a.  Axial  velocity  (w) 

b.  Radial  velocity  (u) 

c.  Tangential  velocity  (v) 

d.  Total  velocity  (V) 

e.  static  pressure  (p$  ) 

f.  Total  pressure  (PT ) 

g.  Axial  momentum  flux  (Mz  or  Mz  ) 

h.  Tangential  momentum  flux  (M#) 


Subroutine  LOCATN  determined  the  location  and  the 
pressure  o£  each  static  pressure  tap,  using  the  coded 
Scanivalve  location  for  each  of  the  74  static  pressure 
taps  (contained  in  FILE  1,  LAYOUT)  and  the  starting 
location  of  each  section  of  the  eight-inch  diameter  outer 
wall  (contained  in  FILE  2,  FINE  DATA  FILE  for  the  test 
run).  The  file  LAYOUT,  containing  the  coded  scanivalve 
location  for  each  of  the  74  static  pressure  taps,  is 
shown  in  Table  E-2.  The  starting  location  of  each  section 
of  the  eight-inch  diameter  outer  wall  for  test  run  #4  is 
shown  in  the  fourth  line  of  the  FINE  DATA  FILE  in  Table 
D-3. 

Subroutine  PRESS  calculated  the  pressure  measured  by  the 
twelve  pressure  transducers,  using  the  calibration 
coefficient  of  the  respective  transducer,  the  measured 
output  voltage  of  the  pressure  transducer  and  the  output 
voltage  of  the  pressure  transducer  for  zero  input 
pressure.  The  output  voltage  of  the  eight  low-pressure 
transducers  had  already  been  corrected  in  the  initial 
data  reduction  program  REDUCE.  The  output  voltages  of  the 
other  four  pressure  transducers  for  zero  input  pressure 
are  contained  in  the  third  line  of  the  FINE  DATA  FILE  as 
shown  in  Table  d-3  for  test  run  #4. 


Subroutine  DENSTY  calculated  the  density  o£  the  air  in 
both  plenums  from  the  plenum's  total  pressure  and 
temperature.  This  subroutine  also  calculated  the  mass 
average  density,  the  axial  velocity  and  the  dynamic 
pressure  based  on  this  axial  velocity  for  the  test 
section.  This  subroutine  used  the  mass  flow  rates  of  the 
two  air  streams  determined  by  the  subroutine  FLOW. 

Subroutine  POSITN  calculated  the  five-hole  probe's  radial 
position  for  each  subrun  based  on  the  voltage  from  the 
Linear  variable  Differential  Transformer  (LVDT) ,  the 
calibration  coefficients  of  the  LVDT  and  the  reference 
voltage  and  position  of  the  probe.  The  reference  radial 
position  and  its  corresponding  voltage  are  the  last  two 
entries  contained  in  the  first  line  of  the  FINE  DATA  FILE 
as  shown  in  Table  D-3  for  test  run  #4. 

Subroutine  FLOW  calculated  the  mass  flow  rates  through 
both  plenums,  using  two  different  methods.  The  first 
method  of  calculating  the  mass  flow  rate  is  based  on 
choked  flow  and  knowing  the  area  of  the  choked  orifice, 
the  total  pressure  upstream  of  the  choked  orifice  and  an 
estimate  of  the  total  temperature  upstream  of  the  choked 
orifice.  The  second  method  of  calculating  the  mass  flow 
rate  is  based  on  a  calibrated  flow  orifice  plate  and 


knowing  the  differential  pressure  across  the  flow  orifice 
plate,  the  total  pressure  upstream  of  the  flow  orifice 
plate  and  an  estimate  of  the  total  temperature  upstream 
of  the  flow  orifice  plate. 

Subroutine  PROBE  calculated  the  flow  properties  at  the 
five-hole  probe's  radial  position,  using  the  probe's 
calibration  coefficients  for  regions  #1  and  #2,  the 
pressures  at  each  of  the  probe's  five  ports,  the  voltage 
from  the  probe's  angular  position  measuring  system  and 
the  method  described  in  the  Chapter  5.  The  following  flow 
quantities  were  determined  for  the  point  in  the  flow 
measured  by  the  probe: 

1.  Axial  velocity  to  total  velocity  ratio  (w/v) 

2.  Radial  velocity  to  total  velocity  ratio  (u/v) 

3.  Tangential  velocity  to  total  velocity  ratio  (v/v) 

4.  Total  velocity  (v) 

5.  Flow  angle  alpha  (a) 

6.  Flow  angle  beta  (3) 

7.  Static  pressure  (pg  ) 

8.  Total  pressure  (pt  ) 

9.  Axial  momentum  flux  (Mz  or  Mz  ) 

10.  Tangential  momentum  flux  (Me) 


198 


Subroutine  INTERGR  used  the  trapezoid  rule  and  data  at 
each  radial  position  to  integrate  flow  properties  between 
the  inner  radius  and  the  outer  radius  of  the  flow 
channel.  This  subroutine  was  used  by  the  main  portion  of 
0ATAR4  to  calculate  the  axial  momentum,  the  angular 
momentum,  the  swirl  number,  the  average  axial  velocity 
and  the  mass  average  total  pressure. 

V.  OUTPUT  FROM  PROGRAM  "DATAR4 " 

Typical  printer  output  from  program  DATAR4  is  shown  in 
Table  e-3  in  Appendix  E  for  test  run  #061.  The  printer 
output  includes  listings  of  the  following  quantities: 

1.  input  data  for  the  test  run 

2.  Wall  static  pressure  versus  axial  position 

3.  Output  data  from  probe  measurements 

4.  printer  plots 

program  DATAR4  also  wrote  the  results  of  each  test  run 
onto  a  separate  disc  file  (Reduced  Data  File  XXX)  as 
shown  in  Figure  6-lb.  The  Reduced  Data  Files  were  used  as 
input  data  for  the  plotting  program  PLDATA  that  is 
discussed  in  the  next  section. 


VI.  PLOTTING  PROGRAM  "PLDATA" 


program  PDATA  plots  the  reduced  data  files  as  shown  in 
Figure  6-lb.  This  program  was  developed  for  plotting 
results  using  the  TEKTRONIX  PLOT-10/Terminal  Control 
System  contained  on  the  PDP-11/23  computer  system  in  the 
Aeronautics  Laboratory  at  the  united  states  Air  Force 
Academy,  program  PLDATA  and  a  description  of  input  and 
output  variables  for  this  program  are  contained  in 
Appendix  F. 

The  main  portion  of  this  program  performed  all  input  of 
data  from  disc  files  and  the  interactions  with  the  user. 
The  subroutine  PVPLOT  accomplished  all  the  plotting  of 
the  results  by  directly  calling  on  PLOT-10  subroutines 
and  by  indirectly  calling  on  PLOT-10  subroutines  in  the 
subroutines  SYMBLS  and  GRID.  Subroutine  SCALE  performed 
the  necessary  scaling  of  the  data  before  plotting. 
Subroutine  SYMBLS  plots  one  of  eleven,  possible  symbols  at 
each  data  point.  Subroutine  GRID  draws  a  two  dimensional 
grid  and  numbers  both  axes. 
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Figure  6-la  Data  Reduction  Flow  Chart,  Part  I 
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CHAPTER  7 


RESULTS  &  CONCLUSIONS 

I.  INTRODUCTION 

This  chapter  discusses  the  experimental  results  of  the 
four  tests  conducted  to  measure  the  turbulent  mixing  of 
coaxial  air  streams  with  swirl  present.  The  data  for  the 
tests  is  presented  in  Appendix  A  and  the  experimental 
measurements,  test  apparatus,  five-hole  probes,  and  data 
reduction  are  discussed  in  Chapters  3,  4,  5,  and  6, 

respectively. 

II.  TEST  #1 

A  limited  amount  of  data  was  taken  for  the  mixing  of  a 
swirling  stream  from  an  annulus  in  a  flat  wall  into  a 
stagnant  environment  (the  test  configuration  is  shown  in 
Figure  3-1)  .  The  reduced  data  for  this  test  is  contained 
in  Figures  A-l  through  A-8  in  Appendix  A.  The  inside 
boundary  of  the  swirling  stream  grows  toward  the 
centerline  and  the  outside  boundary  of  this  stream  grows 
away  from  the  centerline  into  the  stationary  surrounding 
air.  The  growth  of  thre  outside  boundary  of  the  swirling 
stream  is  expected  bicause  this  boundary  is  unstable  due 
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to  the  decrease  in  fluid  angular  momentum  in  the  radial 
direction.  The  inward  growth  of  the  swirling  stream’s 
inside  boundary  can  be  attributed  to  the  zone  of 
recirculation  (sketched  below)  that  is  present  about  the 
centerline  near  the  entrance. 


Outar _ BoundarjL 


III.  TEST  #2 

Data  was  taken  of  the  mixing  of  two  coaxial  air  streams 
(inner  stream  with  swirl)  from  concentric  annuli  in  a 
flat  wall  into  stagnant  surroundings  (see  Figure  3-2  for 
test  configuration) .  The  reduced  data  for  this  test  is 
contained  in  Figures  A-9  through  A-17  in  Appendix  A.  The 
outside  boundary  of  the  inner  swirling  stream  does  not 
grow  as  quickly  when  surrounded  by  a  coaxial  air  stream 
as  was  the  case  for  this  boundary  in  stagnant 
surroundings  (Test  #1).  The  reduced  growth  of  the  outside 
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boundary  of  the  swirling  stream  can  be  attributred  to  the 
presence  of  the  outer  stream  and  the  reduction  in  the 
intensity  of  the  instability  of  this  boundary.  The  outer 
stream  has  a  much  higher  velocity,  axial  momentum,  and 
static  pressure  than  was  the  case  for  Test  Set  #1.  The 
inside  boundary  of  the  swirling  stream  grows  toward  the 
centerline  due  to  the  zone  of  recirculation  that  is 
present  about  the  centerline  near  the  entrance. 

IV.  TEST  #3 

The  test  section  for  the  third  set  of  experiments  is 
shown  in  Figure  3-3.  Two  streams  of  air  (inner  stream 
with  swirl)  from  concentric  annuli  in  a  flat  wall  flow 
into  an  open  region.  This  set  of  experiments  differed 
from  the  second  set  of  experiments  by  the  addition  of  an 
extension  of  the  inner  wall  of  the  inner  annulus 
(centerbody)  into  the  open  region.  The  reduced  data  for 
this  test  is  contained  in  Figures  A-18  through  A-35  in 
Appendix  A.  The  presence  of  the  centerbody  prevents  the 
inward  growth  of  the  swirling  stream's  inside  boundary 
and  the  growth  rate  of  the  swirling  stream's  outside 
boundary  is  increased  over  its  growth  rate  in  Test  #2. 

The  scatter  of  data  for  measurements  at  the  two  far 
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downstream  axial  locations  is  attributed  to  the 
intermittency  of  the  turbulent  flow  at  these  locations. 

V.  COMPARISON  OF  TEST  #1,  #2,  AND  #3 

The  growth  of  the  outside  boundary  of  the  swirling  stream 
is  presented  for  the  first  three  tests  in  Figure  7-1  as  a 
plot  of  the  tangential  velocity's  half-radius  versus 
axial  location.  Both  the  absence  of  a  centerbody  and  the 
presence  of  an  outer  coflowing  stream  retard  the  growth 
of  the  swirl  velocity. 

VI.  TEST  #4 

The  final  and  most  extensive  set  of  experiments  was 
performed  with  the  test  section  in  the  configuration 
shown  in  Figure  3-4.  In  this  configuration,  two  streams 
of  air  (inner  stream  with  swirl)  from  concentric  annuli 
flow  into  a  constant  area  annular  duct.  The  reduced  data 
for  this  set  of  experiments  is  contained  in  Figures  A-36 
through  A-90  in  Appendix  a.  Measurements  were  taken  for 
five  different  mass  flow  ratios,  o's,  or  runs.  Table  7-1 
lists  these  five  mass  flow  ratios,  the  three  character 
abbreviations  used  to  differentiate  the  data  runs  of  each 
mass  flow  rate,  the  Figures  in  Appendix  A  containing  the 


reduced  data  for  the  data  runs,  and  the  swirl  number, 
(defined  by  Equation  G-16  in  Appendix  G),  at  the  first 
axial  location,  z  =*  0.25  inches. 

TABLE  7-1 

MASS  FLOW  RATIOS  FOR  TEST  #4 


Mass  Flow  Ratio 
(  a  =  m0  /ift,  ) 

Abbreviation 
of  Data  Runs 

Figures 

Swirl 

So 

0.00 

XXI 

A- 3  6  -  A- 4  6 

.1762 

0.47 

XX2 

A- 4 7  -  A-57 

.1579 

1.00 

XX3 

A-58  -  A-68 

.1096 

2.13 

XX4 

A- 6  9  -  A- 7  9 

.0379 

3.91 

XX5 

A- 80  -  A-90 

.0125 

In  runs  xxi  (  a  =«  0.00),  a  region  of  recirculation  is 
established  at  the  entrance  of  the  test  section  between 
the  inner  swirling  stream  and  the  outside  wall.  At  about 
four  inches  downstream  of  the  test  section  entrance,  the 
inner  swirling  stream  separates  from  the  centerbody  due 
to  the  large  adverse  pressure  gradient  (see  Figure  7-2) 
and  a  region  of  flow  reversal  is  produced.  Further 
downstream  (at  about  eight  inches  from  entrance) ,  the 
flow  re-attaches  to  the  centerbody.  The  streamline 
pattern  for  runs  XXI  is  sketched  in  Figure  7-3. 

In  runs  XX2  ( a  »  0.47),  the  flow  from  the  outer  annulus 
prevents  the  establishment  of  an  outer  region  of 


recirculation  and  reduces  the  magnitude  of  the  adverse 
pressure  gradient  on  the  centerbody  (see  Figure  7-4) 
preventing  separation.  The  flow  from  the  outer  annulus 
also  reduces  the  rate  of  outward  growth  of  the  tangential 
velocity. 

In  runs  XX3  (  1.00),  the  adverse  pressure  gradient  on 
the  centerbody  is  further  reduced  by  the  flow  from  the 
outer  annulus  (see  Figure  7-5).  The  initial  difference  in 
axial  momentum  flux  between  the  inner  and  outer  streams 
is  lower  than  for  runs  XX2  and  the  rate  of  outward  growth 
of  the  tangential  velocity  is  further  reduced. 

in  runs  XX4  (  a  *  2.13),  the  outside  wall  of  the  test 
section  and  the  higher  velocity  of  the  outside  stream 
cause  the  boundary  between  the  two  streams  to  move  inward 
before  any  significant  outward  growth  of  tangential 
velocity.  The  inward  qrowth  of  axial  momentum  flux 
reduces  the  outward  transfer  of  tangential  momentum  flux 
and  its  associated  tangential  velocity.  The  magnitude  of 
the  adverse  pressure  gradient  present  on  the  initial 
section  of  centerbody  (see  Figure  7-6)  is  reduced  by  the 
inward  growth  of  the  outer  stream. 


In  runs  XX5  (  a=  3.91),  both  walls  of  the  test  section 
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have  adverse  pressure  gradients  (see  Figure  7-7) .  The 
adverse  pressure  gradient  on  the  outside  wall  is  due  to 
the  inward  spread  of  the  outer  stream.  The  inward 
transfer  of  axial  momentum  flux  further  reduces  the 
outward  transfer  of  tangential  momentum  flux  as  compared 
to  runs  XX4. 

The  half-radius  growth  of  the  swirl  (tangential)  velocity 
versus  axial  location  is  plotted  in  Figure  7-8  for  runs 
XXI,  XX2 ,  XX3,  XX4,  and  XX5  of  test  #4.  The  increase  in 
outward  growth  with  increased  swirl  agrees  with  the 
results  for  a  swirling  jet  (Reference  1) .  The  rate  of 
outward  growth  of  the  swirl  velocity  is  reduced  by  a 
flowing  outer  stream  and  the  higher  the  mass  flow  ratio 
(with  its  outer  stream  axial  velocity  and  axial  momentum 


flux)  , 

the  slower  the 

rate 

of 

this 

growth.  When 

the 

initial 

axial  momentum 

flux 

of 

the 

outer  stream 

is 

greater  that  that  of  the  inner  stream  (runs  XX4  and  XX5) , 
the  rate  of  inward  transfer  of  axial  momentum  more  than 
offsets  the  outward  transfer  of  tangential  momentum  in 
the  initial  portion  of  the  test  section  and  the 
half-radius  of  the  swirl  velocity  is  reduced. 

The  variation  of  the  swirl  number  (S) ,  with  axial 
distance  for  runs  XXI,  XX2 ,  XX3 ,  XX4 ,  and  XX5  of  test 
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#4, is  plotted  in  Figure  7-9.  The  swirl  number  has  been 
normalized  by  its  value  at  the  £irst  axial  location  (s  ) . 
This  initial  value  of  the  swirl  number  is  listed  in  Table 
7-1.  The  results  indicate  that  the  decay  of  the  swirl 
number  with  axial  distance  for  runs  XXI,  XX2,  and  XX3  is 
dominated  by  the  decay  of  the  tangential  momentum  with 
axial  distance.  This  agrees  with  the  results  of  swirl 
dominated  flows  (Reference  2).  The  results  for  runs  XX4 
seem  to  indicate  that  this  region  of  flow  is  dominated  by 
the  decay  of  the  axial  momentum  with  axial  distance.  The 


author 

is 

unable 

to 

explain 

the  variation  of 

the  swirl 

number 

for 

runs 

XX5 

except 

to  note  that  the 

five-hole 

probe  was  being  used  in  the  outer  edge  of  its  calibration 
and  small  errors  in  predicted  flow  direction  had  a  much 
larger  effect  on  the  tangential  momentum  flux  than  on  the 
axial  momentum  flux. 

The  variation  of  the  static  pressure  along  the  outside 
wall  of  the  inner  annulus  upstream  of  the  test  section  is 
presented  in  Figures  7-10  and  7-11  for  runs  XXI,  XX2 , 
XX3 ,  XX4 ,  and  XX5  of  test  #4.  The  static  pressures  in 
Figure  7-10  have  been  nondimensionalized  by  the  dynamic 
pressure  associated  with  the  mean  axial  velocity  in  the 
test  section  and  those  in  Figure  7-11  have  been 
nondimensionalized  by  the  dynamic  pressure  associated 
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with  the  mean  axial  velocity  in  the  inner  annulus.  The 
influence  of  the  outer  stream  on  the  pressure  along  the 
outside  wall  of  the  inner  annulus  can  be  seen  for  runs 
XX3,  XX4,  and  XX5  by  comparing  Figures  7-10  and  7-11.  For 
runs  XXI  and  XX2,  the  static  pressure  on  the  outside  wall 
and  at  the  end  of  this  wall  are  only  a  function  of  the 
dynamic  pressure  of  the  inner  stream  and  the  axial 
location  (see  Figure  7-11) .  The  upstream  influence  for 
runs  XX3 ,  XX4 ,  and  XX5  can  be  seen  in  Figure  7-11  (i.e. 
the  higher  the  mass  flow  ratio,  a,  the  further  upstream 
the  influence) . 

The  variation  of  the  static  pressure  along  the  outside 
wall  of  the  inner  annulus  upstream  of  the  test  section  is 
presented  in  Figure  7-12  for  test  sets  #1,  #2,  #3,  and  |4 
(runs  XXI  and  XX2  only).  The  decrease  in  static  pressure 
for  test  sets  #1,  #3,  and  #4  indicates  an  increase  in  the 
fluid  velocity  near  the  outside  wall  at  the  exit. 
However,  the  increase  in  static  pressure  for  test  set  #2 
indicates  a  decrease  in  fluid  velocity  near  the  outer 
wall  at  the  exit  for  this  test  set. 

Static  pressure  measurements  in  test  #4  along  the 
centerbody  and  the  outside  wall  were  used  to  calculate 
the  coefficient  of  friction  (cf  )  for  this  flow.  The 


coefficient  of  friction  for  both  the  inside  wall  and  the 
outside  wall  are  presented  in  Table  7-2  for  each  test. 
Calculation  of  the  coefficient  of  friction  was  based  on 
the  dynamic  pressure  associated  with  the  mean  axial 
velocity  in  the  test  section. 

TABLE  7-2 

COEFFICIENT  OF  FRICTION ,  TEST  #4 


RUNS 

Centerbody 

Outside  Wall 

XXI 

0.0016 

0.0140 

XX2 

0.0031 

0.0108 

XX3 

0.0059 

0.0093 

XX4 

0.0048 

0.0048 

XX5 

0.0043 

0.0043 

The  higher  values  of  the  coefficient  of  friction  at  the 
outer  wall  reflect  the  existence  of  larger  tangential 
velocities  near  this  wall  and  the  higher  turbulence 
present  due  to  this  extremely  unstable  boundary  layer. 
Scott  and  Rask  (Reference  3)  explain  that  this  increase 
in  the  coefficient  of  friction  at  the  outer  wall  is  due 
to  unstable  flow  which  increases  the  turbulent  production 
and  promotes  turbulent  transfer. 

VII.  CONCLUSIONS 

On  the  basis  of  the  experimental  results,  the  following 
remarks  can  be  made  regarding  the  initial  mixing  region 
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(for  the  swirl  rates  tested): 

1.  For  annular  swirling  jet,  the  addition  of  a 
concentric  outside  annular  jet  decreases  the  growth  of 
the  outer  boundary  of  tangential  (angular)  momentum  in 
the  initial  region. 

2.  For  concentric  annular  streams  without  an  outside 
wall,  the  addition  of  a  centerbody  increases  the  radial 
transfer  of  tangential  (angular)  momentum. 

3.  For  concentric  annular  streams  with  a  centerbody, 
the  addition  of  an  outside  wall  decreases  the  radial 
transfer  of  tangential  (angular)  momentum. 

4.  For  concentric  annular  streams  with  a  centerbody 
and  an  outside  wall, 

a.  an  increase  in  the  axial  momentum  of  the  outer 
stream  reduces  the  radial  transfer  of  tangential 
(angular)  momentum. 

b.  an  increase  in  the  mass  flow  ratio  (outside  to 
inside  flow  rates)  increases  the  influence  of  the  outside 
stream  on  upstream  pressures  on  the  outer  wall  of  the 
inner  annulus. 
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APPENDIX  A 


EXPERIMENTAL  RESULTS 


This  appendix  presents  the  experimental  data  gathered  in 
four  sets  of  experiments  using  the  test  apparatus 
described  in  Chapter  4  and  the  five-hole  probes  described 
in  Chapter  5.  The  experimental  data  was  reduced  and 
plotted  in  Figures  A-l  through  A-90  using  the  programs 
described  in  Chapter  6.  The  experimental  data  has  been 
nondimensionalized  for  presentation  and  ease  of  analysis. 
The  quantities  used  to  nondimensional ize  are  shown  in 
Table  A-l.  Table  A-2  lists  the  test  section  configuration 
and  data  plot  Figures  for  each  test  number.  Table  A-3 
presents  the  values  used  to  nondimensionalize  the  data 
and  the  average  reference  values  for  each  test. 

plots  of  the  static  pressure  on  the  4  inch  tube,  the  6 
inch  tube,  and  the  8  inch  tube  for  tests  #3  and  #4  have 
individual  data  points  for  each  run  contained  in  the 
test.  However  the  plots  of  the  static  pressure  on  the  6 
inch  tube  for  tests  #1  and  #2  (Figures  A-l  &  A-9)  have 
data  points  that  are  the  average  of  the  runs  in  the  test. 
The  variations  of  velocity,  pressure,  momentum,  etc.  with 
radial  and  axial  position  for  tests  #1  and  #4  are  plotted 
in  a  format  having  individual  subplots  for  each  axial 


239 


position  measured.  The  variations  of  these  properties  for 
tests  #2  and  #3  are  plotted  on  a  combined  radial/axial 
position  scale  with  error  bounds  (dots)  plotted  for  each 
data  point.  The  horizontal  location  of  the  variable's 
zero  is  repositioned  to  correspond  to  the  axial  location 
of  the  test  run  as  shown  on  the  axial  position  scale  at 
the  top  of  these  plots  (see  Figure  A-21). 


TABLE  A-l 

Dimensionless  Quantities  for  Data  plots 


Dimensional 

Dimensionless 

Symbol 

Quantity 

Quantity 

w 

Axial  velocity 

w/U 

V 

Tangential  velocity 

v/u 

u 

Radial  Velocity 

u/U 

V 

Total  velocity 

V/U 

PT 

Total  pressure 

(!>-%)/ (1/2  PU!  ) 

ps 

Static  pressure 

(Ps-%)/(l/2  0  U2  ) 
(P-R)  /  (1/2  P  U2  ) 

P 

Wall  static  pressure 

M« 

Axial  Momentum  Flux 

Mj  /  (1/2  P  U2  ) 

Mi 

Axial  Momentum  Flux 

Mj/  (1/2  P  U2  ) 

Me 

where 

Tangential  Momentum  Flux 

M#/(l/2  p  U2  ) 

Symbol 

Description 

Test 

Sets 

U 

Mean  Axial  velocity  of  inner  stream 

1-3 

U 

Mean  Axial  velocity  in  Test  Section 

4 

Po 

Atmospheric  pressure 

1-3 

P0 

Static  pressure  on  outer  Wall  of 

Inner  Stream  at  z  ■  -2  inches 

4 

cn  «-3 


21+0 

TABLE  A- 2 

Test  Sets/Configuration/Figures 


Test 

Set 

Number 

Test 

Section 

Figure 

|  Configurat: 

ion 

Data 

Figures 

Outer 

Stream 

Center 

Body 

Outside 

Wall 

1 

3-1 

No 

No 

No 

A- 1  -  A- 8 

2 

3-2 

Yes 

No 

No 

A- 9  -  A- 17 

3 

3-3 

Yes 

Yes 

No 

A- 18  -  A- 3  5 

4 

3-4 

Yes 

Yes 

Yes 

A-36  -  A-90 

TABLE  A- 3 

Average  Reference  values 


Test  Set 

Runs 

• 

m 

(lb  /sec) 

a 

U 

(ft/sec) 

Po 

(psig) 

1 

1-3 

.577 

0 

70 

0 

2 

1-3 

.494 

1.0 

30 

0 

3 

1-5 

.807 

1.0 

49 

0 

4 

XXI 

.652 

0 

33 

-0.010 

4 

XX2 

.652 

0.47 

33 

0.015 

4 

XX3 

.652 

1.0 

33 

0.016 

4 

XX4 

.652 

2.13 

33 

0.014 

4 

XX5 

.652 

3.91 

33 

0.008 

Reference  Density  *  0.079  lbm/ft3 
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Figure  A-l 


Figure  A-2  Dimensionless  Axial  Velocity,  Test  |i 


Figure  a-3  Dimensionless  Tangential  velocity.  Test  fl 


Figure  A-4  Dimensionless  Velocity,  Test  #1 


Figure  A-5  Dimensionless  Total  pressure.  Test  fl 


Figure  a-6  Dimensionless  Static  pressure.  Test  fl 


Dimensionless  Ax 


Figure  A-8  Dimensionless  Tangential  Momentum  Flux 
Test  |1 
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Figure  A-10  Dimensionless  Axial  velocity.  Test  |2 
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Figure  a-26  Dimensionless  velocity.  Test  |3 
Five  Axial  Locations 


Three  Axial  Locations 


-29  Dimensionless  Total  pressure.  Test  |3 
Three  Axial  Locations 


4-50 


Five  Axial  Locations 
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Figure  A-31  Dimensionless  Static  pressure.  Test  #3, 
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|3,  Five  Axial  Locations 
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Figure  A-33  Dimensionless  Axial  Momentum  Flux 
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Figure  A-43  Dimensionless  Total  pressure.  Test  #4 
Runs  XXI 


Figure  A-45  Dimensionless  Axial  Momentum  Flu 
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Figure  A-50  Dimensionless  Axial  velocity.  Test  |4 
Runs  XX2 
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Figure  A-54  Dimensionless  Total  pressure.  Test  |4 
Runs  XX2 
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Figure  A-61  Dimensionless  Axial  velocity.  Test  #4 
Runs  XX3 
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Figure  A-64  Dimensionless  Velocity,  Test  |4,  Runs  XX3 
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Figure  A-70  Dimensionless  Static  pressure  on 
6  Inch  Tube,  Test  |4,  Runs  XX4 
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Figure  A-75  Dimensionless  Velocity,  Test  #4,  Runs  XX4 


Fiyure  A-76  Dimensionless  Total  Pressure,  Test  |4 
Runs  XX4 
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Figure  A-78  Dimensionless  Axial  Momentum  Flux 
Test  #4,  Runs  XX4 
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Figure  A-82  Dimensionless  Static  pressure  on 
8  Inch  Tube,  Test  |4,  Runs  XX5 


Figure  A-86  Dimensionless  velocity,  Test  #4,  Runs  XX5 
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APPENDIX  B 


DEVELOPMENT  OF  A  MODEL  OF  IDEAL  INCOMPRESSIBLE 
MIXER  WITH  SWIRL 


I.  INTRODUCTION 

As  a  model  of  an  ideal  mixer  with  swirl,  consider 
inviscid  incompressibleone-dimensional  flow  in  the 
constant  area  mixer  shown  in  Figure  2-1.  The  outer  stream 
(2)  enters  the  mixer  at  a  uniform  axial  velocity,  w2 ,  and 
uniform  total  pressure,  pT  .  The  inner  stream  (0)  enters 
the  mixer  at  a  uniform  axial  velocity,  wQ ,  and  uniform 
total  pressure,  pT  .  The  inner  stream  then  passes  through 
a  free  vortex  stator  before  entering  the  mixing  region. 
The  mixed  stream  (3)  is  assumed  to  have  uniform  axial 
velocity,  w3  ,  and  an  angular  velocity  profile 
corresponding  to  solid  body  rotation.  The  mixed  stream 
then  passes  through  a  stator  that  removes  this  solid  body 
rotation  and  the  departing  mixed  stream  (4)  has  uniform 
static  pressure,  p4  ,  and  no  angular  velocity.  This 
departing  stream  (4)  is  then  accelerated  through  a  nozzle 
to  an  exit  static  pressure  of  p#  . 


Each  section  o£  the  model  will  be  considered  separately 
in  order  to  develop  the  equations  needed  for  analysis. 
After  all  sections  of  the  ideal  swirl  mixer  are  analyzed, 
the  equations  for  thrust,  using  an  ideal  constant  area 
mixer  and  the  thrust  for  unmixed  flow,  will  be  developed. 
Then  the  equations  for  comparing  the  performance  of  the 
ideal  swirl  mixer  to  the  ideal  constant  area  mixer,  or  to 
the  unmixed  flow  thrust,  will  be  listed  in  the  order  used 
for  analysis. 

II.  UNIFORM  ENTERING  STREAMS 


Both  the  inner  stream  (0)  and  the  outer  stream  (2)  enter 
the  ideal  mixer  with  uniform  axial  velocity  and  total 
pressure.  The  following  equations  apply  to  these  two 
streams: 


A) 

Total  pressure 

%  - 

P0  +  1/2 

po  wl 

(B-l) 

PT2  * 

P2  +  1/2 

p2  w22 

(B-2) 

B) 

Mass  Flow  Ratio 

(  a  ) 

• 

nig 

n  ■  _ s 

P2W2  A2 

/  > 

m0 

Rd  wo  Ao 

) 

C) 

Area  Ratio 

a2 

Rt2  -  rt2 

g2  —  f2 

(  A  \ 

A0 

r  2  _  -2 
rt  rh 

f2-  1 

where 

f  -  rt  /rh 

and  g 

•  R»/rh 
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III.  FREE  VORTEX  STATOR  AND  DEPARTING  SWIRLING  STREAM 


For  an  ideal  free  vortex  stator,  the  total  pressure  is 

constant  on  streamlines.  Because  the  entering  stream  has 

uniform  axial  velocity  (wQ)  and  uniform  total  pressure 

(PT  ) ,  the  departing  stream  will  have  uniform  axial 

velocity  and  uniform  total  pressure.  Thus,  the  entering 

and  departing  total  pressures  are  equal  (pT  »  pT  ) . 

1 1  '0 

Also,  the  entering  and  departing  axial  velocities 
are  equal  (w1  *  wQ)  because  the  flow  is  incompressible 
and  the  area  is  constant.  The  following  equations  apply 
to  the  stream  departing  the  free  vortex  stator: 


A)  Total  pressure 


PTi  -  P*  +  l/2P0(w02  +  v*) 


B)  Angular  velocity  and  Static  pressyre 

vi  -  vi, 

P,  -  Pd  -  J,/2  e„v,2  (rjVr2-  l) 
where  vt^  is  the  value  of  v,  at  r=rt 
and  Pd  is  the  value  of  Pt  at  r=  r, 


(B-5) 


(B-6) 

(B-7) 


The  axial  force  (F,  )  on  the  free  vortex  stator  in  the 
constant  area  annulus  is  determined  by  considering  the 
actuator  disc  limit  axial  force  on  a  blade  row  given  by 


331+ 


v 


Equation  45  in  Reference  l  and  listed  below. 


E,  =  2k  [m#] 


(B-8 ) 


where 


m'-  l(H  - 

*4. 


P0/rt  +  1/2  (w2-v2-u2  ) )  r  dr 


P0  =  Reference  pressure 
H  ■  Stagnation  enthalpy 
and  [  |=  Change  across 


The  axial  force  (F,  )  on  the  free  vortex  stator  using 

Equation  B-8  and  assuming  radial  equilibrium  is  given  by 


or  F  *  -x  pQ  v,J  rt2  In ( f } 


(B-9) 


IV.  MIXED  STREAM 


The  mixed  stream  is  assumed  to  have  uniform  axial 
velocity,  solid  body  angular  velocity  and  no  radial 
velocity.  The  angular  velocity  and  static  pressure  are 
given  by 


\  r  /R, 

(B-10) 

V  +  1/2  °3V>3  (r2/R*>  '  11 

3 

( B— 11) 

where  vt  is  the  value  of  v3  at  r=Rt 

and  p_  is  the  value  of  p,  at  r=R, 

3  1 

The  density  of  the  mixed  stream  can  be  written  in  terms 
of  the  density  of  the  two  entering  streams  and  the  mass 
flow  ratio  (  a )  since  the  volume  flow  rate  is  constant. 

8 

mQ  +  m2  1  +  a 

p  =  - : - -  -  (  b— 1 2 ) 

m0/p0  +  m2/p2  1  +  ap  /p2 

The  axial  velocity  can  be  written  in  terms  of  the 
properties  of  the  entering  streams  by  application  of  the 
conservation  of  mass  to  the  one-dimensional  constant  area 
flow. 
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m0  ^  —  (A^  ^  A2 )  p,,  w. 


0  T  2 '  ^3  3 


or  w3  »  wQ 


1  +  a  Pn/P2 


1  +  a 


( B— 1 3 ) 


A)  Axial  Momentum  Equation 

Application  of  the  axial  momentum  to  the  mixing 

region  gives  the  following  relationship  which  can  be 

solved  to  give  pB  in  terms  of  other  properties 

*3 

/P1  dA  +  P0W0A0  +  (P2  +  P2W2  )A2  =/P3dA  +  P3W32(A0+A2) 


(B-14) 


The  first  integral  is  given  by 


yp,  dA  »  (Pd+pov,V2)  A  -xpov,2rt2ln(f)  (B-15) 

i 

and  the  second  integral  is  given  by 
/p 3dA  -  <PRtJ  P3  vt23/2)  (A0+A2)  +*p3  ~  — -  (B-16 ) 

Substitution  of  Equations  B-15  and  B-16  into 
Equation  B-14  give  the  following  relationship  for  pR 


ill,  poWo  poV' 


f2-l 


2 

+  P3vt2(l-h2  )  / 4 


In ( f ,  +  — «  p  + 


(B-17) 


where  h  *  rh/Rt  a  l/o 
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B)  Conservation  of  Angular  Momentum 

Application  of  the  conservation  of  angular 

momentum  to  the  mixing  region  gives  the  following 

relationship  which  can  be  solved  to  give  v,  in  terms  of 

*3 

upstream  properties 


/p„wovi  r  dA  "  /p3w3v3rdA 


(B-18) 


a0+a2 


Substitution  of  Equations  B-6  and  B-10  into 
Equation  B-18  and  integration  gives  the  following 


relationship  for  v,^ 


‘3  =  Z  V»1 


Pp  w0  rt  Rt  (rt2  -rh2  ) 
p2w2  r  f  -  rh4 


or  v. 


(B-19) 


1+  a  gz  +  1 


V.  SOLID  BODY  ROTATION  STATOR  AND  DEPARTING  STREAM  (4) 


The  ideal  stator  between  stations  3  and  4  is  assumed  to 
remove  all  angular  velocity  from  the  mixed  stream,  since 
the  flow  is  assumed  to  be  inviscid,  the  torque  on  this 
stator  is  equal  to  the  torque  on  the  free  vortex  stator 
but  opposite  in  direction.  Thus,  the  net  torque  on  the 
mixer  is  zero.  To  determine  the  axial  force  on  this 
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stator,  the  axial  velocity  profile  of  the  departing 
stream  (4)  is  required.  Since  the  departing  stream  (4) 
has  no  swirl  and  is  in  radial  equilibrium,  the  static 
pressure  (p4)  will  be  uniform. 

A)  properties  of  Departing  stream 

The  stream  function  (  ijj)  is  defined  as  follows  in 
terms  of  the  axial  velocity  of  the  mixed  stream  (w3 ) . 


dMJ=  - 

r  w3 

dr 

Since  w3 

is 

constant. 

the  above  equation 

integrates  to  give 

the 

radius  in 

terms  of  the  stream 

4 

function. 

r2  ■  R*  +2  (  t|lt  -  MO  /  w3 

where  M»t  is  the  value  of  the  stream  function,  M>, 
at  r  =  R, 

Thus  the  angula  velocity  of  the  mixed  stream, 
given  by  Equation  B-10,  can  be  written  in  terms  of  the 
*  stream  function,  M*  ,  to  give 

=V'3  [  1  " 

Now  the  stagnation  enthalpy,  h  ,  can  be  written 
in  terms  of  the  stream  function  to  give 


33  9 


H-  (ip  )  =  K  - 


W,Rt 


(B-20 ) 


PR  w3z  v{2 

where  K  is  a  constant  (=  — *3+  —  +  — 2) 

P3  2  2 

The  stagnation  enthalpy  is  constant  on  a 
streamline  across  the  ideal  stator 

H3(<4M  =  H4  (ty  )  (B-21 ) 

and  we  have  from  Reference  2  for  radial 
equilibrium 


dH3  dH4 


( B — 22) 


dV  d4* 

where  Equation  B-22  can  be  evaluated  for  the 
mixed  stream  by  using  Equation  B-20  to  get 


dH3 

2  V‘3 

/ R— 0 1 ^ 

d<v 

W3R,2 

VD— Z J ) 

The 

axial  velocity  of 

the 

departing  stream,  w4  , 

is  related 

to  the  stagnation 

enthalpy  and  the  stream 

function  for 

radial  equilibrium 

flow 

with  no  swirl  by 

1  dw4  dH4 


(B-24) 


r  dr  d'V 

Substitution  of  Equation  B-23  into  Equation  B-24 


and  integration  gives 


V 


vt23  r* 

w4  -  —  -r  +  c 

W3  Rt2 


( B— 25 ) 


The  constant,  C,  in  the  above  equation  is 
evaluated  by  application  of  the  conservation  of  mass 
between  stations  3  and  4  to  give  the  following  resulting 
relationship  for  the  axial  velocity  w4. 


vu  *■  2 

w4  =  w3  +  1/2  ~  2  — — -(1+h  ) 

w,  Rt 


(B-26) 


The  value  of  the  static  pressure  at  the  departing 
stream,  P4  ,  can  be  determined  by  evaluating  H3  at  r=Rt  ,  H4 
at  r=»Rt  and  then  equating  h3  and  h4  to  get 


*  V  +  1/2  °3  l“32  +  V,3  '  W'4> 


(B-27) 


where  wt  =  w.  +  1/2  —  (1-h2) 


(B-28) 


B)  Axial  Thrust  on  Solid  Body  Rotation  Stator 

The  axial  thrust  on  the  stator  is  given  by 
Equation  B-8,  rewritten  below  with  the  upper  limit  of 
integration  changed  to  Rt 


F2  =  2*s[m;] 


(B-29) 


31+1 


where 


M3  “  /{(H4+  w42/2)  “(H3+  (w32-v32  )/2  ))r  dr 


Using  Equations  B-10 ,  B-ll,  B-26,  B-27  and  B-28, 
the  two  terms  within  the  integral  can  be  written  as 


h  +  <W3*  -1V3)/2  =  VP3+  ^  +  V»3  (r2/Rt2"  1} 


H4  +  w42/2  =  PR/p  +  W32  +  Vt2  (r2/R2  -  l)/2  +  G 
t3  J  3 

vt2  r  3  r2  I  v  4  r  r2 

where  G  =  — - - (1+h2)  +  -2  2(2 - 1-h2  )  -(1-h2)2 

2  L  Rt  J  8w32  Rt 


(B-30) 


Substitution  of  the  above  relationships  into 


Equation  B-29  gives 


F2  =*  2  *  p 3  I  G  r  dr 


Integration  of  the  above  equation  for  F,  gives 


the  following  relationship 


P3*v»3Rt  U-h2) 


6  (1+h2  )  -  M-  (1-h2)5 


(B-31) 


IV.  NOZZLE  THRUST,  INCOMPRESSIBLE  FLOW 


For  ideal  incompressible  flow  through  a  nozzle  exhausting 
to  a  static  pressure  of  p#  ,  the  velocity  of  the  exiting 
stream  is  given  by 


«.  -  {2  (H4  - 


(B-32) 


Since  h4  varies,  the  above  equation  will  be  applied  to 
each  streamline  and  the  result  used  in  calculating  the 
nozzle  thrust,  FN  .  The  stagnation  enthalpy,  H4  ,  can  be 
expressed  in  terms  of  a  stream  function,  V,  by  defining 


1  dvp 


w4  -  -  “ 


r  dr 


Substitution  of  Equation  B-26  into  the  above  equation  and 


integration  gives 


[yZ  «  yZj^Z 

w3 - ^(1+h  )(Rt2/2)(l-ra/R?)  +  -^-(l-r4/Rt4) 

2w3  JV  4  w3 

where  4*  is  the  value  of  the  stream  function  at  r=R. 


Solving  the  above  equation  for  (r/Rt )  and  substituting 
into  Equation  B-26  yields 


w4  <4») 


v2 

1+  — ^(1-h2) 
2w3 


4Vt2 

w|  Wj  R( 


Using  the  above  equation  for  w4 ,  the  stagnation  enthalpy, 
H4  ,  can  now  be  expressed  in  terms  of  the  stream  function 
as 


H4  (v) 


(1-h2) 


4v2  V-qit 
w|  w3  R,2 


(B-33 ) 


Ndw  substituting  Equation  B-33  into  Equation  B-22  gives 
the  nozzle  exit  velocity,  wa  ,  in  terms  of'  the  stream 
function, 

1/2 

(B-34 ) 

The  mass  flow  rate  for  incompressible  flow  can  be  written 
in  terms  of  the  stream  function  as 


w#  ( m)  ■ 


P4-P. 


r  „  12 

1+  — —  ( 1— h2  )  - 

L  2w2  J 


2  4v,2  4>-4)t 


W3  W3Rt 


dm  =  -  2  x  p3  d  4/ 


(B-35) 


The  nozzle  thrust  is  equal  to  the  momentum  flux  of  the 
exiting  stream 


Using  Equation  B-35,  FN  can  be  written  as 


A 

Fn  »  -  2dp3  w,  d4*  (B-36 ) 

% 

Defining  i|»  »  0,  then  4*  =  w,R?  (l-h2)/2 

t  n  o  I 


Defining  a  new  variable  <p  -  2  4*  vt2/(w3R2)  (B-37) 

Then  <p  =  0,  <p  -  (l-h2)v2  and  d<p  -  2  vtV(w3Rt2)  d4* 
Rewriting  Equation  B-34  in  terms  of  the  new  variable 
gives 


we  (<P) 


P4-P. 

2 -  + 


«?fl*  — (l-h2)]2  -  20 
L  2w»  J 


1/2 


(B-38) 


Substituting  Equations  B-37  and  B-38  into  Equation  B-36 


yields  the  following  equation  for  FN 
F  =  {I  p3  *  R2  { 


(l-h2)vt2 


[w3/vt23)  /  (D  -<p 

Jq 


-<P)  dip 


(B-39) 


where  D  ■  [w#  (<j>t )]  /2 


(B-40 ) 


Integration  of  Equation  B-39  gives  the  following  result 
for  the  nozzle  thrust 


W.  (<Pt)]  -  [we  «|>h ) 


where  w#  (<J>  )  =  2 -  +  w,  1  +  - 2(i-h4) 

‘  (  ^3  3  I  2wf 

!P  — p 

2  -1..  6  +  w2  |x - l2(i-h2  ) 


(B-42) 


(B-43) 


VII.  NOZZLE  THRUST/  COMPRESSIBLE  FLOW 


For  isentropic  flow  of  a  perfect  gas  through  a  nozzle 
exhausting  to  a  static  pressure  of  p#  ,  the  velocity  of 
the  exiting  stream  is  given  by 

w.  =  2,v(h4  -  P./p.)/(7-l)  (B-44) 

Substituting  Equation  B-33  into  the  above  equation  and 
defining  <J>  by  Equation  B-37  we  get 


Substituting  the  above  equation  into  the  thrust  equation. 
Equation  B-36,  and  integrating  gives  the  following  result 
for  the  thrust  for  compressible  flow  through  a  nozzle: 


where  w#  «v 
and  w#  (<J>h) 


2y  IP4 
Jv-1  (p3 


P  W* 

-2+  ^ 

P.  2 


TT  p4  p#  w| 

Jv-1  P3  p.  2 


1+  — *(l-lf ) 

2w? 


1+  — !2(i-h2) 

2w | 


1/2 


(B-46 ) 


1/2 


(B-47) 


VIII.  TOTAL  THRUST  FOR  SWIRL  MIXER  WITH  NOZZLE 


The  total  thrust  for  the  swirl  mixer  with  nozzle,  FSM  ,  is 
the  sum  of  the  thrust  of  each  of  the  two  stators  and  of 
the  thrust  of  the  nozzle. 

fsm  -  F1  +  F2  +  fn  (B-48) 

where  E,  is  given  by  Equation  B-9 
F2  is  given  by  Equation  B-31 

FN is  given  by  either  Equation  B-41  or  Equation  B-45 

IX.  CONSTANT  AREA  MIXER  WITH  NOZZLE 

The  performance  of  constant  area  mixer  with  nozzle  will 
be  developed  in  this  section  for  comparison  with  the 
performance  of  the  swirl  mixer  with  nozzle.  The  physical 
configuration  of  the  constant  area  mixer  with  nozzle  is 
the  same  as  shown  in  Figure  1  without  the  two  stators. 
Thus,  the  flow  properties  are  unchanged  from  station  "0" 


to  station  "1"  and  from  station  "3"  to  station  "4".  The 
Kutta  condition  at  the  splitter  plate  requires  that  the 
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two  entering  static  pressures  be  equal  (PQ  =  P2 )  . 
Equations  B-l  through  B-4  ,  B-12  and  B-13  apply  for  this 
case.  Using  Equations  B-l  and  B-2  and  the  condition 
that  pQ  =  P2  gives  the  following  equation  for  the  axial 


velocity,  w, . 


w2  = 


[<W  +  2  <pt,  -  pt0  )]/p, 


1/2 


(B-49) 


Application  of  the  axial  momentum  equation  across  the 

mixing  region  gives  the  following  equation  for  PT  . 

3 


pT- 


n  w' 
0 


(1+a)  (1+ao/f^) 


-1 


2  a  P0W0  /  (P2W2  ) 
1+«P0W0/(P2W2) 


<W 

(B-50) 


The  thrust  for  ideal  incompressible  flow  through  a  nozzle 
with  uniform  entering  conditions  is  given  by 

F st  58  (m0+*2)w0Jp7  (B-51) 


where  p  =  wz/w2  =  2PT  (1-  P,  /PT  )  /  (  P3  w02) 


(B-52) 


The  thrust  for  ideal  compressible  flow  through  a  nozzle 
with  uniform  entering  conditions  is  given  by  Equation 
B-51  but  with  p.  given  by 


X.  UNMIXED  FLOW  THRUST 


The  thrust  for  two  unmixed  streams  flowing  ideally  and 


incompressibly  through  separate  nozzle  to 

exhaust  pressure,  p#  ,  is  given  by 

the  same 

n 

z 

6? 

mQ+  ra2)w0(>/^1+aN/^)/(i+a  ) 

(B-54 ) 

where 

2(V  '’•)/(powo> 

(B— 55 ) 

and 

f>2=  2lPT2-  p.  >/(  P2W1> 

(B-56) 

The  thrust  for  two  unmixed  streams  flowing  ideally  and 
compressibly  through  separate  nozzles  to  the  same  exhaust 


pressure,  p#  ,  is  given  by  Eqn  B-54  and 

Y  (ptn- 

'  Y-l  p  w2/2 


(B-57 ) 


PT2-(P«)Y  (PT2-^W2/2) 
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XI.  WORKING  EQUATIONS  FOR  SWIRL  MIXER  PERFORMANCE 
COMPARISON 


The  equations  developed  in  the  previous  sections  are 
presented  in  this  section  in  a  form  which  facilitates 
computer  calculations.  The  ratio  of  the  swirl  velocity 
leaving  the  free  vortex  stator  to  the  axial  velocity  at 


r=rt  r  vt, 

/Wo 

,  is  an 

input 

variable  that 

is  used  to  vary 

the  deg  re 

of 

swirl 

during 

a  series  of 

calculations 

to 

determine 

the 

swirl 

mixer 

performance 

over  a  range 

of 

swirl.  This  velocity  ratio  is  defined  by 

St|  =  v,  /w0  ( B— 59 ) 

The  ratio  S,  has  a  minimum  value  when  p_  ,  PT  ,  p„  and  w0 

1  '0  '2  0 

are  defined  because  the  modified  Kutta  condition, 
value  of  P2  /Pd  ,  at  the  splitter  plate  and  the  velocity 
of  stream  2,  w2  ,  must  be  real.  For  w2  =»0 ,  we  get  from 
Equations  B-2,  B-5  and  B-7  the  following  relationship  for 
the  minimum  value  of  St 


PTo  -  PT;(Pd  /p:  ) 
Po  wo  /  2 


(E  -60) 


S,  has  a  minimum  value  of  zero  for  values  of  (S,  )  , 

t  1  min 

less  than  or  equal  to  zero. 
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The  input  variables  were  selected  to  permit  comparison  of 
mixers  and  are  listed  below.  Also  listed  are  the  output 
variables  and  the  equations  in  the  order  of  calculation 


INPUT  VARIABLES:  PTo ,  p0  ,  w0  ,  ,  P2  ,  a  ,  g ,  S,  ,  P#  /P^  , 

Pd  /P,  /  &  Y 

OUTPUT  VARIABLES:  w  /w  ,  A  /A  ,  f,  w  /w  ,  v,  /w  , 

F^/Jm  +  m  )w  >  F2/(m  +  m  )w  ,  and 
the  thrust  ratios  F^/F^,  , 

W  &  fstc/eumc 


WORKING  EQUATIONS 


(S,2) 


min 


PT;  (Pd  /P»  ) 

PoWoJ/2 


1 


if  (s,2 )  ,  is  less  than  zero,  then  St  =  0 
‘  min  ‘min 


A.  Swirl  Mixer  S,  >  St  >0 

*  ~  lmin 

Pd  =  PTo  "  Pa  W0J(l+StJ)/2 
P,  =Pd/(Pd/P3) 
w3  =  [  2(PTj  -  P,  )/P2] 1/2 
A,  /A0  *  aP0  w0  /  (Pj  Wj  ) 

f  =  (g2+  A j /A0  ) / (1  +  A,/A0  ) 
p3/p  =  (1  +  a)/(l  +  a  P0/P,  ) 

Wj  /w  *  (1  +ap0  /p2  )/(l  +  Aj/A0) 


h  =  1/g 


V.  /vt  =  2fg/(l  +  a)  (l  +  h2  ) 

l3 

vt/w3  =  S,  (vtj/vt  )  (w3  /w0  ) 

Wt</w3  ■  1  +  (v,3/w3  )J  (1  -  h5)/2 


(m0  +m3  )  w0 


f;  In (f ) 
(1+a)  (f'-l) 


(m.  +m 


F„  w3  /  vt  \2  1  T  /vt  \J 

- -  =  —  (—  )  —  6  ( 1+h2  )  -(  — )  (1-h2  )2 

^2)w0  w0  \w3  /  24  L  \w3/ 


(m0  +m3  )w0 


Wo  V  PWJ  /2 


.  </2  p  +  3  ,  »»♦«>  %y  «■/ 

(1+a)  (w3/w0  )  \p0w0J  (m0  +m3  )w  J  w„  y  pvi 

+  (  p3  w3J  / 2 )  j  (vtj/w3  )2  (3-h2  )/2  -  1  -  (wt/w3  )2J 

1.  Nozzle  with  Incompressible  Flow 


^-+  1 


w-t  = 


w-h  = 


2(P4  - 

P.  ) 

/  vt3 

\' 

+  w  J  (l  +  — 

(1-h2  ) 

03 

V  2w3j 

/ 

2(P4  - 

P.  ) 

,  /  vt’ 

,  V 

+  w3J  1 - 

(1-h2  ) 

03 

V  2w3 

) 

(mo  +m2  )w0 


rSM 

(m0  +m  )  w 


w>t  -  < 

3w0  v,’  (1-h2) 


Ej  +  F2+  Fn 
(m0  +ms  )w0 


2.  Nozzle  with  Compressible  Flow 


w„ 


hr— i 


2p<  r  /p.\¥i  2  r  t 

—  1  -  -  +  w32  1  +  — 1 

1  p,  L  W  j  L  2w,J 


w6h=| 


Y 

Y-l 


r  /p.\¥i  , r,  Vi> 
—  i  -  —  + «.  i  — - 

o  L  \p/  J  L  2w32 


Y-l 


w3  -  w3 

W«t  W«h 


(m0+m2)Wo  3  wQvtJ(l-hJ) 


rSMC 


R  +  F,  +  Fn 


(m0+m2)w0  (m0+mJ)wQ 


B.  Constant  Area  Mixer  (St  =  0) 


P,  =  Ph 


w2  = 


2(Pt-  P,+ P0wo2/2)/P2 


l'/» 


PT  =  PT  - 

'3  l0 


p  w^ 

0  0 


(1+a)  (1+ap/P,  ) 


2  L  (l+aP,w0/(  P2  w,  » 

2  a  p0w0  /  (  p,  w,  ) 


1  +ctp0wo  /(  p,w2  ) 


(PT- 

*  n 


Nozzle  with  Incompressible  Flow 


2.  Nozzle  with  Compressible  Flow 


STC 


P3  ■  ^PTj  -  (P.  &  (PTj  -  P3w’/2)1x^/(  P3  we*/2) 

■v/57 


(m0  +m3 )  wa 
C.  Unmixed  Flow  Thrust  (St =  0) 


1.  Nozzle  with 

/ 

Incompressible  Flow 

*UM  1 

PT0  -  P«  . 

h  -  f. 

(*„+*,  )w0  1+0 

Powol/2  1 

P3  w32  /2 

2.  Nozzle  with  Compressible  Flow 

0-—  [pT  -  (P.JY1  (PT  -  p0w’/2)VY  ]/(  P0 w02 / 2 ) 
A  Y-l  l  0  0  J 

[pTi  -  <P,  )V  (PTj  -  P,w1’/2)'/V  ]/(  P,w„,/2) 

w  . 

(m0+m3  )w0  1  +  « 
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APPENDIX  C 
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CONSTANT  AREA  MIXER  WITH  SWIRL  COMPUTER  PROGRAM  -  "SWIRL1 


This  appendix  contains  a  listing  of  the  computer  program 
SWIRL.  This  program  uses  the  equations  developed  in 
Appendix  B  and  listed  at  the  end  of  that  appendix.  Table 
C-l  contains  a  list  of  the  variables  used  in  program 


SWIRL. 


TABLE  C-l 


Variables  used  in  Program  SWIRL 

A)  INPUT  VARIABLES  (listed  in  the  order  used) 

NDR  -  Number  of  density  ratios  for  this  run 
RHOZ  -  Density  of  stream  "0" 

PTZ  -  Total  pressure  of  stream  " 0 " 

DRZ2  -  Ratio  of  the  density  of  stream  ”0"  to  the 
density  of  stream  ”2” 

RTHR  -  over  11  tip  to  hub  ratio  of  mixer 
SN  -  Maximum  tip  swirl  ratio 

PREZ  -  Ratio  of  static  pressure  of  exhaust  stream 
to  the  total  pressure  of  steam  *0* 

SPRD2  -  Ratio  of  the  static  pressure  at  point  "d" 
to  the  static  pressure  in  stream  ”2" 

GAM  -  Ratio  of  specific  heats 

NW  -  Number  of  axial  velocities,  w,  of  stream  "0" 

NAL  -  Number  of  mass  flow  ratios,  AL 

NP  -  Number  of  total  pressure  ratios,  p 
W  -  Axial  velocity  of  stream  "0" 

AL  -  Ratio  of  the  mass  flow  rate  of  steam  ”2” 

to  the  mass  flow  rate  of  stream  "0". 
p  -  Ratio  of  the  total  pressure  of  stream  ”0” 
to  the  total  pressure  of  stream  ”2” 

B)  OUTPUT  VARIABLES  (listed  in  the  output  order) 

PT2  -  Total  pressure  of  stream  ”2” 

RH02  -  Density  of  stream  " 2" 

1)  CONSTANT  AREA  MIXER  WITHOUT  SWIRL 
PDP  -  Static  pressure  at  point  "d" 

P2P  -  Static  pressure  of  stream  ”2” 

W2P  -  Axial  velocity  of  stream  *2" 

A2Az  -  Ratio  of  cross-sectional  area  of  stream  "2" 
to  that  of  stream  *0" 

W3  -  Axial  velocity  of  mixed  stream 
F  -  Tip  to  hub  ratio  of  inner  annulus 
FUM  -  Thrust  of  unmixed  flow  with  incompressible 
flow  through  the  nozzles 

FUMC  -  Thrust  of  unmixed  flow  with  isentropic  flow 
through  the  nozzles 

PT3P  -  Total  pressure  of  the  mixed  stream 
PTAM  -  Mass  average  total  pressure  of  the  two 
incoming  streams 

RD2  -  Location  of  dividing  streamline  at  the 
entrance  to  the  mixer 
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RD3  -  Location  of  dividing  streamline  at  the 
exit  of  the  mixer 

PSTUM  -  Ratio  of  the  thrust  of  mixed  flow  with 
incompressible  flow  through  the  nozzle 
to  the  thrust  of  unmixed  flow,  FUM 
FSTUMC-  Ratio  of  the  thrust  of  mixed  flow  with 
isentropic  flow  through  the  nozzle  to 
the  thrust  of  unmixed  flow,  FUMC 


2)  CONSTANT  AREA  MIXER  WITH  SWIRL  (values  of  the 
following  for  each  iteration  of  s) 

S  -  Tip  swirl  ratio  leaving  stator  "1* 

VT3W3  -  Tip  swirl  ratio  of  mixed  stream  at  ”3" 

W2WZ  -  Ratio  of  the  axial  velocity  of  stream  "2” 
to  that  of  stream  ”0” 

W3WZ  -  Ratio  of  the  axial  velocity  of  stream  "3" 
to  that  of  stream  "0" 

A2AZ  -  Ratio  of  the  cross-sectional  area  of 
stream  "2"  to  that  of  stream  *0" 

F  -  Tip  to  hub  ratio  of  inner  annulus 

FI  -  Thrust  on  stator-1 

F2  -  Thrust  on  stator-2 

FN  -  incombressible  nozzle  thrust 

FNC  -  isentropic  nozzle  thrust 

FSWFST  -  Ratio  of  the  thrust  of  swirl  mixed  flow 

with  incompressible  flow  through  the 
nozzle  to  that  of  mixed  flow  without 
swirl 

FSWSTC  -  Ratio  of  the  thrust  of  swirl  mixed  flow 
with  isentropic  flow  through  the  nozzle 
to  that  of  mixed  flow  without  swirl 
FSW  -  Thrust  of  swirl  mixer  with  incompressible 
flow  through  the  nozzle 

FUM  -  Thrust  of  unmixed  flow  with  incompressible 
flow  through  the  nozzles 

FSWUM  -  Ratio  of  the  thrust  of  swirl  mixed  flow 
with  incompressible  flow  through  the 
nozzle  to  that  of  unmixed  flow 
FSWUMC  -  Ratio  of  the  thrust  of  swirl  mixed  flow 
with  isentropic  flow  through  the  nozzle 
to  that  of  unmixed  flow 
SW  -  Swirl  number 
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APPENDIX  D 

INITIAL  DATA  REDUCTION  PROGRAM  -  "REDUCE" 


This  appendix  contains  a  listing  of  the  computer  program 
REDUCE  as  described  in  Chapter  6.  Table  D-l  is  a  list  of 
the  variables  used  in  the  program  REDUCE ,  Table  d-2  is  a 
partial  listing  of  a  RAW  data  file,  and  Table  d-3  is  a 
partial  listing  of  a  FINE  data  file. 


TABLE  D— 1 


Variables  used  in  program  REDUCE 

A)  INPUT  VARIABLES  (listed  in  the  order  used) 

NSRUNS  -  Number  of  subruns  within  this  run 
(the  following  is  input  for  each  subrun) 

NSUB  -  Subrun  number 

NSCANS  -  Number  of  data  scans  for  this  subrun 
(the  following  is  input  for  each  scan  "K") 

NHR  -  Hour  of  scan 

NMIN  -  Minute  of  scan 

NSEC  -  Second  of  scan 

DATA(l)  -  Scanivalve  position  from  controller  voltage 

staircase 

DATA ( 2 )  -  Alpha,  voltage  corresponding  to  the  angle  of 

the  probe 

DATA ( 3 )  -  Voltage  from  supply  line  pressure  transducer 

DATA (4)  -  Voltage  from  supply  tank  pressure  transducer 

DATA ( 5 )  -  Voltage  from  LVDT,  radial  position  of  probe 

head 

DATA ( 6 )  -  Temperature  of  inside  plenum  (c) 

DATA (7)  -  Temperature  of  outside  plenum  (C) 

DATA ( 8 )  -  Voltage  from  P(l)  pressure  transducer 

DATA ( 9 )  -  voltage  from  P(2)  pressure  transducer 

DATA (10)  -  Voltage  from  P(3)  pressure  transducer 

DATA (11)  -  Voltage  from  P(4)  pressure  transducer 

DATA (12)  -  voltage  from  p(5)  pressure  transducer 

DATA (13)  -  voltage  from  P(l)-P(3)  pressure  transducer 

DATA (14)  -  Voltage  from  outside  flow  orifice  plate 

differential  pressure  transducer  . 

DATA (15)  -  Voltage  from  inside  flow  orifice  plate 

differential  pressure  transducer 
VPRESA(K)  -  Voltage  of  scanivalve  pressure  transducer  A 
VPRESB(K)  -  Voltage  of  scanivalve  pressure  transducer  B 


B)  OUTPUT  VARIABLES  (listed  in  the  order  used) 

L  -  Subrun  number 

MHR  -  Hour  of  subrun 

MMIN  -  Minute  of  subrun 

VPRESA(l)  through  VPRESA(49)  -  Input  value  less  the 
transducer  zero  value 

VPRESB(l)  through  VPRESB(49)  -  Input  value  less  the 
transducer  zero  value 

DATA ( 2 )  through  DATA (15)  -  Average  of  scan  #6  through 

scan  #45  (transducer  zero  has  been  subtracted 
for  DATA (8)  through  DATA (15)) 
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APPENDIX  E 


EXPERIMENTAL  DATA  REDUCTION  PROGRAM  -"DATAR4 1 


This  appendix  contains  a  listing  of  the  computer  program 
DATAR4  as  described  in  Chapter  6.  Table  E-l  contains  a 
list  of  the  variables  used  in  this  computer  program. 
Table  E-2  is  a  list  of  the  file  LAYOUT  that  is  the  coded 
connection  of  the  wall  static  pressure  taps.  Table  E-3  is 
a  sample  output  of  the  program  DATAR4  for  run  #061. 


d 


TABLE  E-l 


Variables  Used  in  program  DATAR4 
A)  INPUT  VARIABLES  (From  Files  #1,  12,  #3  &  #4) 

1)  From  File  #1  (Layout  o£  pressures  on  scanivalve) 

NTAPS  -  Number  of  static  pressure  taps  used 

N(l)  -  Number  of  static  pressure  taps  on  4”  OD  tube 

N(2)  -  Number  of  static  pressure  taps  on  6”  OD  tube 

N (3)  -  Number  of  static  pressure  taps  on  8”  ID  tube 

IP(1)  through  ip(74)  -  scanivalve  port  connections  for 
the  74  possible  static  pressures (Can  have  a 
value  of  -1,  1  through  47  and  49  through  95. 

If  =-l,  then  this  static  pressure  is  not 
connected) 

2)  From  File  #3  and  #4  (probe  calibration  coefficients) 

KCX (K , J , I )  -  probe  calibration  coefficients 

I  ■  1,2  -  Corresponds  to  calibration  regions 

#1  &  #2 ,  respectively 

J  »  1,2,3, 4, 5  -  Corresponds  to  c  ,  C  ,  v/v,  u/v  & 

w/V,  respectively 

•K  ■  1  to  15 •  -  Corresponds  to  the  15  calibration 

coefficients 

3)  From  File  #2  (FINE  DATA  FILE  XXX) 


NRUN  -  Run  number 

MONTH  -  Month  of  run 

NDAY  -  Day  of  run 

NYR  -  Year  of  run 

NTCONF  -  Test  configuration  of  run 

N0RF(1)  -  Inside  choked  orifice  size 

N0RF(2)  -  Outside  choked  orifice  size 

NSVANE  -  Swirl  vane  setting 

PATM  -  Atmospheric  pressure  ("  Hg) 

PLINEI  -  Supply  line  pressure  setting  (psig) 

ALPHO  -  Calibration  probe  angle 

VALPHO  -  Voltage  from  potentiometer  for  probe  angle 
ALPHO 


RDO  -  Calibration  radial  position  of  probe  head 

VRDO  -  Voltage  from  LVDT  for  probe  position  RDO 

NTIMEH(l)  -  Hour  at  start  of  run 

NTIMEM(l)  -  Minute  at  start  of  run 

NTIMEH(2)  -  Hour  at  end  of  run 

NTIMEM(2)  -  Minute  at  end  of  run 

TRMS  -  Room  temperature  at  start  of  run  (C) 


TRME  -  Room  temperature  at  end  of  run  (C) 

PTANKS  -  Tank  pressure  at  start  of  run  (psig) 

PTANKE  -  Tank  pressure  at  end  of  run  (psig) 

NSRUNS  -  Number  of  subruns 
VP0(7)  -  Zero  of  tank  pressure  transducer 

VP0(8)  -  Zero  of  line  pressure  transducer 

VPO ( 9 )  -  Zero  of  outside  flow  orifice  plate 

differential  pressure  transducer 
VPO (10)  -  Zero  of  inside  flow  orifice  plate 

differential  pressure  transducer 
Z(l)  to  Z(6)  -  Axial  position  of  probe  section  and 
8"  ID  test  section  A,  B,  C,  D  &  E, 
respectively 

For  each  subrun  "i" 

NSRUN(I)  -  Subrun  number 

LTIMEH (I)  -  Hour  at  start  of  subrun 

LTIMEM(I)  -  Minute  at  start  of  subrun 

VPSV(l)  to  VPSV(96)  -  Scanivalve  pressure  transducer 

output  voltages 

VA  -  Voltage  from  potentiometer  corresponding 

to  probe  angle 

T(l)  -  Temperature  of  inside  plenum  (C) 

T(2)  -  Temperature  of  outside  plenum  (C) 

VP (1)  -  voltage  from  P(l)  pressure  transducer 

VP (2)  -  Voltage  from  p(2)  pressure  transducer 

VP (3)  -  Voltage  from  P(3)  pressure  transducer 

VP (4)  -  Voltage  from  p(4)  pressure  transducer 

VP (5)  -  Voltage  from  p(5)  pressure  transducer 

VP (6)  -  voltage  from  p(i)  -  p(3)  differential 

pressure  transducer 

VP (7)  -  voltage  from  tank  pressure  transducer 

VP (8)  -  Voltage  from  line  pressure  transducer 

VP (9)  -  voltage  from  outside  flow  orifice  plate 

differential  pressure  transducer 
VP(10)  -  Voltage  from  inside  flow  orifice  plate 

differential  pressure  transducer 


B)  OUTPUT  VARIABLES  (To  Files  #6  and  #7,  the  printer 
and  REDUCED  DATA  FILE  XXX,  respectively) 

1)  The  following  input  variables  are  used  in  the 
output  unchanged  (see  INPUT  VARIABLES  for  a 
description  of  each) 

VRUN ,  MONTH,  NDAY ,  NYR,  NTCONF,  N0RF(1) ,  N0RF(2), 

MS VANE,  PATM,  PLINEI,  NTAPS ,  N(l),  N(2),  N(3), 
NTIMEH(l) ,NTIMEM(1) ,  NTIMEH (2) ,  NTIMEM(2) ,  TRMS , 
TRME,  PTANKS,  PTANKE  &  NSRUNS 
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2)  Wall  Static  pressures  (averaged  over  all 
subruns)  and  their  locations 
PS4(J)  -  Static  pressure  on  4"  OD  centerbody  at 
location  "j"  (J  *  1  to  N(l)) 

ZS4 ( J)  -  Axial  location  of  static  pressure  PS4(J) 
PS6 (J)  -  Static  pressure  on  inside  wall  of  6" 

tube  at  location  "j"  (J  *  1  to  N(2)) 

TH6 (J)  -  Angular  location  of  static  pressure  PS6(J) 

ZS6  (J)  -  Axial  location  of  static  pressure  PS6(J) 

NPA  -  Number  of  averaged  static  pressure 
locations  on  6”  tube 

PS6A(J)  -  Average  static  pressure  on  inside  wall 

of  6"  tube  at  location  "j"  (J  *  1  to  NPA) 
ZS6A(J)  -  Axial  location  of  static  pressure  PS6A(j) 
PS8(J)  -  Static  pressure  on  8"  ID  outer  wall  at 
location  "j*  (j  ■  1  to  N ( 3 ) 3 
ZS8  (J )  -  Axial  location  of  static  pressure  PS8(J) 


PS8 (J) 


3)  For  each 
LTIMEH(I)  - 
LTIMEM(I)  - 
PNULL(I)  - 
PTANK(I)  - 
PLINE(I)  - 
FR1(I) 

FR2  (I) 

TAI(I) 
TAO(I) 
PTI(I) 
PTO(I) 
RADIAL (I)  - 
AXIAL ( I )  - 

ALPHA(I)  - 
BETA (I) 
VEL(I) 

UV(I) 

W(I) 

WV(I) 

PT(I) 

PS(I) 

MZ  (I) 

MR(I) 

FLOWID(I)  - 
FL0W2D (I)  - 
FLOWIM(I)  - 


subrun  "I" 

Hour  at  start  of  subrun 

Minute  at  start  of  subrun 

pressure  difference  p(l)  -  p(3)  (psi) 

Tank  pressure  (psig) 

Line  pressure  (psig) 

Inside  flow  orifice  plate  pressure  drop 
(inches  of  water) 

Outside  flow  orifice  plate  pressure  drop 
(inches  of  water) 

Inside  plenum  air  temperature  (C) 

Outside  plenum  air  temperature  (C) 

Inside  plenum  total  pressure  (psig) 

Outside  plenum  total  pressure  (psig) 

Radial  position  of  probe  head  (inches) 

Axial  location  of  probe  head  (inches) 

Flow  angle  alpha  (  ) 

Flow  angle  beta  (  ) 

Total  velocity  (ft/sec) 

Ratio  of  radial  to  total  velocity  (u/v) 
Ratio  of  tangential  to  total  velocity  (v/v) 
Ratio  of  axial  to  total  velocity  (w/v) 

Total  pressure  (psi) 

Static  pressure  (psi) 

Axial  momentum  flux  (psi) 

Tangential  momentum  flux  (psi) 

Mass  flow  rate  of  inside  flow  choked 
orifice  plate  (lbm/sec) 

Mass  flow  rate  of  outside  flow  choked 
orifice  plate  (lbm/sec) 

Mass  flow  rate  of  inside  flow  meter 


,s 
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orifice  plate  (lbra/sec) 

FL0W2M(I)  -  Mass  flow  rate  of  outside  flow  meter 
orifice  plate  (lbra/sec) 

VELAVG (I)  -  Average  axial  velocity  in  test  section 
( ft/sec) 

PQAVG(I)  -  Dynamic  pressure  corresponding  to  VELAVG (I) 
WI(I)  =  WZ(I)  *  RADIAL (I) 

RAD ( I )  =  RADIAL ( I } 

TI(I)  -  MR (I)  *  RADIAL (I) 

X1(I)  =  WV(I)  *  VEL(I)  *  RADIAL (I) 

PAI(I)  ■  PT (I)  *  XI(I) 

4)  Output  data  for  the  run  obtained  by  integration 
of  data  from  the  subruns 

W  -  Axial  momentum  at  axial  location  of  run  (lbf) 

TO  -  Angular  momentum  at  axial  location  of  run  (in-lbf) 
S  -  Swirl  number  at  axial  location  of  run 
WAVG  -  Average  axial  velocity  of  run  from  probe 
measurements  (ft/sec) 

PTM  -  Mass  average  total  pressure  at  axial  location 
of  run  (psi) 
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TABLE  F-l 


variables  used  in  program  "PLDATA" 

A.  INPUT  VARIABLES 

1.  From  File  ONAME  (INPXXX.DAT) 

NN  -  Number  of  input  data  files 

FNAME  -  Name  of  input  data  file  (PLTXXX.DAT) 

N1 ( I )  -  Number  of  data  files  for  plot 

1=1,  Static  pressure  on  4"  tube 
1-2,  Static  pressure  on  6"  tube 
1=3,  Static  pressure  on  8"  tube 
I  =  4f  Axial  velocity 
1=5,  Radial  velocity 
1=6,  Tangential  velocity 
1=7,  Total  velocity 
1=8,  Static  pressure 
1=9,  Total  velocity 
I  =  10,  Tangential  momentum 
I  *  11,  Axial  momentum 

N2  -  Number  of  last  plot  requested,  i 

PLIMIT(I)  -  Lower  limit  on  x  axis 

PLIMIT ( 2)  -  Upper  limit  on  x  axis 

PLIMIT(3)  -  Lower  limit  on  y  axis 

PLIMIT (4)  -  Upper  limit  on  y  axis 

ND ( 1 )  -  Number  of  divisions  on  x  axis 

ND (2)  -  Number  of  divisions  on  y  axis 

2.  From  File  FNAME  (PLTXXX.DAT) 

NSRUNS (II)  -  Number  of  subruns  in  file  FNAME 
N (II , 1)  -  Number  of  data  points  for  static 

pressure  on  4"  tube 

NPA(II)  -  Number  of  data  points  for  static 
pressure  on  6”  tube 

N (II ,3)  -  Number  of  data  points  for  static 

pressure  on  8"  tube 

ZS4  (II ,J)  -  Axial  location  on  4"  tube, 

J  =  1  to  N  ( 1 1 , 1 ) 

PS4(II,J)  -  Static  pressure  on  4"  tube, 

J  =  1  to  N  ( 1 1 , 1 ) 

ZS6A(II,J)  -  Axial  location  on  6"  tube, 

J  »  1  to  NPA(II) 

PS6A ( I I , J )  -  Static  pressure  on  6"  tube, 

J  »  1  to  NPA(II) 

ZS8  (II ,J)  -  Axial  location  on  8"  tube, 

J  =  1  to  N (II  ,3) 

PS8(II,J)  -  Static  pressure  on  8”  tube, 

J  =  1  to  N ( II ,  3) 


For  each  data  point  I,  I  =  1  to  NSRUNS(II) 


RADIAL (II , I) 
VEL ( 1 1 , 1 ) 
UV(II,I) 

W (II , I) 

WV ( 1 1  /  I ) 
IPOINT ( 1 1 , 1) 
PT(II,I) 
PS(II,I) 
MZ(II,I) 

MR (II /I) 
VELAVG (II , I) 
PQAVG (II/I) 


Radial  location 

Total  velocity 

Radial  velocity 

Tangential  velocity 

Axial  velocity 

probe  calibration  region 

Total  pressure 

Static  pressure 

Axial  momentum 

Tangential  momentum 

Test  section  average  axial  velocity 

Dynamic  pressure  of  VELAVG(II,I) 


B.  OUTPUT  VARIABLES  TO  SUBROUTINE  "PVPLOT" 


X (II ,K) 
Y(II,K) 

NA (II) 

IB 

ERR(II,K,1) 

ERR ( II ,K , 2) 

PLIMIT(J) 

ND  (I) 

NC 


-  value  of  X  for  data  point  K  in  run  II 

-  value  of  Y  for  data  point  K  in  run  II 

-  Number  of  data  points  in  each  run, 

II  =  1  to  NC 

-  Plot  format,  IB  »  1  to  ll 

-  X  lower  limit  for  data  point  K 
in  run  II 

-  x  upper  limit  for  data  point  K 
in  run  II 

-  Lower  and  upper  plotting  limits 

for  x  and  Y 

-  Number  of  divisions  on  x  and  y  axes 

-  Number  of  runs  to  be  plotted 
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APPENDIX  G 


DEVELOPMENT  OF  INTEGRAL  RELATIONSHIPS  FOR  AXISYMMETRIC 
FLOW  WITH  SWIRL  IN  AN  ANNULAR  CHANNEL  FORMED  BETWEEN 
CONCENTRIC  CONSTANT  AREA  CYLINDERS 

I.  INTRODUCTION 

The  annualr  channel  of  interest  and  nomenclature  for 
ax i symmetric  flow  are  shown  in  Figure  1.  Integral 
relationships  for  momentum  of  the  mean  velocities  in 
turbulent  flow  within  this  channel  are  developed  for  use 
in  analysis  of  experimental  data  obtained  for  this  flow 
configuration. 

II.  DEVELOPMENT  OF  INTEGRAL  MOMENTUM  EQUATIONS 

For  circular  jets  with  or  without  swirlf  the  mean  axial 
velocity  (w)  is  much  greater  than  the  mean  radialm 
velocity  (u)  and  gradients  in  the  radial  direction  are 
much  larger  than  those  in  the  axial  direction.  Becuase 
the  flow  of  Interest  is  for  a  nozzle  Reynolds  number 
greater  than  a  few  thousand,  the  viscous  stresses *are 
assumed  to  be  much  smaller  than  the  corresponding 
turbulent  shear  stresses.  Under  the  above  conditions,  the 
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time  averaged  Navier  stokes  equations  yield  the  following 
relationships  for  incompressible  flow: 
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where  the  prime  (•)  denotes  the  turbulent  fluctuation  and 
the  overbar  denotes  the  time  average  of  the  quantity. 


Because  a(v'w' ) /  3z  is  very  small  in  comparison  with 

3(u* v' )/  ar,  it  can  be  neglected  in  Equation  g-2  above. 
Neglecting  the  turbulent  normal  stresses  (u'2  ,  v'  and  w'2) 
and  denoting  the  turbulent  shear  stresses  as 


»  -  p  u'w' 
T§  -  -  p  u'v' 


(G-4) 

(G-5) 


then  Equations  G-l  through  G-3  become 
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3W  aw  1  3P  1  a 

u—  +w  —  - - + —  —  (rTz) 

3r  3z  p  3z  pr  3r 


Integration  of  Equation  g-6  with  respect  to  the  radius 
from  the  inner  radius  (r,  )  to  a  radius  (r)  in  the  flow 
yields 


Multiplying  Equation  G-7  by  pr*  and  integrating  with 
respect  to  the  radius  from  the  inner  radius  (r,  )  to  the 
outer  radius  (ra)  yields 

I  r° 

)  (G-10) 
1*1 

Multiplying  Equation  G-8  by  pr  and  integrating  with 
respect  to  the  radius  from  the  inner  radius  (  r,  )  to  the 
outer  (r0)  yields 


Equations  G-10  and  g-11  car  be  s*  ylified  further  by 


defining  the  terns  within  each  integral  as  follows 


jo rlwv  dr  »  j  Me 


r  dr 


(G-12) 


o 

w  *  y (p  - 


Pa  +  pw2)  r  dr 


r  dr 


(G-13) 


where  Po  is  a  reference  pressure 


Using  the  above  definitions  of  G8and  W,  Equations  C-10 
and  G-ll  sinplify  to 


dG« 

<r'T,> 

dz 

dw 

1 

dz 

Equations  G- 

14  and 

that  can  also  be 

(G-14) 


(G-15) 


(Reference  1)  between  the  change  in  the  axial  direction 
of  the  radially  integrated  angular  momentum  to  the 
tangential  shear  stress  at  the  walls  and  the  change  in 
the  axial  direction  of  the  radially  integrated  pressure 
plus  axial  momentum  (total  "stream  thrust”)  to  the  axial 
shear  stress  at  the  walls. 
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III.  SWIRL  NUMBER 


The  swirl  number  for  a  flow  is  the  ratio  of  the  angular 
momentum  to  the  total  "stream  thrust"  times  an 
appropriate  radius,  it  is  a  simensionless  number  that  is 
used  to  characterize  flows  containing  swirl.  The  exact 
form  of  the  equation  used  to  calculate  the  swirl  number 
varies  between  references  (see  References  2  and  3) .  The 
swirl  number  (S)  can  be  defined  as 

S  =  Ge/(^W)  (G-16 ) 

where  Ge  and  W  are  defined  by  Equations  G-12  and  G-13, 
respectively,  and  r0  is  the  outer  radius  of  the  annular 
channel.  As  defined  above,  the  swirl  number  is  constant 
in  the  axial  direction  when  the  shear  stresses  at  the 
wall  are  negligible. 


An  alternate  definition  of  the  swirl  number  (S)  can  be 
obtained  by  noting  that  Equation  G-13  can  be  integrated 
l^y  parts  to  give 


O 

w-  J 


(P  -  P9+  pw!)  r  dr 


w 


(p  -  Pe  )  rJ 


r0  r0 


+ 

*i  ri 


J P(w2- 


v2/2)  r  dr  (G-17) 


Defining  the  axial  momentum  flux,  Gz ,  as 


Gz 


f  P(w2- 


vJ/2)  r  dr  ■  j  M'r  dr  (G-18) 


then  Equation  G-17  becomes 


W  =  Gz  + 


(P  -  Pn  )rJ 


(G— 19 ) 


The  swirl  number  (S)  can  then  be  defined  as  follows 

Ge 

s  »  - - - -  (G-20) 

r.G.t  (r/2)  [<Pf(>-  PJ  lj  -  <Pf)-  P.  )  tf] 

where  Ge  and  Gz  are  defined  by  Equations  g-12  and  G-18, 
respectively,  p  and  p  are  the  wall  static  pressures  at 

■O  I 

r*r0  and  r=r,  ,  respectively,  p0  is  a  constant  reference 
pressure,  and  r,  and  r0  are  the  inner  radius  and  the 
outer  radius  of  the  annular  channel,  respectively. 


Equation  G-20  for  the  swirl  number  (S)  is  more  desirable 
to  work  with  than  Equation  G-16  because  the  intergrals 
for  determining  Gg  and  Gz  contain  only  the  axial 
velocity,  w,  and  the  tangential  velocity,  v.  whereas 
Equation  G-16  for  the  swirl  number  (S)  requires  that  the 
static  pressure  (p)  in  the  flow  be  integrated  between  the 
inner  and  outer  wall  for  the  calculation  of  W  (Equation 


G— 13 ) .  Since  present  measurement  methods  permit  a  more 
accurate  measurement  of  the  wall  static  pressure  than  the 
static  pressure  in  the  flow.  Equation  G-20  yields  a  more 
accurate  value  of  the  swirl  number  (s). 

Equation  G-20  for  the  swirl  number  (s)  can  be  simplified 
for  the  cases  of  a  swirling  stream  flowing  into  a  space 
of  uniform  pressure,  p0  ,  where  a  centerbody  exists 
(Figure  3-3)  and  where  a  centerbody  does  not  exist 
(Figures  3-1  and  3-2) .  The  swirl  number  for  the  case  when 
a  centerbody  does  exist  is  given  by 

Ge 

S  - -  (G-21) 

ro<3*  +  (^2)(PC-  Pr))r, 

where  r0  is  the  outer  radius  of  the  stream  at  its  exit,  r, 
is  the  radius  of  the  centerbody,  p0  is  the  ambient 
pressure,  and  pf|  is  the  wall  static  pressure  on  the 
centerbody.  The  swirl  number  for  the  case  when  no 
centerbody  exists  is  given  by 

S  -  G§ / ( re Gz  )  (G-22) 

where  re  is  the  outer  radius  of  the  stream  at  its  exit. 
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